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Abstract 

With over 1800 planets discovered outside of the Solar System in the past two decades, 
the field of exoplanetology has broadened our perspective on planetary systems. Research 
priorities are now moving from planet detection to planet characterization. In this context, 
transiting exoplanets—planets that cross in front of their star from our point-of-view—are 
of special interest due to the wealth of data made available by their orbital configuration. In 
this thesis, I introduce two methods, and their Markov chain Monte Carlo implementations, 
to gain new insights into the atmospheric and interior properties of exoplanets. 

The first method aims to map an exoplanet’s atmosphere based on the eclipse scanning— 
which is obtained while a planet is occulted by its host star. Ultimately temperature, compo¬ 
sition, and circulation patterns could be constrained in three-dimensions from these maps, 
a significant asset for informing atmospheric models. I introduce the basics of eclipse map¬ 
ping, its caveats (particularly, the correlation between the planet’s shape, brightness distri¬ 
bution, and four system parameters), and a framework to mitigate the caveats’ effects via 
global analyses including transits, phase curves, and radial velocity measurements. I use 
this method to create the first two-dimensional map and the first cloud map of an exoplanet 
for the hot-Jupiters HD 189733b and Kepler-7b, respectively. 

The second method, MassSpec, aims to determine transiting planet masses and atmo¬ 
spheric properties solely from transmission spectra, i.e. the starlight filtered by a planet’s 
atmosphere during transits. Determination of an exoplanet’s mass is key to understanding 
its basic properties, including its potential for supporting life. To date, mass constraints 
for exoplanets are predominantly based on radial velocity measurements, which are not 
suited for planets with low masses, large semi-major axes, or those orbiting faint or active 
stars. I demonstrate that a planet’s mass has to be accounted for by atmospheric retrieval 
methods to ensure unbiased estimates of atmospheric properties. Utilizing MassSpec, the 
James Webb Space Telescope (launch date: 2018) could determine the mass and atmo¬ 
spheric properties of half a dozen Earth-sized planets in their host’s habitable zones over 
its lifetime, which could lead to the first identification of a habitable exoplanet. 
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Professor of Physics 
Class of 1941 Professor 
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Chapter 1 


Introduction 


Planets outside of the Solar System, known as exoplanets, have been an objeet of inves¬ 
tigation for more than two centries (e.g., Jaeob[ [T855 ). Yet, the first exoplanet around a 
main-sequenee star was deteeted less than two deeades ago ( Mayor and Quelo"^ 19951. 
Sinee then, the field of exoplanetology has signifieantly and irreversibly broadened our 
perspeetive on planetary systems with more than 1800 planets found in over 1100 systems 
different from our^ We ean now plaee observational eonstraints on the planet popula¬ 
tion in our galax)|^ and eontemplate the results of our seareh for habitable world^ whieh 
reeently reaehed a climax with the first Earth-sized exoplanet found in a habitable zonej^ 
( [Quintana et al.[ 2014| ). A significant amount of theoretical, observational, and engineer¬ 
ing efforts is still required to characterize distant worlds, understand their climates, assess 
their habitability, and, hence, identify habitable worlds. While the next generation of ob¬ 
servatories is on its way, it is not clear yet what are the next scientific achievements future 
observatories will bring within our reach and, most importantly, how to reach them. 


'For an up-to-date list, refer to the Extrasolar Planets Encyclopaedia (http://exoplanet.eu 
Schneider et al. 2011|l and the Exoplanet Data Explorer (http: //exoplanets . org/, Wright et al. 


20111 . 


Wt least 100 billion planets orbit stars in the Milky Way (Cassan et al. 2012| l. As many as 40% of those 
planets are Earth-sized planets orbiting in the habitable zones of Sun-like stars and stars with later spectral 


type (Petigura et al. 2013 Swift et al. 
^ Discovery publications include 


20131. 


Vogt k al.|(2010l; 


hi; 


Wordsworth et al. 


Barclay et al. (20131; 


Anglada-Escude et al. (2013|l; Tuomi et al. (2013a 

“'The habitable zone is the region around a star within which planetary' 
mospheric pressure can support liquid water at their surfaces. 


(|201l|l; 


Borucki et al. 


Borucki et al. 


(2012i; 


(20131. 


objects with sufficient at- 
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1.1 First Steps towards Worlds beyond the Solar System 


1.1.1 Detections 


Planet detections rely largely on indirect detection methods because most exoplanets cannot 
be directly detected with current observatories. Exoplanets are too faint and/or too close 
to their stars to be imaged. Over the first decade following [Mayor and Queloz (1995 )’s 
discovery, most of the exoplanets found were hot Jupiter^ (see Figure flO^ and e.g., Marcy 
et al.[ [2005 [ and references therein). The large mass and the proximity of hot Jupiters 


to their hosts favored their indirect detection by the radial velocity (RV) methocQ which 
searches for the Doppler shift of a star’s spectrum as it orbits the planet-star common center 
of mass. The RV method contributed ~ 30% of the planet discoveries to date. The most 
prolific planet discovery technique is the transit metho(0(~ 60% of the planets discovered 
to date), which searches for the brightness drop of a star as a planet transits, i.e., passes 
in front of it. One mission, the Kepler mission, has changed the planet-detection game 


as it detected 961 confirmed exoplanets and more than 2900 planet candidates (Borucki 
et al. 2010 2011a|b ; Batalha et al.[|2013 ; Burke et al. 2014| via monitoring over 145,000 


main-sequence stars for over 3 years. 

The transit and the RV methods are complementary in two ways: (1) in terms of detec¬ 
tion biases and (2) for characterization purposes. (1) The sensitivity of the transit method 
depends mainly on the planet-to-star area ratio, while the sensitivity of the RV method de¬ 
pends on the planet-to-star mass ratio and distance. In other words, the detection of smaller 
planets is easier with the transit method. (2) Transit and RV measurements constrain the 


planet-to-star area and (minimum) mass ratios, respectively (e.g., Winn 20101. Hence, the 
observation of a planet with both techniques enables the first step towards its characteriza¬ 


tion via its bulk density. A planet’s density can yield its bulk composition (Seager et al. 


1 2007 [[Fortney et al.[|2007[ ), whether it is a gas giant or a rocky exoplanet potentially suitable 
for life, as we know it. 

^Hot Jupiters are Jupiter-sized planets heated by their close stars—within 0.1 astronomical unit (AU). 

^The leading initiatives for planet discovery using RV are the Geneva Extrasolar Planet Search (ELODIE, 
CORALIE, HARPS, and HARPS-N spectrographs) and the Lick-Carnegie Exoplanet Survey (HIRES spec¬ 
trometer). The first initiative led to the first planet discovery around a main sequence star, the second con¬ 
tributed over 70% of the discoveries performed before 2010. 

^Leading initiatives for planet discovery using transits include HATNet Project, CoRoT, and WASP. 


18 




































10^ 


10'^ 


_® 10 


w 

? 10^ 


o 

c 

n3 0 

E 10 


10" 


10" 



. 

.* ‘ M 


, 



• • 
• h , * • 

' • 



. • s 

•• 



• 


• 

U N 



LU 

> 



Ma 



Me 


10 ' 


10“ 


U N 


Ma 


10 "" 10 “ 10 " 10 ’ 
Orbital period [days] 


10 “ 


10" 


Me 

-2 


10 10 10 
Orbital period [days] 


10” 


Figure 1-1: Planets detected as of 2005. The planets of the Solar System are shown in red 
and the exoplanets as black dots. The black line on the period-mass plot represents the 
brown-dwarf transition, the mass around which a planet-sized body can fuse deuterium. 




10 10 10 10 
Orbital period [days] 

Figure 1-2: Planets with constrained bulk density as of May 2014. The exoplanets with 
constrained bulk density are shown as green dots and the other exoplanets known as black 
dots. The panels show the same data type as Figure [TT| 

1.1.2 Characterization 

The community determined the bulk composition of ~ 20% of the exoplanets known (Fig- 
1-2|). Simultaneously, complementary techniques were developed to further the char- 


ure 


acterization of exoplanets. As an example, one of these techniques lead the detection of 
atmospheric components—atomic and/or molecular—for about half a percent of the de- 
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A^: High 


A2: Low 

atmospheric 


atmospheric 

absorption 


absorption 


tected exoplanets (i.e., 10). Now, how ean we leam so mueh about a planet when its whole 
planetary system appears to us as a single bright point in the sky? 

Transiting exoplanets are, among the exoplanets deteeted to date, the ones that ean 
be most extensively eharaeterized with eurrent teehnologjj^ Beyond derivation of the ex¬ 
oplanet orbital inelination and density (e.g., |Winn[ |2010[ ), transiting exoplanets are key 
objeets beeause their atmospheres are observationally aeeessible through transit transmis¬ 
sion and oeeultation emission speetrophotometry (see e.g., [Seager and Deming[ [SOTOj and 
referenees therein) and phase eurve measurements (Figure fT^ . 

Transmission Spectroscopy studies the wavelength- 
dependent flux drop during a planetary tran¬ 
sit ( [Seager and Sasselov[ |2000| ). While a 
planet transits, part of its star’s light is going 
through its atmosphere (Figure fF^j ). Henee, this 
light ean show absorption features from the at- 
mospherie eomponents, whieh affeet the atmo¬ 
sphere’s opaeity and its wavelength dependenee. 

At a wavelength with high atmospherie absorp¬ 
tion, Ai, the transit depth is slightly deeper than 
at a wavelength with lower atmospherie absorp¬ 
tion, A 2 —due to the more extended opaque at¬ 
mospherie annulus. Also, the planet appears 
larger as a relative flux-drop, ^(A), is assoei- 
ated with an effeetive planetary radius, Rp{X) = 

Figure 1-3: Basies of transmission speetroseopy. Transit-depth variations, ^(A), indueed 
by the wavelength-dependent opaeity of a transiting planet’s atmosphere. The stellar disk 
and the planet are not resolved; the flux variation of a point souree is observed. 




^Non-transiting exoplanets can also be characterized, but to another degree. In particular, the cross¬ 
correlation technique successfully implemented by Snellen et al. ( 2010) l allows for the detection of atmo¬ 
spheric components and the measurement of atmospheric winds and orbital inclination, which yield the planet 
mass when combined with RV—see, e.g., applications to the hot Jupiters r Bootis b (Brogi et aT] [20 1 2[ [Rodlerj 


et al. 


2012 1 , HD 179949 b (Brogi et al. 20141, and the young /3 Pictoris b (Snellen et al. 


20141. 
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Figure 1-4: Basics about transiting exoplanets. During a transit, the flux observed drops 
because the planet blocks a fraction of the stellar disk—whence, the flux drop relates to 
the stellar brightness and the planet-to-star area ratio. Then the flux rises as the planets 
dayside comes into view—most exoplanets detected so far always face their host with the 
same hemisphere (dayside), they are tidally-locked as the Moon is to Earth. During an 
occultation, the flux drops relates to the emerging planetary flux. Then the flux drops as the 
planets nightside comes into view. Figure from I Winnl (120 1 01) . 
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Emission Spectroscopy studies the wavelength-dependent flux drop during a plane¬ 
tary oeeultation. Observations of the oeeultation depth at different wavelengths eonstrain 
the hemisphere-integrated emission speetrum of an exoplanet. The emission speetrum of 
an exoplanet eonstrains its brightness temperature, but also its atmospherie eomposition. 
As atmospherie eomponents affeet the wavelength dependenee of the atmospherie optieal 
depth|^ different wavelengths probe different atmospheric layers with different brightness 
temperatures. In other words, the emission spectrum of an exoplanet translates primarily 
into constraints on the temperature and composition structure of its atmosphere. 

Phase curve measurements relate to the observations of a planetary system’s flux vari¬ 
ations with the planetary phase, which provide information about the planets’ brightness 
distribution. For most exoplanets detected, the rotation rate is synchronized with the or¬ 
bital perioc^ because close-in exoplanets are easier to detect and subject to strong tidal 
friction. Hence those planets always show the same face to their stars (dayside), similarly 
to the Moon with the Earth. As a result, along their orbits different planetary hemispheres 
are visible, which modulate the planetary flux contribution (Figure 1-4| ). Hence, monitoring 
the system flux modulation over an orbital period for synchronized planets can be translated 
into constraints on their longitudinal brightness distribution. 


Overview of Exoplanet Characterization Achievements 


The Hubble Space Telescope (HST) and the Spitzer Space Telescope {Spitzer, Werner et al. 


20041 —HSHs cousin in the infrared—have played major roles in furthering the degree to 
which exoplanets can be characterized. HST yielded the first detection of an exoplanet at¬ 


mosphere by revealing its sodium absorption feature in transmission (Charbonneau et al. 


2002 1. Spitzer yielded the first (broadband) emission spectrum of an exoplanet (Deming 
et al.[ 2005 [ [Charbonneau et alT| 2005), which indicated a temperature inversion in the tar¬ 


get’s atmosphere ( [Knutson et al.[ |2008| ). Spitzer also led to the first map of an exoplanet 
based on its phase curve measurement ( [Knutson et al.[[2O07] ). HSTs observations revealed 


®The optical depth expresses the quantity of light removed from a radiation beam by extinction as it 
passes through a medium. In particular, the ratio of incoming intensity to emerging intensity is where 

T is the optical depth. 

'°Close-in rocky planets can be locked in configurations different from 1:1 (e.g., 

Makarov and Berghea 2014[l. 


Makarov et al. 


2012 
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for the first time the atmospherie eseape from an exoplanet (|Leeavelier Des Etangs et aF 


20101 and the interaetion between an exoplanet’s atmosphere and stellar aetivity (Leeave 


her des Etangs et am2012| ). Reeently, the Wide Field Camera 3 (WFC3) on HST provided 
the most exquisite exoplanetary transmission speetra between 1 and l.V/rm revealing the 
1.4/im absorption feature of water for a few hot Jupiters ( Doming et al.[ 2013[ Huitson 
et al4 2013[ Mandell et al. 20131 and a flat signal within the measurement unoertainties for 


others (Berta et al. 2012 Knutson et al. 2014[ Ranjan et al. 2014; Kreidberg et al. 20141. 


1.2 Context 


By the time this thesis started, the quality of exoplanetary datasets was in rare oases suffl- 
oient to yield longitudinal maps of exoplanets and deteotions of atmospherie eonstituents 
using transmission speetroseopy. The perspeetives of new observatorie^ and unpreee- 
dentedly informative datasets resonated with the field’s goals to eharaeterize exoplanets, 
probe their atmospheres, and assess their habitability. Yet, it was not olear what are the 
new insights into the atmospherie and interior properties of exoplanets that will soon be 
within our reaoh and how to reach them. 

As an example, the assessment of a planet’s habitability neeessitates the determina¬ 
tion its surfaee eonditions—i.e., temperature, pressure, eomposition, and gravity. In most 
eases, the atmospherie layers at the surfaee level eannot be probed. Henee, atmospherie 
models are needed to extrapolate the atmospherie eonditions from the layers probed by 
transmission and/or emission speetroseopy to deeper layers, possibly ineluding the sur¬ 
faee. Transmission speetroseopy eonstrains the atmospherie eonditions at relatively low 
pressures (typieally below 0.1 bar) and over the planetary limbs—i.e., not a planet’s 3D 
elimate. Similarly, emission speetroseopy eonstrains the hemisphere-integrated eonditions 


^'including the James Webb Space Telescope (JWST; launch date 2018, |^ampin 20101, the Exoplanet 
Characterisation Observatory (£’C/!(9|^ the European Extremely Large Telescope (E-ELT), the Thirty Meter 
Telescope (TMT), and the Giant Magellan Telescope (GMT). 

^^EChO was a M3 mission candidate of the European Space Agency (Tinetti et al. 20121. While EChO 
was not funded, such a mission could still revolutionize our understanding of exoplanet atmospheres. For 
that reason, we will emphasize in this thesis the benefit of an EChO-class mission in the context of exo¬ 
planet atmosphere characterization. For brevity, we henceforth refer to such a future FC/zO-class mission as 
“ECha\ 
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of a planet’s atmosphere—i.e., degenerate constraints on atmospheric models. In other 
words, 3D maps of exoplanetary atmospheres would be a valuable asset to extrapolate the 
temperature, pressure, and composition of unprobed layers. In addition 3D maps would 
provide feedback to atmospheric models and may reveal atmospheric processes or forcing 
factors yet unmodeled. The first method introduced in this thesis aims to address that need. 
The determination of a planet’s surface gravity is also required to assess its habitability. 
The surface gravity is derived from the planet mass, which is traditionally determined by 
the RV method. Yet, the RV method is mainly effective for massive planets around rela¬ 
tively bright and quiet stars. Hence the masses of small—e.g.. Earth-sized—planets could 
remain out of reach, preventing from their in-depth characterization, including the assess¬ 
ment of their habitability. The second method introduced in this thesis aims to determine 
the planetary mass in an alternative way to the RV method. 
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1.3 Thesis Overview 


The continual increase in data quality of exoplanet observations drives the development of 
new methods to further characterize distant worlds. The main objective of this thesis is to 
investigate theoretically and practically the question ""What are the new insights into the 
atmospheric and interior properties of exoplanets within reach and how to reach them?" 

In Chapterj^ I introduce how transiting exoplanets can be mapped while their stars can¬ 
not be resolved. I present the degeneracies intrinsic to exoplanet mapping and the frame¬ 
work proposed to mitigate those degeneracies while performing consistent mapping. I 
apply my mapping method to the hot Jupiters HD 189733b and Kepler-7 b, respectively in 
the infrared and in the visible. 

In Chapterj^ I introduce MassSpec, a method to determine the mass of transiting ex¬ 
oplanets based solely on their transmission spectra. MassSpec simultaneously and self- 
consistently constrains the mass and the atmospheric properties of an exoplanet, and al¬ 
lows mass measurements for transiting planets for which the radial velocity method fails. I 
demonstrate MassSpec's feasibility theoretically and highlight the fact that a planet’s mass 
has to be accounted for by atmospheric retrieval methods to ensure retrieval quality. I 
also show the good agreement between the mass MassSpec retrieves for HD 189733b from 
transmission spectroscopy with that from RV measurements. 

In Chapter]^ I investigate the prospects of the methods introduced in Chapters]^ and|^ 
I introduce short-term goals such as mapping for the first time an exoplanet’s atmosphere 
in 3D with the Spitzer Space Telescope and weighing the hottest exoplanet known with the 
Hubble Space Telescope. I discuss long-term goals for the next decade such as obtaining 
time-dependent 3D maps of exoplanetary atmospheres and assessing the habitability of 
Earth-sized planets with the James Webb Space Telescope and/or an ECfiD-class mission. 

In Chapter]^ I summarize the thesis’ results and conclusions. Appendix[A| introduces 
the basics of rational functions. In Appendix]^ I present the first step towards a general¬ 
ization of the transmission spectrum equations. I use these new equations in Appendices]^ 
andj^ to introduce a simple validation test and an efficient implementation procedure for 
numerical transmission spectrum models, respectively. 
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Chapter 2 


Mapping Planets’ Atmospheres While 


Their Stars Cannot Be Resolved 


Until recently, observations of extra-solar planetary systems primarily constrained their 
orbital configurations and, to a lesser extent, some planetary atmospheres—accessible 
through transmission and emission spectrophotometry (see e.g., Seager and Demin^ 2010 


and references therein). Although significant theoretical developments have been achieved 
in modeling exoplanet atmospheres by combining hydrodynamic flow with thermal forcing 
Showman et al.| ( |2009| ); |Rauscher and Menou| ( |2010| ); [Dobbs -Dixon et al.| ( |20T0| ); |Rauscher 


and Menou (20121 and/or with ohmic dissipation (e.g., Batygin et al.[|20lT Heng 2012 


Menou 20121, other potential forcing factors are yet to be modele(^^] —or observed. The 


observation of specific spatial features within an exoplanet atmosphere, such as hot spots 
or cold vortices, is essential for constraining its temperature, composition, and circulation 
patterns. Therefore, mapping exoplanet atmospheres is key for gaining further insights into 
their physics. 

Eclipses have proven to be powerful tools for “spatially resolving” distant objects, in¬ 
cluding binary stars (e.g., [Warner et al.[|1971[ ) and accretion disks (e.g., |Horne[|1985| ). Pre¬ 
vious theoretical studies investigated the potential of occultations in order to disentangle 

( | 2012 | , 


*Work published in Astronomy and Astrophysics, Volume 548, A 128 (2012), see 


de Wit et al. 


and in The Astrophysical Journal Letters, Volume 776, Issue 2, L25 (2013), see Demory et al. (2013 1 . 

'As an example, the magnetic star-planet interactions can also have a signihcant role in this matter. 
Nevertheless, magnetic interactions have to date been only observed at the stellar surface, in the form of 
chromospheric hot spots rotating synchronously with the companions (e.g., [Shkolnik et al.| [2005 1 [Lanzaj 


20091. 
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exoplanetary atmospherie eireulation regimes (e.g., Williams et al. 2006 Rauseher et al. 


2007b I. 


In this Chapter, we introduee how transiting exoplanets ean be mapped while their 
stars eannot, we present the degeneraeies intrinsie to exoplanet mapping and the frame¬ 
work proposed to mitigate those degeneraeies and perform mapping eonsistently. We also 
present two applieations of our method in this Chapter: first, the 2D map of the hot-Jupiter 
HD 189733b (based on de Wit et al. 2012[ ); seeond, the ID map of the hot-Jupiter Kepler-7 


b (based on the author’s eontribution to Demory et al. 20131. 


2.1 Concept 


The monoehromatie light eurve of a transiting and oeeulted exoplanet with a non-zero im- 
paet parameter ean enable a 2D surfaee brightness map of its day side. Sueh an exoplanet 
is seanned through several proeesses that provide flux measurements of different plane¬ 
tary sliees, whieh ean be reeombined into a map (see Figure 2-1) . First, the exoplanet 
is gradually masked/revealed by its host star during oeeultation ingress/egress (“eelipse 
seanning”). Seeondly, the exoplanet rotation provides its phase-dependent hemisphere- 
integrated flux—i.e., its phase eurve (see Seetion T.1.2| ). A planet’s phase eurve eonstrains 
its brightness distribution (hereafter, BD) in longitudinal shoes—as long as the exoplanet 
spin is olose to the projeetion plane, e.g., for a transiting and synohronized exoplanet. The 
three seanning proeesses (ingress, egress, and phase eurve) provide eomplementary pieoes 
of information that eonstrain the target’s BD over a speoifio “grid” (see the oomponent 
labeled ""combined” in Figure|27T]) beoause they soan the planet over different direotions. 


2.2 Degeneracies 

Eelipse mapping is subjeot to two levels of degeneraoy. 
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Figure 2-1: Schematic description of the different scanning processes observable for an 
occulted exoplanet. The green dotted lines indicate the scanning processes during its oc- 
cultation ingress/egress. The red dashed lines indicate the scanning process that results 
from the exoplanet rotation. The component labeled “combined” shows the specific grid 
generated by these three scanning processes. 


2.2.1 Eclipse Mapping: an Intrinsically-Underdetermined Problem... 


First, eclipse mapping is an underdetermined problem. The reason is that, conceptually, 
the map is composed of “iV^”-grid cells constrained by the measured fluxes in 2N slices 
during the target’s ingress and egress (compare the component labeled ""combined” to the 
components labeled ""ingress” and ""egress” in Figure [?T] ). Therefore, different brightness 
maps can be solutions to the same observations, hence the intrinsic degeneracy. To account 
adequately for this intrinsic degeneracy, mapping methods have to be able investigate the 
different type of maps that can fit a given dataset. For that purpose, the analysis method we 


introduce in Section 2.3.1 use multiple types of brightness models to reveal different types 


of brightness patterns. 


2.2.2 ... with Multiple Non-Intrinsic Contributing Factors 

Secondly, eclipse mapping is based on the shape of a target’s occultation ingress/egress. 
However, the shape of occulation ingress/egress is not solely affected by its brightness 
distribution. To propose a consistent methodology for mapping exoplanets, it is important 
to investigate the possible contributors to the shape of an occulation ingress/egress. Those 
possible contributors are: 
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• the planet’s projected shape at conjunctions, 

• the planet’s brightness distribution, 

• four of the planetary system parameter^ 

- the planet’s orbital eccentricity, e, 

- the planet’s periapsis argument, u, 

- the planet’s impact parameter, b = a/R^ cosi (where a is the exoplanet semi¬ 
major axis and i the orbital inclination), 

- the stellar density, p*. 

We present in Figure [2^ a schematic description of the dependency of the occultation 
shape on the planet’s shape (yellow) and brightness distribution (red). Both synthetic sce¬ 
narios show specific deviations from the occultation photometry of uniformly bright disk 
(black curve) in the occultation ingress/egress (see Figure [2^ bottom panel). 

We present in Figure [23] the dependency of the occultation shape on e, u, b, and p*. The 
top panel of Figure [2^ presents a simulated occultation using the system parameters of the 
hot Jupiter HD 189733b (black curve) and the simulated occultations when perturbing each 
of the system parameters. Although the effects may be considered as of second order—the 
deviations from the reference simulation are barely visible on the top panel, they actually 
corresponds to typical deviations expected in ingress/egress in order to perform mapping 
(compare with the bottom-right panel of Figure [27T0l ). 

Each of these deviations can be explained by one of the following effects, e and u 
enable both to shift the occultation by 

2P 

At —ecoscj, (2.1) 

71 

^In each set of parameters sufficient to model the planetary system orbits, four of them can affect the 
occulation shape. The four presented here relates to the parameter set introduced in Equation|2.3| 
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Deviation 


Figure 2-2: Schematic description of the effect of the shape or the brightness distribution of 
an exoplanet on its occultation shape. The red curve indicates the occultation photometry 
for a non-uniformly bright disk (hot spot in red). The yellow curve indicates the occultation 
photometry for an oblate exoplanet (yellow ellipse). Both synthetic scenarios show specific 
deviations from the occultation photometry of uniformly bright disk (black curve) in the 
occultation ingress/egress. 




Phase 


Figure 2-3: Effects of system parameters on an occultation ingress/egress shape. Top 
panel: Simulated occultation using the system parameters of the hot Jupiter HD 189733b 
(black curve) and the simulated occultations when perturbing system parameters (6 in blue, 
y/e cos oj in green, y/e sin oj in red, and p* in cyan). Bottom panel: Deviations in occultation 
ingress/egress from the reference model. 
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Figure 2-4: Schematic description of the effect of ecosce on an occultation shape. The 
green curve indicates the occultation photometry for a decreased e cos oj, which translates 
into a decreased timing-offset of the occultation (At). 
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Figure 2-5: Schematic description of the effect of esince on an occultation shape. The red 
curve indicates the occultation photometry for an increased e since, which translates into an 
increased occultation duration Woe —in particular an increased ratio occultation-to-transit 
durations 

Wtr 


32 






























Figure 2-6: Schematic description of the effect of the impact parameter (b) on an occul- 
tation’s shape. The red curve indicates the occultation photometry for a decreased b. A 
smaller impact parameter yields a longer eclipse duration and less gradual ingress/egress. 



*■ 


► t 


Figure 2-1 \ Schematic description of the effect of the stellar radius (i?*) on an occultation 
shape—representative here of the effect of the stellar density as the stellar mass is not 
independently constrained by a system light curve. The red curve indicates the occultation 
photometry for an increased i?*. A larger stellar radius yields a longer eclipse duration and 
more gradual ingress/egress. 
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(hereafter “uniform time offset’]^ and to ehange its duration 




Wf 


1 + 2e since, 


( 2 . 2 ) 


tr 


(Figures [2^ and [23j respeetively). Henee inereasing ecosce primarily leads to an earlier 
oeeultation (negative deviation in ingress and positive deviation in egress) and inereasing 
e sin ce, to an inerease of the oeeultation duration (negative deviation in ingress and nega¬ 
tive deviation in egress), b affeets the orbital inelination, henee, both the time of first and 
last eontaets and the oeeultation duration (Figure [2^ ). Deereasing b primarily leads to a 
negative deviation in ingress and positive deviation in egress. Finally, p* affeets direetly 
the stellar radius and henee the time of first and last eontaet (Figure [?7] ). Deereasing p* 
primarily leads to a positive deviation in ingress and positive deviation in egress. Note 
that also b and p* inerease the oeeultation duration when deereased, their effeets on the 
ingress/egress are different. While the ingress and egress slopes are steeper than the refer- 
enee ease (blaek) when deereasing b, their are more gradual when deereasing p*. Henee, 
eaeh of the four system parameters has a unique effeet on the oeeultation shape. 


2.2.3 Steps Towards Mitigating Degeneracies 

The various shapes of the deviations and their symmetry/asymmetry reveal that eombi- 
nations of perturbations of the system parameters ean adequately mimie the effeet of a 
non-uniformly bright exoplanet (see Figure [2^ . In other words, e, u, p*, and b are eorre- 
lated with the projeeted shape and brightness distribution of an exoplanet in the framework 
of oeeulation dataset. 


Avoid Unassessed Assumptions 


A first step toward mitigating the effeet of the degeneraeies is to avoid assumptions that 
relate to one of the eorrelated parameters mentioned above. An assumption relative to one 
of those parameters would affeet direetly the estimates of the other parameters, possible 


^The uniform time offset measures the time lag between the observed secondary eclipse and that predicted 
by a planet with spatially uniform emission Williams et al. (2006 1 . 
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biasing these estimates (see applieation in Seetion [2.4.3[ ). In other words, assumptions sueh 
as “the target’s orbit is eireularized” (i.e., e = 0) or “the target is a sphere” will prevent 
analysis methods from the proper exploration of the parameter spaee. We highlight here 
that assuming a planet to be uniformly-bright while fitting high-SNR oeeultations is also 
inadequate. For that reason, it is erueial to make a parsimonious use of assumptions unless 
they ean be assessed. 


Perform Global Analysis 


Assumptions regarding the planet’s shape, brightness, or system parameters ean be assessed 
using eomplementary data to oeeultation light eurves. Note that eomplementary oeeultation 
data obtained in different speetral bands ean also help in disentangling the eontributions of 
the system parameters to the oeeultation shape—as these are wavelength-independent. 

Transit Light Curves provide eomplementary eonstraints on the system parameters 
(ineluding b, a/c since, and p*) and on the projeeted shape of a eompanion at eonjunetion 


Barnes et al. (20091; Carter and Winn (2010). Furthermore, the signifieanee of shape- 
indueed deviations in transit is about one order of magnituded larger than the shape-indueed 
deviations in ooeultatioii|^ Note that eonstraining an exoplanet shape requires a priori 
knowledge on the host star limb-darkening (e.g.. Claret and Bloem^ 201 1[ ) and on the 
transit timing to avoid overfitting its possible signature in the transit ingress/egress—the 
transit timing is independently eonstrained by RV measurements. 

Phase Curves provide eomplementary eonstraints on a eompanion brightness and shape. 
Henee phase eurves are benefieial to mitigate both levels of degeneraey: ( 1 ) it provides 
eomplementary eonstraints on the planet’s brightness distribution (different seanning di- 
reetion, Figure |2^ and ( 2 ) it eonstrains the planet’s shape. A planet’s shape affeet its 
phase eurve in a way similar to a tidally distorted star induees ellipsoidal light variations 
( [Russell and MerriH] 1952; Kopak 1959). The eontributions of the shape and the BD of 
an exoplanet to its phase eurve may be disentangle beeause of their different periods (see 
e.g., Faigler and Mazeh] 20111. As an example, for a synehronized exoplanet, the shape- 


'^The effect ratio of shape-induced signals in oeeultation to in transit is proportional to the planetary 
hemisphere-averaged relative brightness < Ip >/</*> in the spectral band covered 
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induced modulation has a period of P/2—twice the same projection an orbital period, 
while the brightness-induced modulation as a period of P. 

RV Measurements provide complementary constraints on the orbital configuration, 
including on e and u (see Section [Z3.4[ ). 


Coping with Multiple Solutions 

Global analysis of a complementary dataset is key to constrain consistently the possible 
contributing factors to the shape of an occultation. However, the complex multidimensional 
correlation of the problem may produce multiple solutions—their disentanglement depend¬ 
ing on the data SNR. Additional arguments could then be used a posteriori to disentangle 
further between the different solutions proposed by the retrieval method. For example, con¬ 
sideration from atmospheric physics or statistics (e.g., the Bayesian Information Criteria, 


Gelman et al.l 20041 should be investigated. 


2.3 Methodology 

We present in this section the retrieval method developed to map a transiting planet, details 
regarding the data reduction are provided in Section [2At The method can also be used to 
analysize traditionnal datasets such as RV measurements and/or individual eclipses. 


2.3.1 Analysis Method 


We use an adaptive Markov Chain Monte Carlo (MCMC; see e.g. Gregory 2005 Ford 


20061 algorithm. MCMC is a Bayesian inference method based on stochastic simulations 


that sample the posterior probability distribution (PPD) of adjusted parameters (x) for a 
given fitting model, /sjm(x). More specifically, our implementation uses Keplerian orbit^ 


and models photometry using the light-curve models introduced in Section 2.3.2 and RV 


data using the model introduced in Section 2.3.4 In addition, each photometric time-serie 


is corrected for systematics (see Section 2.3.3) 


’Our orbit model also includes the host-star’s movement around the system center of gravity 
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Basics about MCMC Algorithms MCMC algorithms generate a ehain in the parameter 
spaee based on a sueeession of stoehastie jumps submitted to rules that ensure eonvergenee 
of the ehain towards the PPD. The basies of a ehain generation ean be deseribed by the 


following steps (Ford 20051: 


1. Initialize the ehain with some x„ (n = 0). 


2. Evaluate fsimi'^n) and the assoeiated ehi-squared (x^(x„)). 

3. Generate a trial state x' aeeording to a transition kernej^g (x'|x„) (e.g., a Gaussian 
funetion eentered on x„). 


4. Evaluate /sim(x') and the assoeiated x^(xO- 

5. Estimate the probability of transition from x^ to x' cx exp (-i (x^(x') - X^(Xn)))- 

6. Draw a random number from a uniform distribution between 0 and 1. 

7. If this number is below the probability of transition x„+i = x', else x„_|_i = x„. 

8. n = n + 1 and baek to step 3 unless n = numit- 


In ease of a priori knowledge eoneeming the distribution of parameters, the expression- 

whose evaluation eonditions the aeeeptanee of the attempted jump at step 5—is modified 
to inelude this a priori information. As an example, if Mgtar is known a priori to be of 
M* 0 ± (Tm* under normal distribution—i.e., N(M*^o, then the ehi-squared funetion 
ineludes this a priori knowledge in the form of the following additional term: I 


“Conventional” Part of our Parameter Set 


In order to model “eonventionaf flight eurves, the required parameter set (x) has to eon- 
strain the area and brightness ratios of the eompanions as well as their orbits. Several 


ehoiee of parameter sets are possible for modeling system light eurves, but it is neeessary 


®Refer to e.g., Gregory ( 2005 1 and reference therein for further details on the transition-kernel properties 
required to ensure the good behavior of the MCMC chain. 

^“Conventional” refers here to models or methods that assume uniformly-bright planets. 
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to aim for one with a mitigated eorrelation between the parameters to enhance the method 
efficiency (see discussion in Ford[ 2006). The conventional part of our implementation uses 


j^2 < / > 

6, Tt, P,^/ecos{uJ), ^/esm{uJ),p^, M*, / 

Ri < h> 


(2.3) 


where Rp/Rl is the planet-to-star area ratio, b is the impact parameter (= a/i?*cosz 
where a is the exoplanet semi-major axis and i is the orbital inclination), Tt is the time 
of minimum light (i.e., the transit center), P is the orbital period, e and u are the eccen¬ 
tricity and the periapsis argument, p* is the stellar densit)|^ M* is the stellar mass, and 
< Jp > / < /* > planetary hemisphere-averaged relative brightness. 


We used the following jump parameters 


i?; 


b, Tt, P, ^/ecos{u), ^/esm{u), and p* 


and assumed a uniform prior distribution for all these jump parameters. Because the stellar 
mass is not constrained by the targeted dataset, the method uses a Gaussian prior from the 
literature. 


Specificity of the Linear Parameters 

< Jp >/</*> is not part of our set of jump parameters because it is a linear parameter 
of our model. Hence, it can be determined by linear least squares minimization at each step 
of the MCMC. For this purpose, we employed the Singular Value Decomposition (SVD) 
method ( Press et aL||1992 ). Distinguishing between linear and non-linear parameters en¬ 
ables to increase significantly the efficiency of an MCMC implementation by reducing the 
dimension of the space to probe—i.e., the number of jump parameters. Therefore, the 
linear parameters of our model are not considered as jump parameters. 

We found that in order to use consistently the SVD method within such a stochastic 
framework, it is optimal to perturb the linear-coefficient estimates using the covariance 
matrix derived by the SVD. Perturbing the estimates prevents the SVD from mitigating 
systematically the increase due to inadequate jumps in the parameter space via the use 


^We choose to use the stellar density, p*, instead of the transit duration (from the first to last contact), W, 


as jump parameter; because it relates directly to the orbital parameters (Seager and Mallen-Ornelas 2003 i. 

®Jump parameters are the model parameters that are randomly perturbed at each step of the MCMC 
method. 
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of the best estimates of the linear parameters. 


2.3.2 Light-Curve Models 

Stellar Contribution: Transit Model 


We model the stellar contribution assuming a quadratic limb-darkening law for the star as 
Mandel and Agol|(2002). We draw the limb-darkening coefficients from the theoretical 


m 


tables of Claret and Bloemen (2011 1 based on the host-star spectroscopic parameters— 
effective temperature, Tg//, surface gravity, log ( 7 , and metallicity, We add this a 

priori knowledge as a Bayesian penalty to our merit function, using as additional “conven¬ 
tional” jump parameters the combinations Ci = 2ui + U 2 and C 2 = — 2 ^ 2 , as described 

inlGillon et al.|(|2010a]). 


Planetary Contribution 


We model the planetary contribution by performing a numerical integration of the observed 
exoplanet flux. Model approximations include: ignoring the time variability of the target 
atmosphere (in line with atmospheric models, e.g.. Cooper and Showmanj 20051—because 
it is currently necessary to combine multiple observations to obtain a sufficient signal-to- 
noise ratio (SNR); ignoring the planet limb darkening; and assuming the target’s rotational 
period to be synchronized with its orbital period (synchronization occurs over 10 ® yr for 
hot Jupiters, e.g., [Winn] |2010[ ). 

From the brightness distribution and the orbits, we model the planet’s flux temporal 
evolution—which accounts for the light-travel time across the system—^by sampling its 
surface with a grid of 2A^ points in longitude (0) and N + 1 points in latitude (9). We fix 
N to 100 to mitigate numerical effects up to 10“® the secondary eclipse depth, i.e. at least 
two orders of magnitude below the current photometric precision. 


Non-Uniform Brightness Models As discussed in Section [2^ it is critical to enable the 
mapping method to investigate different types/groups of maps in order to cope consistently 
with the degeneracy intrinsic to eclipse mapping. Therefore, we use three groups of bright- 
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ness models: (1) ID models, (2) spherieal harmonies, and (3) toy models. (1) The ID 
models are preferably used when a planet ean solely be mapped from its phase-eurve—i.e., 
when the available mapping information is mainly longitudinal (see Seetion [23| ). The ID 
models ean be subdivided in two model fa mi lies: one using n longitudinal bands with fixed 


positions on the dayside (“beaeh-ball” similar to the ones introdueed in |Cowan et al.[|2009| ) 
and another using longitudinal bands whose positions and widths are jump parameters in 
the MCMC fit. (2) The spherieal harmonics require two additional jump parameters for the 
direction (longitude and latitude) and a linear parameter for the amplitude, per harmonic. 
The generalized formulation of this brightness distribution model is; 


VsHA(t>^0) = 


(2.4) 


/=o 


where {cf) — — /S.9i) is the real spherical harmonic of degree I and order 0 , di¬ 

rected to A 0 J in longitude and in latitude. Ii is the /*^-harmonic amplitude that can 
be estimated at each step of the MCMC using the perturbed SVD method. As an example, 
the additional parameters for a dipolar fitting model compared to the conventional analy¬ 
sis method are: A 0 i and A 6 *i as jump parameters and Ji as a linear coefficient— Iq is the 
amplitude of the uniformly bright mode, included in the conventional analysis method (as 
< Jp >/</*> in Equation [23] ). 

(3) The toy models are developed specifically to model/retrieve one—or several—cold 
or hot spot. Their expression can be written as; 


ri(0,0) 

T2{4>A) 


/i0"exp [-Ao] cos^^o -t- Jo, 

Jl cos " 0 O COS ^ 6 o + lo if 0 O > 0 
Ji cos cos '’'6*0 + Iq if 0 o < 0 


1 1 exp 
1 1 exp 


A) ^ 

a ) 


exp 


r / X 21 


/ X 2I 

1- 

1 

exp 

1 

1 

1_ 


10 if 0O A 0 
-f Jq if 00 < 0 


(2.5) 

( 2 . 6 ) 


(2.7) 


where 0 o and 6 *o are respectively the longitude and latitude relative to the position of the 
model extremum, i.e. 0o = /(0, A0, a, (3, 7 ) and 60 = f{0, A 6 , a, 0, 7 ). A0 and A 6 * are 
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respectively the longitudinal and the latitudinal shift of the model peak from the substel- 
lar point. a,l3 and 7 parametrize the shape of the hot/cold spot. These toy models add 
to the conventional analysis method five jump parameters (A0, A 6 *, a,/3, 7 ) and a linear 
coefficient (Ji) per spot modeled. 
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2.3.3 Systematics and Correlated Noises 


Observations do not contain solely the signals targeted and described above. As an exam¬ 
ple, the time-dependent sensitivity of the telescope detector can significantly contribute to 
the observations. Figure [2^ presents an example of raw light curves (panel a), binned down 
to reveal an eclipse and the asymptotic behavior of the measured electron flux that is due to 
the time-dependent sensitivity of the detector (panel b), and the binned and corrected light 
curve that contains solely the targeted signals—the stellar and planetary fluxes—(panel c). 


Instrumental Systematics 


The analysis method is currently capable of analyzing reduced data from Infrared Ar¬ 
ray Camera (IRAC: [Fazio et ^ |2004| ) on board of the Spitzer Space Telescope ( jWerner 


et am2004| ). Therefore, the only instrumental systematics currently implemented relate to 
Sp/tzer/IRAC—or are systematics with standard functional forms. Hence, applications to 
datasets from other facilities may require the use of data reduced and corrected for system¬ 
atics (e.g., see Section [23| ). 

IRAC instrumental systematic variations of the observed flux, such as pixel-phase or 


detector ramp, are well-documented (e.g.. Desert et al. 2009, and references therein). 


Intrapixel-Sensitivity Induced Flux Variations are due to change in position of the 
target’s point-spread function (PSF) over a detector with non-uniform intrapixel sensitivity. 
Hence, these flux variations are strongly correlated with the PSF position (see Figure [2^ . 
For IRAC, only the channels at 3.6 and 4.5 /im—both using InSb-hased detectors—have 
shown significant effects of intrapixel-sensitivity induced flux variations. Most analysis 
methods perform a parametric correction based on a polynomial of the stellar centroid 


position ( jCharbonneau et ah]|2005[[2008| [Knutson et al.[|2007[[Gnion et al.[|2010b| ). Here, 


we implement the same “pixel mapping” method as in Lewis et al. (2013b), which was 


introduced in Ballard et al. (2010). The pixel mapping method consists in deriving at each 


step of the MCMC the sensitivity function, W (a;*, yi )—where xfyi are the stellar centroid 
positions in the frame—from the residuals using a Gaussian low-pass filter. In practice. 
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4280.4 4280.45 4280.5 4280.55 4280.6 

BJD - 2450000. 


(a) 




Figure 2-8: Light curve: noise and instrumental systematics. Raw time series (dark sym¬ 
bols) of the AOR r22809088 (see Section [2A| ) and the best fit (red curves), [^presents the 
raw time series obtained using aperture photometry, |(b)| the raw time series binned 64 times 
and fitted, [^presents the “detrended” time series binned in 3.5min intervals and fitted. 


43 













the pixel map is derived from 


W{xi,yi) 


X Fo, 


( 2 . 8 ) 


where 


Ki(j) = exp 


(xj-Xif {yj-ViY 


2^2 


2 ^^ 


(2.9) 


is the Gaussian low-pass filter, Fqj is the flux measured in the frame, and and ay are 
the widths of the Gaussian filter in the x and y direetions. In theory, Ud is the number of 


frame, but we adopt here the same proeedure as in Lewis et al. (2013b I to prevent the pixel 
mapping to be eomputationally prohibitive—see their Appendix B. 

Flux Ramp is due to the time dependenee of a detector response. Most IRAC’s data 
exhibit a ramp-like behavior (see e.g., Figure jl^Sb] ). For IRAC 8 /rm-channel—using Si : 
As-based detector, this effect is due to charge-trapping ([Knutson et al.[ [2007 [ [Agol et ah 


20101. For IRAC 3.6 and 4.5 /rm-channels, the physical origin of the ramp is not well- 


documented yet, although is can be adequately correct (see e.g., [Lewis et al.[[20l3b| ). Our 
method includes several functional forms to correct for flux ramps: linear, quadratic, and 
logarithmic functions of time. 


Possible Additional Systematics 

Our method can also account for possible low-frequency noise sources (e.g., instrumental 
and/or stellar) with time-dependent polynomials. 


Correlated Noise 


We take into account the remaining correlated noise following a procedure similar to Winn 


et al. (2008) for obtaining reliable error bars on our parameters. For each light curve, we 


estimate the standard deviation of the best-fitting solution residuals for time bins ranging 
from 3.5 to 30 minutes in order to assess their deviation to the behavior of white noise with 
binning. For that purpose, the following factor (3red is determined for each time bin: 
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BJD-2455200 


Figure 2-9: Intrapixel-sensitivity induced flux variations (Data from Sp/tzer/IRAC’s AOR 
r39136768). Top panel: Flux measured. Central panel: PSF position in pixel along the 
x-axis of the detector. Bottom panel: PSF position in pixel along the y-axis of the detector. 
The measured flux is strongly correlated with the target PSF position. 
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( 2 . 10 ) 


„ aN N{M - 1) 


where N is the mean number of points in eaeh bin, M is the number of bins, and ai 
and (Jat are respeetively the standard deviation of the unbinned and binned residuals. The 
largest value obtained with the different time bins is used to multiply the error bars of the 
measurements. 


2.3.4 Radial Velocity Measurements 

The retrieval method also enables fitting for radial veloeity (RV) measurements as they pro¬ 
vide additional eonstraints that are key to mitigate the degeneraeies intrinsie to eelipse map¬ 
ping (see Seetion [2!2.3[ ). The method ean fit for RV data separately, using adequate priors 
on the system parameters not eonstrained by sueh dataset (namely, Rp/Rl, b, y/ecos{u), 
-N/esin(ci;), p* and M*, in our parameter set). 

The method fits for the radial velocity projection 

Vr = Vz + Vzt + K{cos{u + f) + ecos{u)), (2.11) 


where Vz is the motion of the system barycenter, Vz is the long term velocity drift—that 
accounts for the possible effect of an yet-undetected companion with a large orbital period, 
K is the velocity semi-amplitude, uj is the argument of the periapsis and / is the true 
anomaly. K can be expressed in terms of the system parameters as follow 


K = 


Mr, 


sinu 


(Mp + M*)i Vl-e^ 


,27rG,i 


( 2 . 12 ) 


where Mp and M* are the planet and host-star masses respectively, i is the orbital incli¬ 
nation, e is the orbital eccentricity, G is the universal gravitational constant, and P is the 


orbital period (Winn 2010). 


There are no additional “conventional” jump parameter required to fit RV data simulta- 
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neously to light curves. The reason is that the other additional parameters V^, and K 
are linear parameters and ean henee be derived using SVD (see Seetion [2.3.1[ ). Note that K 
leads direetly to Mp. 

The eurrent version of our method does not fit for the Rossiter-MeLaughlin effeet 
(Rossiter 1924[ MeLaughfin[ 1924; Queloz et al. 2000 Triaud et al.[ 20091, henee mea¬ 
surements obtained during transit have to be disearded. 


2.3.5 In Practice: Fitting Model & Data Quality 

Fitting Model In praetiee, the fitting model, /sim(x), in our MCMC implementation 
works as follow. At eaeh step of a ehain, (1) the system orbits are simulated, (2) the differ¬ 
ent modes (stellar/planetary eontributions and the systematies) that eompose the simulated 
fight eurve and/or the radial veloeity data are simulated and (3) the mode’s amplitude are 
determined based on their simulated fight eurve by least squares minimization using the 
perturbed SVD method. 


Data Quality Eelipse mapping requires neat light-eurves. The oeeultation shapes need 
to be monitored properly, both in terms of time sampling and measurement preeision. We 
advoeate that the minimum requirements to perform eelipse mapping are: 

1. Time bins with a duration < 5% of the ingress/egress duration (i.e., > 20 bins to 
sample the ingress/egress). 


2. A RMS per time bin < 3% of the target’s oeeultation depth. 


Examples of datasets fulfilling sueh requirements are shown in Eigures [2-10| and 4-2 We 
show in Seetion|4^ that sueh requirements are suffleient to enable the loealization of the 
peak in the day side brightness distribution with an uneertainty < 3° in longitude and <10° 
in latitude and a preeision better than 17° on the latitudinal extent of the hot-spots and 
better than 8° on the longitudinal one, i.e., a preeision relevant for eomparisons with three- 
dimensional hot-Jupiter GCMs (Showman and Polvani 2011[ Heng and Workman 20141. 
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2.4 Application I: HD 189733 b 


The subset of exoplanet infrared (IR) observations aiming at mapping exoplanets is grow¬ 


ing. Currently, the Spitzer Space Telescope (Werner et al. 2004) has observed thermal 


phase eurves for a dozen different exoplanets as well as the IR oeeultations of over fourty 
exoplanets. Among these, HD 189733b ( Bouehy et al.[[2005| ) is arguably one of the most 
favorable transiting exoplanet for detailed observational atmospherie studies; in partieular, 
beeause its K-dwarf host is the elosest star to Earth with a transiting hot Jupiter. This means 
the star is bright and the eelipses are relatively deep yielding favorable signal-to-noise ratio 
(SNR). As sueh, HD 189733b represents a “Rosetta Stone” for the field of exoplanetology 
with one of the highest SNR seeondary eelipses ( |Deming et al. 2006 Charbonneau et al]| 
2008[ ), phase-eurve observations ( |Knutson et al.[ |2007[ |2009[ |2012| ) and, eonsequently, 
numerous atmospherie observational eharaeterizations (e.g., Grillmair et al.[ [2007 ; [Pont 


et al.i|200T[|Tinetti et al.H2007t|Redfield et al.i[2008HSwain et al.H2008HMadhusudhan and 


Seager 

2009 

Desert et al.[ 2009[ Deroo et al.[ 2010[|Sing et al. 

2011 

; Gibson et al.[ 2011 


Hudson et al. 2012). Although HD 189733b’s atmospherie models are in qualitative agree¬ 
ment with observations, important diserepaneies remain between simulated and observed 
light eurves and emission speetra (see e.g., [Showman et aL| |2009[ Figures 8 and 10). In 
addition, diserepaneies exist between several published inferenees—in partieular moleeu- 
lar deteetions—whieh emphasize the importanee of data reduetion and analysis proeedures 
(e.g., [Desert et al.[[2009[ [Gibson et al.[|2011| ). Henee, we undertake a global analysis of all 
HD 189733’s publie photometry obtained with the Spitzer Space Telescope for assessing 
the validity of published inferenees. 


In this Seetion, we present the first seeondary eelipse seanning of an exoplanet showing 
a deviation from the oeeultation of a uniformly bright disk at the 6 a level, whieh we obtain 
from the arehived SpitzerllKKC S-pm data of the star HD 189733. We apply the method 
introdueed in Seetion [23[ to disentangle the possible eontributing faetors to the anomalous 
eelipse shape. As a result, we perform a new step toward mapping distant worlds by eon- 
straining eonsistently HD 189733b’s system parameters, shape, and brightness distribution 
at 8 pm. 
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At the time of submission to Astronomy & Astrophysics, we learned about a similar 
study by Majeau et al. ( 2012| ), hereafter Ml2, using the same data but different frameworks 
for data reduetion and analysis. Our study differs from Ml2 in three main ways. First, we 
find a deviation from the oeeultation of a uniformly bright disk at the 6 a level in eontrast to 


the < 3.5 a level deviation in the phase-folded light eurves from Agol et al. (20101 (see their 
Figure 12) used in Ml2. Seeondly, this deviation has multiple possible eontributing faetors 
and our study provides a framework for eonstraining eonsistently these eontributing faetors 
(see Seetion [23] ). Thirdly, and related to the seeond point, we do not eonstrain a priori the 
system parameters to the best-fit of a eonventional analysis, nor the orbital eoeentrieity to 
zero; instead, we estimate the system parameters simultaneously with the BD. We eompare 


the methodologies and the results in Seetion 2.4.5 


This Seetion begins with a summary deseription of the Spitzer S-pm data. In See¬ 


tion 2.4.1 we present our data reduetion and eonventional data analysis—i.e., for whieh we 
model HD 189733b as a uniformly bright disk. This first analysis reveals that HD 189733b’s 


oeeultation has an anomalous shape. In Seetion 2.4.2, we present and assess the robustness 


of the anomaly in HD 189733b’s oeeultation. In Seetion [2.4.3[ we analyze HD 189733b’s 
S-pm photometry—transits, oeeultations, and phase eurve—to investigate the origin of 
HD 189733b’s oeeultation shape. We purposely foeus on HD 189733b’s photometry to 
point out the eomplementary insights gained from a global analysis ineluding the RV mea¬ 
surements in Seetion [2.4.4[ We diseuss in Seetion l2.4.5l the robustness of our results. 

The prospeets of our mapping method are introdueed in Chapter]^ Seetion[43| 


2.4.1 Data Reduction & Conventional Analysis 

We present here the eonventional data reduetion and data analysis performed for deter¬ 
mining HD 189733’s system parameters (Equarion|2.3|) based solely on the Infrared Array 


Camera (IRAC: Fazio et al. 20041 S-pm eelipse photometry. We emphasize that we si¬ 
multaneously analyze the whole data set, instead of eombining eaeh separately-analyzed 
eelipse events. Sueh approaeh helps to mitigate the effeets of noise by extraeting simul¬ 
taneously the information eommon to multiple light eurves. Therefore, it also enables the 
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detection of previously-unnoticed signals in the data. We begin with a summary description 
of the data sets; then, we introduce the analysis method and the physical models used for 
the parameter determination. 


Data Description and Reduction 

The eight secondary eclipses and the six transits of HD 189733t:[^used in this study are 
described by Astronomical Observation Requests (hereafter AOR) in Table |2.1| The data 
were obtained from November 2005 to June 2008 with IRAC at 8 /xm. These are calibrated 
by the Spitzer pipeline version SI8.18.0. The SI8.18.0 version enables improvements in 
the quality of data reduction over the original published data sets that used older Spitzer 
pipeline version^ Each AOR is composed of frames; each of which corresponds to 64 
individual subarray images of 32x32 pixels. 

The data reduction consists in converting each AOR into a light curve, hence each 
images into a flux measurement. For that purpose, we follow the procedure introduced 

^°Data available in the form of Basic Calibrated Data (BCD) on the Infrared Science Archive : http: 
//sha.ipac.caltech.edu/applications/Spitzer/SHA// 

"http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/ 
iracinstrumenthandbook/73 


Table 2.1: HD 189733 b’s AOR description 


AORKEY(0^) 

16343552(0) 

20673792(P) 

22808832(0) 

22809088(0) 

22809344(0) 

22810112(0) 

24537600(0) 

27603456(0) 

22807296(T) 

22807552(T) 

22807808(T) 

24537856(T) 

27603712(T) 

27773440(T) 


PI 

Publication Ref. 

Data sets'’ (64x) 

D. Charbonneau 

Charbonneau et al. (20081 

1359 

D. Charbonneau 

Knutson et al. (2007 1 

1319 

E. Agol 

Agol et al. (20101 

690 


E. Agol 


Agol et al. (20101 


690 


Exposure time [s] 

Ol 

0.4 


0.4 


0.4 


Aperture [px] 


3.6 

4.8 


4.8 


4.8 


^ AORKEY target: T, O or P respectively transit, occultation or phase curve. 

Present AOR are composed of data sets, each data set corresponds to 64 individual subarray images 
of 32x32 px. 
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in 


Gillon et al. (2010c I, which is performed individually for eaeh AOR. We first eonvert 


fluxes from Spitzer units of speeifie intensity (MJy/sr) to photon eounts; then, we perform 
aperture photometry on eaeh image with the IRAF{^ /DAOPHOT software (Stetson 1987). 
We estimate the PSF eenter by fitting a Gaussian profile to eaeh image. We estimate the 
best aperture radius (see Table |2.l| ) based on the instrument point-spread funetion (PSFpj), 
HD189733b’s, HD189733’s, HD189733B’s and the sky baekground flux eontributions. 
For eaeh image, we eorreet the sky baekground by subtraeting from the measured flux its 
mean eontribution in an annulus extending from 10 to 16 pixels from the PSF eenter. Then, 
we diseard: 


• the first 30 minutes of eaeh AOR for allowing the deteetor sensitivity stabilization, 

• the few signifieant outliers to the bulk of the x-y distribution of the PSF eenters using 
a 10 (T median elipping, 

• and, for eaeh subset of 64 subarray images, the few measurements with diserepant 
flux values, baekground and PSF eenter positions using a 10 a median elipping. 

Finally, the light eurve eorresponding to an AOR is the time series of the flux values av¬ 
eraged aeross eaeh subset of 64 subarray images; while the photometrie errors are the 
standard deviations on the averaged flux from eaeh subset. 


Photometry Data Analysis 


See “eonventional” model in Seetion 2.3.2 We use the following Gaussian prior for HD 189733’s 


mass: M* = 0.84 ± 0.06Mq (Southworth 20101. We first estimate the system parame¬ 
ters setting e = 0 —similarly to Agol et al. ( 2010[ ), based on the small inferred value of 
{e cos u, e since} and theoretieal predietions advoeating for HD 189733b’s orbital eireular- 


ization (e.g., |Fabryeky[|20T^ . We ehose to perform a first analysis under the eireularized- 
orbit assumption to point out qualitatively the influenee of unassessed assumptions in the 
eontext of a eorrelated parameter spaee (see Seetions [2^ and [2A3] ). 


^^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa¬ 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

^^http://ssc.spitzer.caltech.edu/irac/psf.html 
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Figure 2-10: SpitzerllKAC 8-/rm HD 189733b’s transit and occultation ingress/egress pho¬ 
tometry binned per 1 minute and corrected for the systematics (black dots) with their 1 a 
error bars (red triangles) and the best-fitting eclipse model superimposed (in red). The 
green dots present the residuals from |Agol et al.| ( [2010[ ) (their Figure 12), obtained using 
an average of the best-fit models from 7 individual eclipse analyses. Left: Phase-folded 
transits show no significant deviation to the transit of a disk during ingress/egress. Right: 
Phase-folded occultation ingress/egress deviate from the eclipse of a uniformly bright disk, 
highlighting the secondary eclipse scanning of HD 189733b’s dayside. Top: Phase-folded 
and corrected eclipse photometry. Bottom: Phase-folded residuals. 


2.4.2 HD 189733b’s Anomalous Occultation Shape 


Significance 


We detect an anomalous shape for HD 189733b’s occultation (see Figure [27T^ bottom-right 
panel). In particular, we find that its occulation deviate from the one of a uniformly bright 
disk with a 6.2 a significance!^ 

The IRAC S-pm photometry of the eclipse ingress/egress corrected for the systematics 


and binned per 1 minute are shown in Figure 2-10, with the best-fitting eclipse model super¬ 
imposed (solid red line). The error bars (red triangles) are rescaled by [i^red (~l-2) to take 
into account the correlated noise effects on our detection. In addition, we take advantage of 
the MCMC framework to account for the uncertainty induced by the systematic correction 


14 \ 


determine the significance of this structure as ^J2ieingress - Y^ieegress Tfdi, where 
and cTi are the flux measurement residual and its standard deviation at time i. 
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on the phase-folded light eurves. For that purpose, we use the posterior distribution of the 
aeeepted baselines to estimate their median instead of using the best-fit model—whieh has 
no statistieal signifieanee. We also propagate the uneertainty of the systematie eorreetion 
using the standard deviation of eaeh bin from the baseline posterior distribution—errorbars 
inereased by up to 20%. 

We inelude in Figure [2^d0| the residuals from Agol et al. ( 2010[ ) (their Figure 12). These 
are obtained using an average of the best-fit models from 7 individual eelipse analyses— 
not from a simultaneous analysis of all eelipses. For further eomparison, the strueture 
deteeted in HD 189733b’s oeeultation ean be translates into a uniform time offset of 37±6 


see (~6cr), in agreement with |Agol et al.| p010| ) estimate of 38±11 see (~3.5(j), light 
travel time dedueed. 


Robustness 


We test the robustness of our results against various effeets ineluding the baseline mod¬ 
els, the limb darkening, HD 189733b’s eireularized-orbit assumption and the AORs. Al¬ 
though below we mainly diseuss the robustness of the strueture in oeeultation, the ro¬ 
bustness of system parameter estimates is assessed simultaneously, but taekled in detail 
in Seetion l2!4.31 


As mentioned in Seetion 2.3.3 , Agol et al. (20101 advoeate using a double-exponential 
for modeling the deteetor ramp; but we use the quadratie funetion of log{dt) introdueed 


by Charbonneau et al. (20081. The reason is that by taking advantage of our Bayesian 


framework we show that the quadratie funetion of log{dt) is the most adequate for eor- 
reeting the present AORs. In partieular, we use two different information eriteric[^ that 
prevent from overfitting based on the likelihood funetion and on a penalty term related to 
the number of parameters in the fitting model. We obtain both a higher BIG (A BIG ~90) 
and a higher AIG (A AIG ~1.3) with the double-exponential; this means the additional 
parameters do not improve the fit enough, aeeording to both eriteria. Henee, the most ade¬ 
quate ramp model to eorreet these AORs is the less eomplex quadratie funetion of \og{dt). 


use as information criteria the BIC and the AIC (see, e.g., Gelman et al. 20041. The penalty term 
for additional model parameters is larger in the BIC than in the AIC. 
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Nevertheless, we assess the robustness of our results to different baseline models ineluding 
the double-exponential ramp, phase-pixel eorreetions, sinusoidal terms, and polynomials. 
These different MCMC simulations do not signifieantly affeet the anomalous shape found 
in oeeultation ingress/egress to within 0.5 a. The reason is that the typieal time seale of the 
baseline models is mueh larger than the timeseale of the struetures deteeted in oeeultation 
ingress/egress. Therefore, baseline models have a mitigated influenee on the ingress/egress 
shape (and henee on the retrieved maps), and vice versa. 

We assess the influenee of the priors on the limb-darkening eoeffieients; for that pur¬ 


pose, we perform MCMC simulations with no priors on ui and U 2 (see Seetion2.3.3). 
Again, we observe no signifieant influenee to within 0.5 a. We emphasize the neeessity of 
preeise and independent eonstraints on the limb-darkening eoeffleient, even if these might 
appear to be well-eonstrained by high SNR transit photometry. 

We assess the influenee of assuming HD 189733b’s orbital eireularization. For that 
purpose, we perform MCMC simulations with a free eeeentrieity. We observe no signifieant 
influenee to within 0.5 a for the jump parameters, exeept for the stellar density, the impaet 
parameter, ^/ecosu, and a/c since (see Seetion [2^ . In partieular, we observe a net drop 
of the anomalous shape signifieanee beeause of the potential for ^/e cos u and -y/e sin u to 
eompensate partially for deviations in ingress/egress. 

Finally, we also validate the independenee of our results to the set of AORs ineluded 
by analyzing different subsets of seven out of the eight eelipses. We observe relative sig¬ 
nifieanee deerease of ~ y/T/S; whieh are eonsistent with a signifieanee drop due to a re- 
duetion of the sample size. In partieular, it shows that the anomalous shape in oeeultation 
ingress/egress is not due to one speeifie AOR. 


Possible Contributing Factors: HD 189733b’s shape 

HD 189733b’s oeeultation shape is in agreement with the expeeted signature of the offset 
hot spot indieated by HD I89733b’s phase eurve ( {Knutson et al.[ 20071. However, it eould 
also be due to HD 189733b’s shape or biased estimates of b and p* resulting from the 
eireularized-orbit assumption (see Seetion [X2| ). 

We ean assume that HD 189733b’s projeetion at eonjunetions is a disk for two reasons. 
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(1) The transit residuals show no anomalous strueture (Figure [2-10[ bottom-left panel). This 
is in agreement with eurrent eonstraints on the HD 189733b oblateness (projeeted oblate¬ 
ness below 0.056, 95% eonfidenee, Carter and Winn[ 2010) and wind-driven shape (ex- 


peeted to introduce a light curve deviation below 10 ppm, see |Barnes et al.[[2009] ). In partic¬ 
ular, our transit residuals constrain HD 189733b’s oblateness below 0.0267 in case of a pro¬ 
jected obliquity of 45° and below 0.147 in case of a projected obliquity of 0°—95% confi¬ 


dence, Figure 1 of Carter and Winn (2010). No significant deviation in transit ingress/egress 
means that HD 189733b’s shape-induced effects in occultation ingress/egress are negligi¬ 
ble. Indeed, shape-induced effect in occultation ingress/egress are expected to be about one 
order of magnitude lower than in transit (effect ratio a< Jp > / < /* >, see Section [Z2.3[ ). 
(2) HD 189733b’s phase curve shows mainly a P-modulation—while shape-induced mod¬ 
ulation are expected to have a period of P/2; what is in agreement with our previous con¬ 
straint on HD 189733b’s shape. We therefore assume HD 189733b to be spherical, for a 
further analysis of the present data. 


2.4.3 Non-Conventionnal Analysis of HD 189733b’s Photometry 

Our second analysis aims to investigate the origin of HD 189733b’s occultation shape us¬ 
ing its IRAC 8-/im photometry—transits, occultations, and phase curve. We discuss in 


Section 2.2.3 the need for global analyses to consistently probe the parameter space while 
relaxing unverified assumptions possibly biasing the parameter estimates. However, we 
purposely focus the following analysis on HD 189733b’s photometry to emphasize latter 
the complementary insights gained while performing a global analysis including the RV 
measurements. 


HD 189733b’s Photometric Data 

In this second analysis, we use HD 189733b’s corrected and phase-folded eclipses from 
Section [Z4.2| (Figure |2-10| ) and its phase curve corrected for stellar variability ( |Agol et aT 


|2010| Figure 11). As discussed in Section 2.4.2 the transit and the phase curve enable to 
constrain respectively HD 189733b’s shape at conjunctions and its brightness distribution. 
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Note that we diseard the first third of the phase curve since it is strongly affected by the 
detector stability. 


To constrain the system parameters, we choose to use HD 189733b’s phase-folded tran¬ 
sits instead of using the system parameter estimates from Section 2.4.1 (see Table |2.2[ col¬ 
umn 2) for two reasons. First, these estimates are under the form of ID Gaussian dis¬ 
tribution while light curves provide a complex posterior probability distribution (PPD) 
over the whole parameter space. Secondly, these estimates are affected by HD 189733b’s 
crrcularized-orbit assumption (see implications in Section [2.4.3| ). For those reasons, it is 
relevant to simultaneously analyse HD 189733b’s scans and estimate the system parame¬ 
ters, to be consistent with our methodology of a “global” analysis and to avoid propagation 
of bias through inadequate priors, i.e., inadequate assumptions (Section[2.2.3[). 


Analysis Method 


See Section 2.3.1 and note that because we use the phase-folded light curves, the main 
constraint on the orbital period is missing. We, therefore, use a uniform prior on P based 
on the conventional-analysis PPD for P. In particular, we use a large prior (with an arbitrary 
symmetric extension of 10 up, where ap is the estimated uncertainty on P, see Table |2.2| ) to 
prevent our results to be affected by the assumption underlying this estimate. Note however 
that this has no influence on our further results—the reason is that the orbital period is 
highly-constrained by the transit epochs and, therefore, is not affected by secondary effects 
on the occultation shape. 


Results 

For HD 189733’s system, we find that relaxing the eccentricity constraint and using more 
complex brightness distributions (BDs) lead to lower stellar—and, hence, planetary—density 
and a more localized and latitudinally-shifted hot spot. We find that the more complex 
HD 189733b’s brightness model, the larger the eccentricity, the lower the densities, the 
larger the impact parameter and the more localized and latitudinally-shifted the hot spot es¬ 
timated. As discussed in Section [2^ the correlation between the planet brightness(/shape) 
and the system parameters e-oj-b-p^ is of primary importance for data of sufficient quality. 
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Table 2.2: Fit properties for different fitting models of HD 189733’s photometry in the 
SpitzerllKAC 8-/im channel 


Parameters (units) 

Uniform brightness 

Unipolar brightness 

Multipolar brightness 

e = 0 

e free 

^SH,l 

r2 

^SH,2 

^SH,3 

b{R,) 

n f?c: 7 fi_L 0.0021 
U.D0(0=tQ 0021 

o.6579±8:88ii 

o.6598±8:88it 

o.67i9±8:88?i 

o.6609±8:88i? 

o.6683±8:88?i 

y/ecosu! 

- 

o.oo43±8:™i^ 

o.ooo7±8:88?g 

-o.ooo 2 ± 8 : 888 i 

o.oooi± 8 : 88 }l 

coo 

oo 

oo 

oo 

-H 

o 

O 

o 

o 

^/esmu! 

- 

-o.oo 8 ± 8 : 8 i 

o.oi 6 ± 8 : 8 i! 

o.i42±8:8I8 

o.o46±8:8i 

n 1 oi _|_0.036 

U.lZizto.o64 

P*{Pq) 

1 Q1fi_L0.015 

i.yioiCo 016 

1 QQO + O.OIS 
i.yozzco 015 

1 qiq_l0.018 

l.yiOiCo 031 

i. 8 i 6 ± 8 : 8 !i 

1 Qf)q-L0.023 
i.yuyiCo 051 

1 OdC-LO.Oei 

i.O^Oito 058 

ABIC/AAIC 

0/0 

0.9/-0.4 

-170.1/- 173.3 

-167.2/- 172.3 

-167.4/- 172.0 

-170.0/- 175.8 


Rl/Rl = 0.024068±O:“ ; P{days) = 2.2185744±0:— ; 

Toihjd - 2453980) = 8.803352±[j:™o[]5f and = 0.0034117±[j:™oo37 


We present in this section our results for increasing model complexity to gain insight into 
the influence of the model underlying assumptions—e.g., the circularized-orbit assumption 
and the uniformly-bright exoplanet assumption. 

We gather the system parameter estimates for different fitting models in Table |2.2[ it 
shows the median values and the 68% probability interval for our jump parameters. We 
compute our estimates based on the posterior probability distribution (PPD) of global 
MCMC simulations, i.e., not as a weighted mean of individual transit or eclipse analy¬ 
ses. Our conventional analysis estimates are in good agreement with previous studies (e.g., 

. We discuss further the system 

parameter estimates obtained from our second analysis. 


Winn et aT| |2007[ [Triaud et~aL| |2009[ |Agol et aL| [2010 


Assuming HD 189733b to Be Uniformly Bright We show here the similar effects of 
the circularized-orbit assumption and the uniformly-bright exoplanet assumption. Both 
assumptions prevent from exploring dimensions of the parameter space and, therefore, lead 
to more localized, and possibly biased, PPD. 

We first present an unusual correlation between ^/e cos u and ^/e sin u (see Figure]^ 
1 la| ). It emerges from the partial compensation of the anomalous occultation enabled pri¬ 
marily by adequate combinations of y/ecosu and a/c since. As introduced in Section [2^ 
these parameters enable both to shift the occultation and change its duration (see Equa- 
tions |2.1| and |2.2[ respectively). 

Similarly, p* and b are also affected by the occultation shape. We emphasize this point 
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Figure 2-11: Effect of assuming an exoplanet to be uniformly bright and its orbit to be cir¬ 
cularized. |(a)| Marginal posterior probability distribution (PPD, 68%- and 95%-confidence 
intervals) of ^/ecosoj and i/esino; that shows an unusual correlation. |(b)| Marginal PPD 
(95%-confidence intervals) of p* and b that highlights the increase of adequate solutions 
enabled by the additional dimensions of the parameter space probed when relaxing the 
circularized orbit assumption. 



Figure 2-12: Deviations in occultation ingress/egress from the median-fit model for the 
individual perturbations of b, ^/ecosu, y/esmuj and p* by their estimated uncertainty (see 
Table |2.2[ column 2). It outlines that the system parameters enable compensation of an 
anomalous occultation that emerges from, e.g., a non-uniformly-bright exoplanet (see Fig- 
2-2[) possibly leading to biased estimates of the system parameters (Section[2^. 


ure 
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in Figure [2^T^ (to be compare with Figure [2^ . Figure [2^T^ shows the deviations in occul- 
tation ingress/egress induced by individual perturbations of b, y/ecosu, ^/esmu and p* 


by their estimated uncertainty (see Table. 2.2, column 3). Because these parameters affect 
an occulation shape, there exist adequate combinations of these parameters that can mimic 
an anomalous occultation within the measurement precision (see Figure[22T0l bottom-right 


panel, and Figure 2-12). Therefore, the circularized-orbit assumption and the uniformly- 
bright exoplanet assumption also affect the marginal PPD of {p*, b} by inhibiting the explo¬ 
ration of dimensions of the parameter space. We present in Figure [2^1 lb the extension of 
the {p*, b} marginal PPD that results from relaxing the circularized-orbit assumption. Sim¬ 
ilarly, we expect that the relaxation of the uniformly-bright exoplanet assumption would 
significantly affect the system-parameter PPD—compare Figures 2^ and 2-12 That is 


why we advocate in Section 2.2.3 for performing global analyses that constrain simultane¬ 
ously the possible contributing factors to the shape of an occultation and, thus, prevent the 
influence of unverified assumptions. 


Assuming HD 189733b to Be Non-Uniformly Bright We emphasize here the effect of 
relaxing the uniformly-bright exoplanet assumption on the fit improvement and the system- 
parameter PPD. We present in Figure |2^T3 the phase-folded IRAC 8-pm photometry of 
HD 189733b, corrected for the systematics with the best-fitting eclipse models for a uni¬ 
formly (blue) and a non-uniformly (green) bright exoplanet superimposed. The best-fitting 
non-uniformly bright eclipse model is shown for the F sh,i model—chosen arbitrarily, as 
the non-uniformly-bright models provide similar fits (see Table |2.2[ ). In particular, these 
models are significantly more adequate according to both the BIC and the AIC, see Ta- 
ble |2.2| (odds ratio: ~ 10^®). 

First, we introduce the results obtained using unipolar (i.e., with one spot on the plan¬ 
etary dayside) BDs to gain insight into the influence of relaxing the uniformly-bright exo¬ 
planet assumption. Then, we introduce the results obtained using multipolar BDs to assess 
the validity of trends observed in the unipolar-model results. 
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Figure 2-13: Phase-folded IRAC 8-/rm HD 189733b’s photometry binned per 1 minute and 
corrected for the systematics (black dots) with their 1 a error bars (red triangles) and the 
best-fitting eclipse models superimposed. Left: Phase-folded transits. Right: Phase curve 
and phase-folded occultations that show the benefit of using non-uniform brightness model. 


Unipolar Brightness Distribution: We first present the effect of relaxing the uniformly- 
bright exoplanet assumption on the system parameters. For that purpose, we show in Fig¬ 
ures : 


2-14a and 2-14c respectively the marginal PPDs of { y/e cos u , -y/e sin u } and {p*, ^/e sin a;} 


for the i model, and the ones for the r 2 model in Figures 2-14b and 2-14d We 


superimpose in Figures 2-14c and 2-14d the 95% uncertainty interval obtained for the 
uniformly-bright model to extend our previous observations regarding the effect of relaxing 
unassessed assumptions to probe a correlated parameter space, here the circularized-orbit 


assumption (see Section2.4.3, Figure2-1 lb). We observe the increases of y/esinu and b 
and the decrease of p* while y/e cos u is constrained closer to zero (see Table |2.2| ). 

The reason is that the compensation of HD 189773b’s anomalous occultation is now 
also enabled by the non-uniform brightness models; which provide a better compensation 
than e solely—with y/e cos uj for conventional analysis. Therefore, numerous combinations 
of e- and brightness-based compensations are adequate, in agreement with the correlation 
introduced in Section [2^ as highlighted by the PPD in Figure 2-15a[ Finally, we note a 
progressive evolution of the system-parameter PPD with the brightness model complexity 
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Figure 2-14: Influence of the brightness model complexity on the system-parameter pos¬ 
terior probability distribution (PPD), using unipolar models, (a) & |(b)| Marginal PPDs 
(68%- and 95%-confidence intervals) of \fe cos uj and a/c sin uj for two unipolar brightness 
models, respectively, with a fixed and large structure (F sh,i) and with free-confinement 
structure (r 2 ). A comparison with Figure 2-1 la shows that ^ecosa; is constrained closer 


to zero. The reason is that non-uniform brightness models enable the exploration of addi¬ 
tional dimensions of the parameter space and, therefore, provide additional adequate com¬ 
binations of the other contributing factors to compensate an anomalous occultation. In 
particular, the uniform time offset is now mainly compensated by a non-uniformly-bright 
model, rather than by y/ecosoj as in conventional analysis. It shows also the evolution of 
the marginal PPD toward larger y/e sin u) when using more complex brightness models— 
which enable more localized brightness structure, (c) & |(d)| Marginal PPDs of p* and 
/esincc for, respectively, Tsh,i and r 2 . These show the effect of relaxing unassessed as¬ 


sumptions to probe consistently a correlated parameter space; in particular, it outlines the 
possible overestimation of p* by 5% (i.e., at 6 a of the conventional estimate) when relaxing 
the conventional assumption of uniformly-bright disk to model occultations. 
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Figure 2-15: Dependence and significance of the brightness peak localization. |(a)|Marginal 
PPD (68%- and 95%-confidence intervals) of the brightness peak longitude, A0, and 
ecoscu for the Tsh,i brightness model—the simplest non-uniform brightness model used 
in this study. shows the correlation between the brightness model and the uniform time 
offset. This correlation emerges from enabling compensation of HD 189733b’s anomalous 
occultation with a larger set of contributing factors, here by relaxing the uniformly-bright 
planet assumption. In particular, the occultation shape is now mainly compensated by the 
brightness distribution rather than by y/e cos cu as in the conventional analysis (see Fig¬ 
ure 2^ Equation |2.1| and Figure [2H^ . |(b)| Marginal PPDs (68%-confidence intervals) of 
the brightness peak localization for the Tsh,i and r 2 brightness models. It shows that the 
brightness peak localization is model-dependent—mainly as a result of the fact that the 
mapping problem has degenerate solutions that are favored differently by each model. For 
example, the longitudinal F sh,i peak localization is constrained by the phase curve because 
of its large and constant extension; while the free extension of the F 2 model relaxes this 


longitudinal constraint (see Section 2.4.3). 


(from uniform to r 2 ). We assess further the validity of these observations, using spherical 
harmonics of higher degree. 

We now turn to HD 189733b’s brightness distribution. We show the dayside estimates 
for the Tsh,i and r 2 models with their corresponding uncertainties in Figures 2-16| and 
2- 17| respectively; we focus on HD 189733b’s dayside as it is effectively constrained by 


the combination of the phase curve and the secondary eclipse scanning. In particular, 
note that Figures |2- 1 6| and [2- 17 1 present HD 189733b’s brightness relative to HD 189733’s 
hemisphere-averaged brightness in the IRAC 8-/im channel (i.e., /p(0, 9)/ < J* > Ispm)- 
In addition, the figures are time-averaged; our estimates aim to approach the global pattern 
of HD 189733b’s BD based on eight snapshots taken from November 2005 to June 2008. 
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Figure 2-16: Estimate of HD 189733b’s brightness distribution in the IRAC 8-/rm channel 
using the Tsh,i brightness model. Left: Relative brightness distribution of HD 189733b’s 
dayside. Right: Dayside standard deviation. Because of its fixed and large structure, the 
Tsh,i brightness model is well-constrained in amplitude (by the occultation depth) and in 
longitudinal localization (by the phase curve). However, it is less constrained in latitude 
(by the secondary eclipse scanning) than more confined model such as r 2 (see Figure lTHT] ). 
These model-induced constraints are observable on the dayside standard deviation; which 
is significantly lower than for more complex brightness models (see Figures[2- 17y2-20a|and 


2-20b). In addition, the standard deviation distribution for the F sh,i model is related to its 


gradient (with a larger variation from the brightness peak localization along the latitude 
axis than along the longitude axis), because the brightness distributions accepted along the 
MCMC simulations differ from each other mainly in latitudinal orientation (see Figure[2^ 


18a). 



Figure 2-17: Estimate of HD 189733b’s brightness distribution in the IRAC 8-/im chan¬ 
nel using the r 2 brightness model. Left: Relative brightness distribution of HD 189733b’s 
dayside. Right: Dayside standard deviation. Because of its increased complexity, the F 2 
brightness model enables more localized structure that are less constrained in amplitude 
(by the occultation depth) and in longitudinal localization (by the phase curve) than large- 
and-fixed-structure models such as, e.g., the F sh,i model (see Figure [27T^ . However, it is 
well-constrained in latitude by the secondary eclipse scanning that is sensitive to confined 
brightness structure. These model-induced constraints are observable on the dayside stan¬ 
dard deviation; which shows a maximum at the brightness peak localization and extended 
wings towards west and east along the planetary equator. The reason is that the bright¬ 
ness distributions accepted along the MCMC simulations mainly affect the former by their 
amplitude changes and the latter by their structure changes (see Figure|2-18b[). 
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Finally, these estimates eorrespond to the median and standard deviation of the map trials 
aeeepted along the MCMC simulations, similarly to our approaeh for the eorreeted and 
phase-folded light eurves (see Seetion |Z4.2| ). 

Both models retrieve a spatial feature in HD 189733b’s BD that eorresponds to a hot 
spot. The Tsh,i model retrieves a hot spot shifted to the east of the substellar point, see 
Figure [27T^ The F 2 model retrieves a hot spot shifted to the east of the substellar point but 
also away from the equator, see Figure [27T7| However, we eannot diseuss the direetion of 
this latitudinal shift due to a North-South ambiguity (E. Agol, private eommunieation). 

The BD estimates shown in Figures |2-16| and [2-17| are signifieantly different both in pat¬ 
tern and in intensity. These differenees are due to the intrinsie degeneraey of the mapping 
problem (see Seetion[2^; whieh requires the use of different mapping models to enable a 
thorough diseussion. For example, brightness models with non-eonstant strueture (“eom- 
plex”, i.e., in opposition to a dipole) are less eonstrained by a phase eurve that is only 
dependent on the hemisphere-integrated brightness. To emphasize these model-indueed 
eonstraints, we present in Figure [27T^ animations showing the eompilations of dayside BDs 
aeeepted along the MCMC simulations for the Tsh,i and r 2 models. These eompilations 
show that ( 1 ) the amplitude and ( 2 ) the longitudinal loealization for the Tsh,i model are 
more eonstrained than for the F 2 model (by the oeeultation depth and by the phase eurve, 
respeetively) beeause of its fixed and large strueture. However, (3) the Tsh,i model is 
less eonstrained in latitude (by the seeondary eelipse seanning) than the more eomplex F 2 
model whieh enables more eonfined struetures that induee larger deviations in oeeultation 
ingress/egress. For that reason, the brightness peak loealization for the r 2 model is well- 


eonstrained in latitude (see Figure 2-15b I, while for the F sh,i model it is well-eonstrained 
in longitude. 

These observations reeall that the problem of exoplanet mapping has degenerate so¬ 
lutions, henee, the brightness peak loealization eannot be eonstrained uniquely without a 
priori assumption (e.g., assuming a dipolar BD). Therefore, we will further refer to our 
brightness-distribution estimates instead of the brightness peak loealization; whieh is not 
representative of eomplex BDs, in addition to being model-dependent. Nevertheless, we 
propose in Seetion[2.4.5| another unidimensional parameter to replaee the brightness peak 
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(a) 



(b) 


Figure 2-18: Insight into the model-dependence brightness distribution estimates. An¬ 
imations showing the compilations of HD189733b’s dayside brightness distributions 
accepted along MCMC chains for the Tsh,i and r 2 models, [(a)|and|(b)|respectively. These 
animations show that (1) the amplitude and (2) the longitudinal localization for the Tsh,i 
brightness model is more constrained than for the T 2 model (by the occultation depth and 
by the phase curve, respectively) because of its fixed and large structure. However, (3) 
the Tsh,i model is less constrained in latitude (by the secondary eclipse scanning) than 
the more complex r 2 model, which enables more confined structures that induce larger 
deviation in occultation ingress/egress (schematic description in Figure[2^. 


{Videos available in the electronic version.) 
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localization. 

Finally, note that these model-indueed eonstraints are also observable on the dayside 
standard deviation; whieh is signifieantly lower for T sh,i model than for more eomplex 
models. In partieular, the standard deviation distribution for the Tsh,i model is related to 
its gradient—with a larger variation from the brightness peak loealization along the latitude 
axis than along the longitude axis, because the BDs aeeepted along MCMC ehains differ 


from eaeh other mainly in (latitudinal) orientation, see Figure 2-18a This is in eontrast 
with the standard deviation distribution for the r 2 model that shows a maximum at the 
brightness peak loealization and extended wings towards west and east along the equator; 
beeause the BDs aeeepted along the MCMC simulations mainly affeet the former by their 


amplitude ehange and the latter by their strueture ehange (see Figure 2-18b I. 
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Multipolar Brightness Distribution: We observe an evolution of our inferenees when 
inereasing the eomplexity of our fitting model. To assess the validity of this observation, 
we present here the results obtained when using spherieal harmonies up to the degrees 2 
(quadrupole) and 3 (oetupole). 


We present in Figures 2-19a and 2-19e respeetively the marginal PPDs of {-y/e cos u , ^/e sin a;} 
and {p*,-y/esince} for the Tsh ,2 rnodel, and in Figures 2-19b and 2-19d for the Tsh ,3 
model. These PPDs appear as intermediate steps between the results obtained with the 
Tsh,i and r 2 models (see Figure [22T4] ). This gradual evolution reemphasizes the strong 
influenee of underlying model assumptions in the eontext of eorrelated parameter spaee. 

We show the day side brightness estimates for the Tsh ,2 and T sh ,3 models with their 
uneertainty in Figures [2-20a and 2-20b As observed for the system-parameter PPDs, the 
brigthness distributions also appear as intermediate steps. The major evolutions are the 
shrinking of the strueture retrieved and its shift away from the equator. 

These progressive evolutions of both the system-parameter PPD and the retrieved bright¬ 
ness strueture show that, for HD 189733’s system, relaxing the eeeentrieity eonstraint and 
using more eomplex BDs lead to a lower stellar density and more loealized and latitudinally- 
shifted hot spot. In partieular, we find that the more eomplex HD 189733b’s brightness 
model, the larger the eeeentrieity, the lower the densities, the larger the impaet parame¬ 
ter and the more loealized and latitudinally-shifted the hot spot estimated. We diseuss the 
signifieanee of these findings in the next seetion. 
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Figure 2-19: Influence of the brightness model complexity on the system parameter pos¬ 
terior probability distribution (PPD), using multipolar models, (a) & |(b)| Marginal PPDs 
(68%- and 95%-confidence intervals) of ^/e cosuj and -\/e since for the quadrupolar and 
octupolar brightness models, respectively. |^ & [^ Marginal PPDs of p* and a/c since 
for the quadripolar and octopolar brightness models, respectively. These confirm the trend 
observed in Section [2.4.3| (Figure |2-14| ) towards larger -y/esince and lower stellar density 
when increasing the complexity of the brightness distribution. 
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(b) 


Figure 2-20: Estimates of HD 189733b’s brightness distribution in the IRAC 8-/im channel 
using multipolar brightness models. Left: Relative brightness distribution of HD 189733b’s 
dayside. Right: Dayside standard deviation, (a) Estimate using the E sh ,2 brightness model. 
(b)| Estimate using the E sh,?, brightness model. These estimates confirm the trend toward 


a more localized and latitudinally-shifted hot spot when increasing the brightness-model 
complexity. Together with Figure [27T9| these estimates outline that the more complex 
HD 189733b’s brightness model, the larger the eccentricity, the lower the densities, the 
larger the impact parameter and the more localized and latitudinally-shifted the hot spot 
estimated. 
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“The Most Adequate Model” 


We present in Seetion 2.4.3 the results obtained using several fitting-models for assessing 


the model-dependenee of our inferenees in the eontext of underdetermined problem—i.e., 
degenerate solution (Seetion [2!2] ). As a result, we quantitatively outline the eorrelation be¬ 
tween the system parameters and the BD of an exoplanet introdueed in Seetion[2^ In 
partieular, we show a progressive evolution of the system-parameter PPD and the BD esti¬ 
mate when inereasing the eomplexity of the fitting model. We diseuss below the relevanee 
of our model-eomplexity inerease, as it may ultimately lead to data overfitting. 

We take advantage of our Bayesian framework using the BIC and the AIC. Both in¬ 
formation eriteria are in favor of models that relax the assumptions of a eireularized orbit 
and a uniformly bright HD 189733b. In partieular, the Tsh,i and Tsh ,3 models are favored; 
the BIC insignifieantly favors Tsh ,3 (odds ratio ~1.02) while the AIC signifieantly favors 
r SH ,3 (odds ratio ~3.5). Beeause not deeisiv^^ the information eriteria only suggest that 
the Tsh ,3 model provides the most adequate eonstraints on HD 189733’s system. 

In partieular, the information eriteria suggest that HD189733b’s hot spot is shifted 
both east of the substellar point and away from the equator and HD 189733b’s density has 
been overestimated by 3.6%. Furthermore, it suggests that HD 189733b’s orbit is possi¬ 
bly not fully eireularized (e = 0.015±g ;q®2)> although its eoeentrieity is eonsistent with 
zero. This reealls that the assumption of eireularized orbit has to be eontinuously assessed 
with eonstantly-inereasing data quality; even for old hot Jupiters that may show a hint of 
eoeentrieity (e.g., CoRoT-16b, see |011ivier et al.[[2012| ). Finally, for data-quality reason, 
the interpretation of HD 189733b’s BD has to foeus on global trends: the presenee of an 
asymmetrioal hot spot. 

HD 189733b’s dayside presents a shifted hot spot. The eastward shift is in agreement 


with the literature: (1) with previous derivations, from HD 189733b’s phase ourve (Knutson 


et al. 20071 and an eolipse timing eonstraint (Agol et al. 20101; and (2) atmospherio mod¬ 


els suggesting a super-rotating equatorial jet (e.g., [Showman et al.[ |2009[ ). In opposition, 
the suggested shift away from the equator is new. The small-seale origin of this latitudinal 
asymmetry remains uneonstrained beeause we use large-seale brightness models to be eon- 


16 


'Note that theoretical studies favor the AIC, e.g., 


Burnham and Anderson 


(2002 i; Yang (2005 i. 
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sistent with the data quality. For that reason, additional observations would be required to 
improve our understanding of HD 189733b’s atmosphere (see Chapter^; in partieular, its 
interaetions with HD 189733 eould induee unexpeeted thermal patterns, e.g., asymmetrie 


patterns in its BD (Leeavelier des Etangs et ah, 20121. For example, magnetie star-planet 


interaetions may lead to energy dissipation due to the stellar field penetration into the exo¬ 
planet envelope (e.g., |Laine et al.[[2008| ) and to extensive energy injeetions into the auroral 


zones of the exoplanet from magnetie reeonneetions (e.g.. Ip et al. 2004)—similarly to the 


Jupiter-lo flux tube (e.g.. Bigg 1964). However, sueh magnetie reeonneetions have so far 


been only observed at the stellar surfaee, in the form of ehromospherie hot spots rotating 
synehronously with the eompanions (e.g., Shkolnik et al.[|2005 ; Lanza[ 2009). 


2.4.4 Global Analysis: Complementary Insights from RV Data 

We emphasize in our study the neeessity of global analyses to fit high-SNR oeeultation (see 
Seetion [Z2l ). For that reason, we purposely foeus in the previous Seetion on the photome¬ 
try (transits, oeeultations, and phase eurve) in order to emphasize here the further insights 
gained while performing a global analysis ineluding the RV measurements. RV measure¬ 
ments may provide a eomplementary insight into the e-a;-6-p*-BD eorrelation. Therefore, 


we analyse first HD 189733’s out-of-transit RV data (Winn et al., 2006 Boisse et al. 20091 


separately to assess the eonstraint derived solely from the RV data on e and oj. We present 
in Figure|2j^ our overall Keplerian fit and the marginal PPD of the parameters \fe cos u 
and since. This shows that the RV data does not eonstrain HD 189733b’s eeeentrieity 
further than the Spitzerl\K\C 8-/im photometry does for low-eomplexity brightness mod¬ 
els (see Figures 2-1 la[ 2-14a and 2-19a). However HD 189733’s RV measurements may 
affeet our inferenees for more eomplex models that favor a loealized hot spot and larger 


eeeentrieity (see Figures 2-19b and 2-14b), by rejeeting the solutions involving 


e since 


> 


0.15. 


We observe no ehange of the system parameters PPDs for simple brightness models, 
as expeeted from the eomparison of Figure 2^2 lb| and Figures 2- 14a| and 2-19a However, 
we observe the effeet of rejeeting the solutions involving a/c since > 0.15—in the eon- 
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Figure 2-21: Influence of HD 189733’s RV measurements on inferences obtained from 
the Spitzerll^hC 8-/im photometry. [^Overall Keplerian fit (green) of the out-of-transit 
RV data (black dots) with their 1 a error bars (red triangles) from [Winn et al. ( |2006 1 and 
Boisse et al. ( 2009| ). |(b)| Marginal PPD (68%- and 95%-confidence intervals) of y/ecosu 
and ^/e sin u obtained solely from the RV data. The RV data does not constrain further 
HD 189733b’s eccentricity than the 5pzYzer/IRAC S-fim photometry for low-complexity 
brightness models (see Figures 2-1 la[ 2-14a[ and 2-19a). However, for more complex mod¬ 
els that favor a localized hot spot and larger y^sincu, HD 189733’s RV measurements may 
affect our inferences by rejecting the solutions involving a/b since > 0.15 in the context of 


the e-u-b-pi,-BD correlation (see Section 2.2 j 











































Table 2.3: Parameter estimates obtained from the photometry and the RV measurements. 


Parameters (units) 

Unipolar brightness 

Fsr/.i r2 

Multipolar brightness 

^SH,2 ^SH,3 

HR.) 

^/e cos uj 
^/e sin Lu 

pApo) 

n ac;qi -l0.002S 

U.OOyi=to.o023 

0 001 

U.UUIZZCq QQ22 
f) Qr)7_L0.050 

U.UU(=to.o30 

1 qoo_L0.0i7 
i.yzzzCQ Q22 

o.eeoiiO;™!? 

-0.0002±O;OOi| 

0.031±O:OIi 

i-9i2±g:8i? 

o.659o±8:88ii 

o.ooo3±8:88i 

o.oo6±8:8i 

i.922±8:8^? 

o.6592±8:88i 

o.ooo3±8:88t8 

o.oo9±8:8il 

i.922±8:8JI 


text of the e-a;-6-p*-BD eorrelation—for eomplex brightness models (rsH ,3 and r 2 ). In 
partieular, the evolution of the parameter estimates—toward lower densities, larger impaet 
parameter and a more loealized and latitudinally-shifted hot spot—providing larger eeeen- 
trieity are mitigated. We present the influence on the system-parameter PPD in Figure 2-22[ 


The {y/ecosoj^ A/esincej-PPD (Figures 2-22a and 2-22b) conceptually corresponds to the 


marginalized product of the PPDs estimated using solely the photometry (see Figures 2-19b 


and 2-14b) and using solely the RV data (Figure2-21a). These highlight the redistribution 
of the probability density that results from the rejection of solutions involving ^/esinoj 
> 0.15 by HD 189733’s RV data. As a consequence of this probability redistribution in 
{-n/ccoso;,- y/e since}, the set of acceptable combinations to compensate HD 189733b’s 
anomalous occultation is reduced too. In other words, the b, p* and BD estimates are af¬ 
fected. On the one hand, the system parameter estimates are consistent with those obtained 
for less complex brightness models (see Table|2.3| and compare Figures 2-22c| and 2-22d 


to Figures 2-14c and 2-19c). On the other hand, HD 189733b’s dayside brightness esti¬ 
mates present less confined patterns (compare Figures 2-20b and 2-17 with Figures [2-23a 
and|2-23b] respectively). 


Conclusion 

We estimate that HD 189733b’s brightness distribution in the IRAC 8-/im channel is best 
represented by Figure jlH^ based on our global analysis of the Spitzer/IRAC S-pm pho¬ 
tometry and the RV measurements of HD 189733. In addition, we refine the constraint on 
HD 189733b’s orbital eccentricity to e < 0.011 (95% confidence). 
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Figure 2-22: Influence of the RV measurements on the system-parameter PPD estimated us¬ 
ing complex brightness models that suggest a non-zero eccentricity from the Spitzer/IRAC 
8-/im photometry separately (see Figures [T-19b| and |2-14b[ ). [(a)] & Marginal PPDs 
(68%- and 95%-confidence intervals) of y/e cos ui and 


/e sin uj for the F sh ,3 and the F 2 
brightness models, respectively. Conceptually these correspond to the marginalized prod¬ 
uct of the PPDs estimated using solely the photometry (see Figures |2-19b and 2-14b I and 
using solely the RV data (Figure 2-21 aj ). In particular, these highlight that the solutions 
involving y/esinuj > 0.15 are strongly rejected by HD 189733’s RV data, leading to a re¬ 
distribution of the probability density. |(c)| & |(d^ Marginal PPDs of p* and a/c since for the 
Tsh ,3 and r 2 brightness models, respectively. Th ese show the impact of RV data in the 
context of the correlation introduced in Section |2^ RV data favors solutions with y/e sin u 
<0.15 and, hence, solutions consistent with the inferences obtained with less complex 
brightness models (see Fig|2-14[). 
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(b) 


Figure 2-23: Influence of HD 189733’s RV measurements on HD 189733b’s dayside bright¬ 
ness distribution, estimated using complex brightness models that suggest non-zero eccen¬ 
tricity from the SpitzerllKAC 8-/rm photometry (see Figures IT- 1 9b| and |2-14b[ ). Left: Rel¬ 
ative brightness distribution of HD 189733b’s dayside. Right: Dayside standard deviation. 
[^Estimate using the F sh, 3 , brightness model, (b) Estimate using the F 2 brightness model. 
The influence of HD 189733’s RV data is an extension of the brightness patterns (compare 
with Figures 2-20b and 2-17, respectively) associated with a decrease of the brightness 
peak—conservation of the hemisphere integrated flux. These show that localized bright¬ 
ness pattern that are favored by the photometry are rejected by the RV data, because these 
are associated with ^fesmuj > 0.15. 
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2.4.5 Discussion 


Adequacy of Conventional Analyses 

We show based on high-SNR data the limitation of eonventional analyses to interpret light 
eurves beeause of they assume the planet to be uniformly bright despite the e-a;-6-p*-BD 
eorrelation. The significanee of this limitation is related to an oeeultation SNR—which is 
particularly high for HD 189733b’s oeeultation in the Spitzerl\KhC 8-/im channel. A “suf¬ 
ficient” SNR on the oeeultation requires resolving the oeeultation ingress/egress; therefore, 
it has to be about one order of magnitude less than the oeeultation depth, for a time bin 
about one order of magnitude less than the oeeultation ingress/egress duration. In other 
words, conventional analyses are typically adequate when the ratio eclipse depth to photo¬ 
metric precision on a ~l-min bin is ^ 10. For occultations with significantly lower SNR, 
the underlying assumption of conventional analyses have a negligible effect on the system 
parameter PPD, and hence conventional analysis are adequate. 

We briefly outline that the conventional assumption of a uniformly-bright exoplanet 
could also affect the inferred planetary interior models; because p* is possibly affected and, 
therefore, pp is too (e.g., suggested 3.6%-overestimation for HD 189733b’s). 


Reducing Brightness Distributions to Unidimensional Parameters 


We show in Section 2.4.3 that the light curve of an exoplanet does not constrain uniquely 
its brightness peak localization. We discuss here the reasons why we strongly advocate 
discussing the BD estimates and, if necessary, using with care the dayside barycenter as 
an alternative possibility for representative 1-D parameter. The reason is that the dayside 
brightness barycenter weights the BD according to the geometrical configuration at supe¬ 
rior conjunction (i.e., it contains partial 2D information). 

We show in Figure|2^^ the marginal PPDs (68%-confidence intervals) of the bright¬ 
ness peak localization for the Tsh,i and r 2 brightness models. A comparison with the 
marginal PPDs of the brightness peak localization (Figure |2-15b ) shows the reduced model- 
dependence of the dayside barycenter. In particular, it shows a less-extended PPD for the 
r 2 -model barycenter; because this extension for the PPD of the brightness-peak-localization 
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Figure 2-24: Reducing brightness distributions to unidimensional parameters. Marginal 
PPDs (68%-confidence intervals) of the dayside barycenter brightness peak localization for 
the Tsh,i and r 2 brightness models. A comparison with the marginal PPDs of the bright¬ 
ness peak localization (Figure 2-15b I shows the reduced model-dependence of the day- 
side barycenter. In particular, it shows a less-extended PPD for the F 2 -model barycenter; 
because this extension for the brightness-peak-localization PPD emerges from the model 
wings, weighted by the dayside barycenter. In addition, it shows the shift and slight shrink¬ 
ing of the r sH,i PPD that reflect the barycenter weighting according to the geometrical 
configuration at superior conjunction; map cells closer to the substellar point have more 
weight. 


emerges from the model wings—weighted by the dayside barycenter. In addition, it shows 
the shift and slight shrinking of the F sh,i PPD that reflects the barycenter weighting ac¬ 
cording to the geometrical configuration at superior conjunction; map cells closer to the 
substellar point have more weight. This emphasizes the primary drawback of the dayside 
barycenter that is to attenuate the offset of BDs. This recalls that unidimensional parameters 
cannot stand adequately for complex BDs and, therefore, have to be used complementary 
to BD estimates. 


Effect of HD 189733’s activity 


HD 189733 presents high-activity levels that may affect the transit parameters—inch, R^/Rl, 


due to occulted/unocculted star spots ([Pont et ar||2007[|Sing et al.[|2011[). Therefore, treat¬ 


ing coherently spots of active stars is key. However, while important in the optical, the 
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Figure 2-25: Influence of HD 189733’s variability on individual transit-depth estimates 
in the Sp/tzer/IRAC 8-/xm channel. Transit-depth estimated for six transits in individual 
(markers) and global (solid line, the dashed lines refer to the \-a error bar) analyses. Our 
estimates (green) show no significant transit-depth variation with the relative eclipse phase 
(blue numbers) in opposition to[Agol et~aL (2010) (red). 


stellar activity may be negligible at 8 /im with the current data quality. We assess this state¬ 
ment performing individual analysis of the 6 transits used in this study (see Table |2.l| ). We 
present in Figure 2^25] our individual transit-depth estimates. These estimates show no sig¬ 


nificant temporal variation (in opposition to Agol et al. 2010 who attributed these to stellar 
activity). Similarly, we observe no significant variation of the transit parameters from one 
individual analysis to another. In addition, we observe no pattern specific to the occultation 
of a star spot (i.e., similar to, e.g., “Features A and B” in the Figure 1 of [Pont et~am2007|). 


Therefore, we consider that our time-averaged inferences (see Section 2.4.3) are not biased 
by HD 189733’s activity. 


Comparison with Majeau et al. (2012) 


We obtain qualitatively similar results as M12: an offset hot spot. Nevertheless, significant 
differences exist between both studies. (1) The starting point of our mapping analysis is 
our detection of HD 189733b’s eclipse scanning at the 6a level (in contrast to their ~3.5 a). 
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Our study provides consistent estimates of the contributions of the multiple possible factors 
to HD 189733b’s occultation shape. (2) In addition, and related to the second point, we do 
not constrain a priori the system parameters to the best-fit of a conventional analysis nor 
the orbital eccentricity to zero; instead, we estimate the system parameters simultaneously 
with the BD. (3) We also investigate the model-dependence of our inferences because the 
mapping problem has degenerate solutions, in contrast to Ml2 who focused on a dipolar 
brightness model (see related discussion in Section[2.4.3|). In particular, we show their 


impact on the system-parameter PPD (see Section 2.4.3). 

As a consequence. Ml2 estimate the brightness peak localization with narrow error bars 
(21.8 ± 1.5°east and 3.1 ± 9.4° away from the equator) for a circularized HD 189733b’s 
orbit; while we show that the brightness peak localization, as well as the system-parameter 
PPD, are model-dependent because of the e-a;-6-p*-BD correlation. Nevertheless, for a di¬ 
rect comparison to M12’s estimate of the brightness peak localization for a dipolar bright¬ 
ness model, we estimate it to 11.5 ± 4.3°east and 3.1 ± 11.4° away from the equator, see 


Section 2.4.5 Figure 2-15b 
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2.5 Application II: Kepler-7 b 


The mapping work introduced in this Section was performed for and published in Demory 


et al. (2013). The overall analysis presented in Demory et al. (20131 not only presents the 
first evidence of optically thick clouds in an exoplanetary atmosphere but also the first cloud 
map of an exoplanet (see Figure [2^27] ). Kepler-7 b’s cloud map provides unprecedented 
insights into its climate: our current understanding is that its equatorial jet brings hot gas 
from its dayside to its nightside where the temperature drops enough to form clouds— 
possibly silicate-dominated clouds—that are then brought to the dayside where they reflect 
the host-star light until their properties change as they encounter higher temperature around 
the substellar point. 


Introduction 


Kepler-7b is a hot Jupiter orbiting a subgiant G star in 4.89 days (Latham et al. 20101. Its 
relatively low mass Mp = 0.44 ± 0.04Mj„p and large radius Rp = 1.61 ± 0.02i?j„p result 


in a very low density pp = {i.lAg.cm^ (Demory et al. 20111. Remarkably, Kepler-7b has 
a significant geometric albedo Ag ~ 0.35 and exhibits a clear phase-curve modulation in 
the Kepler bandpass ( jPemory et al.[ |2011[ [Kipping and Bakos] |2011[ [Coughlin and Ldpez^ 
Morales [ 2012). Kepler-7b’s effective temperature places it in an exceptionally rich region 
of condensation phase space. Because of the extreme difference between its equilibrium 
temperature and the brightness temperature as derived from its occultation in the Kepler 
bandpass, the origin of Kepler-7bs albedo has been attributed to the presence of a cloud or 


haze layer in its atmosphere or to Rayleigh scattering (Demory et al. 2011). 


Demory et al. (2013) presents a detailed study of Kepler-7 b’s observation using Spitzer 


and Kepler. (1) Description, reduction, and analysis of Spitzer' 3.6 and 4.5 pm photometry 
and (2) Description, reduction, and analysis of Kepler's photometry spanning over three 
years of quasi-continuous monitoring—which includes the mapping work based on the 
phase-folded light curve. Finally, combining the insight gained from infrared and visible 
dataset, Demory et al.| ( 20T^ constrain the origin of Kepler-7 b’s visible flux as scattered 
light from optically thick clouds and discussed the possible properties of those clouds. Only 
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Figure 2-26: Kepler-7b’s phase eurve based on Kepler Q1-Q14 data. Data are binned per 5 
minutes. The Lambertian sphere (red), l-free-band (blue), and 3-fixed-band (green) best-fit 
models are superimposed. 


the work eoneerning Kepler-7 b’s mapping is part of this thesis and is introdueed hereafter. 


For further details regarding the overall and multidiseiplinary analysis yielding to the first 


map of eloud in an exoplanet atmosphere, please refer to Demory et al. (2013). 


Mapping using Kepler^s Dataset 


We map Kepler-7b using its eorreeted and phase-folded phase eurve obtained in the visi¬ 
ble using the Kepler Space Telescope over more than three years (Figure[?-26|). The data 


reduetion and eorreetion is detailed in Demory et al. (2013), together with a diseussion of 
the robustness of Kepler-7b’s phase eurve over the different quarters. 

Kepler-7b’s phase eurve deviates from a pure Lambert-law phase-dependent behavior 
(e.g., |Sobolev[ |1975| ) expeeted for isotropie seattering alone (Figure [2^26| red). The main 
feature of Kepler-7b’s phase eurve is a delay of 13 ± 3.5 hr of the phase-eurve’s peak from 
the oeeultation eenter. This delay implies that the hemisphere-integrated flux is maximum 
to the west of Kepler-7b’s substellar point. We further measure a phase-eurve amplitude 
of 50 ± 2 ppm and an oeeultation depth of 48 ± 3 ppm, eorresponding to a geometrie 
albedo Ag = 0.35 ± 0.02. This oeeultation depth translates to a brightness temperature of 
2645±|o K in the Kepler bandpass, whieh is 1000 K and 800 K larger than the infrared 
brightness temperatures upper limits measured at 3.6 and 4.5/im respeetively (see Seetion 
2 of [Demory et al.[ 2013| ). The phase-eurve amplitude and oeeultation depth are in agree¬ 
ment with previous analyses (Demory et al.[ 201 1[[Kipping and Bakos[|2011 ; Coughlin and 
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Lopez-Morales 20121. The key features of Kepler-7b’s phase-eurve translate direetly into 


eonstraints on maps ( Cowan and Agol[ 20081 assuming a tidally looked planet on a oiroular 
orbit. Indeed, a planetary phase-eurve Fp/F^, measures the planetary hemisphere-averaged 
relative brightness < Ip > / < /* > as follows: 


a = 


Ip ^ (tt) / Rp 

< F> \R^ 


(2.13) 


where a is the orbital phase. 

We first notioe that Kepler-7b’s planetary flux oontribution starts from phase 0.18±0.03, 
when the meridian eentered 25 ±12° east of the substellar point appears. Seeond, the phase- 
ourve’s maximum is loeated at phase 0.61 ± 0.03, implying that the brightest hemisphere is 
eentered on the meridian loeated 41 ± 12° west of the substellar point. Third, the planetary 
flux oontribution vanishes around the transit, implying that the “bright” area extends up to 
the western terminator, while its extension to the east of the substellar point is nominal. We 
finally note that the phase-ourve’s amplitude of 50 ± 2 ppm oonverts into an hemisphere- 
averaged relative brightness 74 ± 2 x 10^ (Equation |2.13| ). 

We longitudinally map Kepler-7b using the MCMC implementation presented in See- 
tion|2.3.1j We use the ID models as Kepler-7b’s brightness distribution is mainly eon- 


strained longitudinally. We ohoose the two simplest models from these families: a 3-fixed- 
band model and l-free-band model so as to extraet Kepler-7b’s longitudinal dependenee 
of the dayside brightness as well as the extent of the “bright” area. For both models, we 
eompute eaeh band’s amplitude from their simulated light eurve by using a perturbed sin¬ 
gular value deeomposition method. The eorresponding median brightness maps are shown 
on Figure 2-27 The l-free-band model (Figure |2-26[ blue) finds a uniformly bright lon¬ 
gitudinal area extending from 105 ± 12° west to 30 ± 12° east with a relative brightness 
78 ± 4 X 10^ (Figure 2-27a). The 3-fixed-band model (Figure [2-26[ green) finds bands of 
relative brightness deereasing from the west to the east with the following values: 100 to 


68 and 3 ± 6 x 10^ (Figure 2-27b). We finally note that the 1-freehand model finds a bright 
seetor extending to the night side, due to the sharp flux inerease observed around transit 
(Figure [2-26[ ). 
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Figure 2-27: Longitudinal brightness maps of Kepler-7b. Kepler-7b’s longitudinal bright¬ 
ness distributions /p(a)//* as retrieved in Kepler’s bandpass using the l-free-band model 


(a) and the 3-fixed-band model (b) It shows that Kepler-7b’s is bright in Kepler's band 
in an area that extends from its dayside western terminator to ~ 30° east of its substellar 
point. 
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Chapter 3 


Constraining Exoplanet Mass from 
Transmission Spectroscopy^] 


The mass of a planet is a fundamental parameter beeause it is eonneeted to its internal and 
atmospherie strueture. A planet’s mass affeets basie planetary proeesses sueh as interior 
eooling and plate teetonies ( jStamenkovic et al^ |2012[ jStamenkovic and Brener] |2014| ), 
magnetie field generation (|Stevenson||2003[|Stamenkovic et al.[[2011|), outgassing ([Elkins 


and Seagerj [2008 ), and atmospherie eseape (Lammer et ah, 20091. Measurement of a 


planetary mass ean in many eases reveal the planet bulk eomposition, allowing to determine 
whether the planet is a gas giant or is roeky and suitable for life as we know it. 

Planetary mass is traditionally eonstrained with the radial veloeity (RV) teehnique, us¬ 
ing single-purpose dedieated instruments. The RV teehnique measures the Doppler shift of 
the stellar speetrum to derive the planet-to-star (minimum) mass ratio as the star orbits the 
planet-star eommon eenter of mass. Although the RV teehnique has a pioneering history 
of sueeess laying the foundation of the field of exoplanet deteetion, it is mainly effeetive 
for massive planets around relatively bright and quiet stars. Most transiting planets have 
host stars that are too faint for preeise RV measurements. For suffieiently bright host stars, 
stellar perturbations may be larger than the eompanion-indueed signal, possibly preventing 


a determination of the planet mass with RV measurements even for hot Jupiters (e.g.. Col¬ 


lier Cameron et al. 20101. In the long term, the limitation due to the faintness of targets 


‘Work published in part in Science, Volume 342, Issue 6165, pp. 1473-1477 (2013), see 


de Wit and 


Seager (20131. 
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will be reduced with technological improvements. However, host star perturbations may 
be a fundamental limit that cannot be overcome, meaning that the masses of small planets 
orbiting relatively quiet stars would remain out of reach. 


Current alternative mass measurements to RV are based on modulations of planetary- 
system light curves ( Mislis et al.[|2012 l or transit-timing variations ( Fabryck^ 2010). The 
former works for massive planets on short period orbits and involves detection of both 
beaming and ellipsoidal modulations (van Kerkwijk et aL||2010; Faigler and Mazeh 2011 


Shporer et al.[ 2011[). The latter relies on gravitational perturbations of a companion on 


the transiting planet’s orbit ( |Miralda-Escude[|2002^|Agol et al.[[2005| [Holman and Murray [ 


2005). This method is most successful for companions that are themselves transiting and 


in orbital resonance with the planet of interest (Carter et al., 20121. For unseen companions 
the mass of the transiting planet is not constrained, but an upper limit on the mass of the 


unseen companion can be obtained to within 15 to 50% (Steffen et al. 20131. 

Transiting exoplanets are of special interest because their densities can be estimated 
providing insights into their internal structures. Furthermore, the atmospheric properties 
of a transiting exoplanet can be retrieved from the host-star light passing through its at¬ 
mosphere when it transits. However, the quality of atmospheric retrieval is reduced if the 


planet’s mass is inadequately constrained ( [Barstow et al.[|2013a] ). 

In this Chapter, we show that a planet’s mass is a key parameter of its transmission 
spectrum that has to be accounted for by atmospheric retrieval methods to avoid the deteri¬ 
oration of the retrieval quality introduced by Barstow et al.| (2013a). Therefore, we develop 
MassSpec, a method for simultaneously and self-consistently constraining the mass and 
the atmospheric properties of an exoplanet based solely on transit observations, thereby 
enabling mass measurements for transiting planets for which the RV method fails. 


3.1 Mass Spec’s Concept and Feasibility 

The mass of a planet affects its transmission spectrum through the pressure profile of its 
atmosphere (i.e., p{z) where z is the altitude), and hence its atmospheric absorption profile. 
For an ideal-gas atmosphere in hydrostatic equilibrium, the pressure varies with the altitude 
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as dln(p) = —j^dz, where H is the loeal atmospherie seale height defined as 


H = 


kT 


(3.1) 


where k is Boltzmann’s eonstant and T, /i and g are the loeal (i.e., altitude dependent) 
temperature, mean moleeular mass, and gravity. By expressing the loeal gravity in terms 
of the planet’s sizes (mass, Mp, and radius, Rp), Equation|3.1|ean be rewritten as 


Mp 


kTRl 

fiGH' 


(3.2) 


Thus, our method eoneeptually requires eonstraining the radius of the target as well as its 
atmospherie temperature, mean moleeular mass, and seale height. Can these independent 
eonstraints be provided solely by a transmission speetrum? 


A planet transmission speetrum ean be seen as 
a wavelength-dependent drop in the apparent bright¬ 
ness of the host star when the planet transits (Fig¬ 
ure!; 


3-11. At a wavelength with high atmospherie ab- 


High 

atmospheric 

absorption 


A2: Low 
atmospheric 
absorption 


sorption, Ai, the planet appears larger than at a wave¬ 
length with lower atmospherie absorption, A 2 — 
beeause of the deeper transit depth due to the larger 
opaque atmospherie annulus. In partieular, a rela¬ 
tive flux-drop, ^(A), is assoeiated with an effeetive 
planetary radius, Rp{X) = Ri,\J ^(^. Transmis¬ 
sion speetroseopy mainly probes low-pressure lev¬ 
els; therefore, the mass eneompassed in the sphere 
of radius Rp{\) (Equation |3.2| ) is a good proxy for 
the planetary mass. 


Figure 3-1: Transit-depth variations, ^(A), indueed by the wavelength-dependent opaeity 
of a transiting planet’s atmosphere. The stellar disk and the planet are not resolved; the flux 
variation of a point souree is observed. 
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The atmospheric properties of a planet directly determine its effective radius in the 
following manner: 

poo 

TrRl{X) = tt[Rp^o +heffiX)f = / 27rr (l - dr, (3.3) 

Jo 

where Rp^, heff{X), and are respectively a planetary radius of reference[^ the ef¬ 

fective atmosphere height, and the planet’s transmittance at radius r (Figure [3^ . r(r, A) is 
the slant-path optical depth defined as 


/ -^oo 

r(r,A) = 2/ ^ni(r') x [T(r'),p(r'), A] dx, 


(3.4) 


where r' = \/r‘^ + (Figure [3^ B), and nj(r') and Ui [T(r'),p(r'), A] are the number 
density and the extinction cross section of the atmospheric component at the radial 


distance r' (Seager 20101. A planet’s atmospheric properties [ni{z),T{z), andp(^)] are 
embedded in its transmission spectrum through r(r, A) (Equations|3.3|and3.4). 



annulus area — 27 rrdr 


A 

dr = dz 


Stellar surface 



C 


Transmittance 


Figure 3-2: Basics of a planet’s transmission spectrum (planetary atmosphere scaled up to 
enhance visibility). (A) In-transit geometry as viewed by an observer presenting the areas 
of the atmospheric annulii affecting the transmission spectrum. (B) Side-view showing 
the flux transmitted through an atmospheric annulus of radius r. (C) Transmittance as a 
function of the radius at wavelengths with high and low atmospheric absorption—Ai (solid 
lines) and A 2 (dash-dotted lines), respectively. Due to higher atmospheric absorption at 
Ai, the planet will appear larger than it does at A 2 , because of the more-extended opaque 
atmospheric annulus [heff{Xi)) > /fe//(A 2 )] that translates into an additional flux drop. 


'a planetary radius of reference is any radial distance at which the body is optically thick in limb-looking 
over all the spectral band of interest, i.e., T(i?p qj A) ^ 1 VA S the band. 
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The integral in Equation |3.3| ean be split at the radius of referenee (beeause the planet 
is opaque at all A at smaller radii), and thus Equation |3.3| beeomes 

[-Rp,o + ^e//(A)] = + (3.5) 

noo I 

c = 2 (1 + r/)(l — dr/, ^ 

Jo ]_ 

y = z/Rpfl, 

leading direetly to the expression of the effeetive atmosphere height 

^e//(A) = Rpfl ^—1 + Vl + • (3.6) 

The embedded atmospherie information ean be straightforwardly aeeessed for most opti- 
eally aetive wavelength ranges using 


^e//(A) — i?p,o-B(7EM + In Aa)) 


(3.7) 


where 'fEM is the Euler-Maseheroni eonstant (Euler 17401 —'^em = linin-s.+oo YJJk=i I ~ 
Inn 0.57722 (see Seetion [T^ . In the above equation, B is a multiple of the dimension¬ 
less seale height and A\ is an extended slant-path optieal depth at referenee radius. The 
exaet formulation of B and Ax depends on the extinetion proeesses affeeting the transmis¬ 
sion speetrum at A (Table [3T| ). As an example, for Rayleigh seattering 


B = 


H 


R. 


and 


'P,0 




(3.8) 

(3.9) 


where nsc,o and crsc(A) are the number density at Rp^ and the eross-seetion of the seatter- 
ers. Coneeptually, Equation |3.7| tells us the height where the atmosphere ean be eonsidered 
transparent in limb-looking. As an example, if Ax is 10^ then the atmosphere ean be eon¬ 
sidered transparent at 9 (~ ln(10"‘)) seale heights above the referenee radius. 


^Edited from the version published in de Wit and Seager ( 2013| to correct Equation |3.6| and have consis¬ 
tent definitions of c in Equations|3.5|and 3.16 based on Stephen Messenger’s comments. 
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Interestingly, we also show that Equation|3.7|ean be rewritten as 


heffW = {z ■ t{z,X) = e , 


(3.10) 


meaning that the slant-path optieal depth at the effeetive height is the eonstant e zem inde¬ 
pendently of A (see Figure [3^ C). In other words, r(he//(A), A) ^ Teq = e ZEM _ Therefore 




lim n, I I ( 

n—>-1-00 J- 


-1/fc 


0.56146, 


(3.11) 


k=l 


whieh analytieally extends previous numerieal observations that Teq ~ 0.56 in some eases 
(jLeeavelier Des Etangs et aL| 2008|). Equations|3.7|and[3.10|ean be rewritten as 


Rp{\) = Rpfi [1 + B {jem + In Aa)] and 
RpW = : T(r, A) = , 


(3.12) 

(3.13) 


providing a simple analytieal way to implement and validate transmission speetrum models 
(see respeetively, Seetion [3^2.4| and Appendix]^ and Appendix]^. 

Most importantly, Equation |3.7| shows the dependeney of a transmission speetrum on 
its key parameters: in partieular. Ax is dependent in unique ways on the seale height, the 
referenee pressure, the temperature, and the number densities of the main atmospherie eon- 
stituents (for a summary of these dependeneies, see Seetion [3.2.4 ), whieh lead to the mean 
moleeular mass. The uniqueness of these dependeneies enables, in theory, the indepen¬ 
dent retrieval of eaeh of these key parameters. Therefore, a planet’s mass ean, in theory, 
be eonstrained uniquely by transmission speetroseopy (Equation |3.2| ). In praetiee, the de¬ 
termination of an exoplanet’s mass and atmospherie properties using solely transmission 
speetroseopy requires a high signal-to-noise ratio (SNR) over an extended speetral eover- 
age (e.g., 0.2 to 5 /im) adequately sampled (A/AA A 300)—see Seetions [33| and |A3^ 

The next seetion introduees the analytieal derivations that identify the key parameters 
of an exoplanetary transmission speetrum and its dependeneies on these. By identifying 
those dependeneies, we aim to find out if independent eonstraints on an atmosphere’s tem¬ 
perature, eomposition, and pressure ean be obtained solely from transmission speetroseopy. 
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3.2 Demonstration: Key Parameters of and their Effects 


on a Transmission Spectrum 


Here, we derive analytieally the dependeneies of a transmission speetrum on its main 
parameters for different extinetion proeesses. In partieular, we show that the effeetive 
height takes the form i?p o-B(7£;M+ln ^a) for extinetion proeesses sueh as Rayleigh seatter- 
ing, eollision-indueed absorption (CIA), and moleeular absorption for most optieally aetive 
wavelength ranges. Rp^B is a multiple of the seale height and Ax is an extended slant-path 
optieal depth at Rp^. B and Ax summarize how a planet’s atmospherie properties are em¬ 
bedded in its transmission speetrum. In partieular, the formulations of B and Ax identify 
the key parameters behind a planet’s transmission speetrum. The formulations of B and Ax 
(i.e., the way the atmospherie properties are embedded by transmission speetroseopy) de¬ 
pend on the extinetion proeesses. Therefore, we first approaeh in detail the ease of Rayleigh 
seattering—or any proeesses with an extinetion eross-seetion independent of p. Then we 
extend our demonstration to other proeesses sueh as CIA and moleeular absorption. 

For the eoming demonstrations, we will use the following assumptions: (al) the extent 
of the optieally aetive atmosphere is small eompared to the planetary radius {z Rp{X)), 
(a2) the atmosphere ean be assumed isothermal {dzT{Rp{X)) ~ 0), and (a3) the atmosphere 
ean be assumed isoeompositional, dzXi{Rp{X)) ~ 0 (where Xi is the mixing ratio of the 

atmospherie eonstituent). For a later generalization (see Appendix]^, we speeify for 
eaeh ease at whieh step these assumptions are used. 

3.2.1 Dependency for Rayleigh Scattering 

For extinetion proeesses like Rayleigh seattering, the eross seetion is independent of the 
pressure—i.e.,crA ^ fx{p)- Using the assumptions al, a2, and a3, the slant-path optieal 
depth (Equation|3.4[) ean be formulated as 



(3.14) 
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where the last term eomes from the integral over dx ([Fortney 


2005 j), or as, 


T 


{y,X)^Axe where < 


y = z/Rpfl 

v4a = Yji , (3.15) 

B = H/ Rpfi 


i.e., y and B are the dimensionless altitude and atmospherie seale height, respeetively, and 
Ax is the slant-path optieal depth at Rp^o —we reeall that the referenee radius is any radial 
distanee at whieh the body is optieally thiek in limb-looking over all the speetral band of 
interest. Therefore, Equation|3.5|ean be rewritten as 


/*00 , . 

ylffW + ^VeffW =c= 2 {1 + y) dy. 


(3.16) 


By solving Equation, 3.16 one finds j/e//(A) = —l + x/l + c. The integral in Equation |3. 16| 
evaluated analytieally over r is 


c 

2 


fO 


(1-0 (-|) dr 


'Ax 

B ^ Be-^ 


B^ In ^ B^ In 

Ax Ax 


dr. 


IA 


\ 


T T 


An evaluation of the integral of eaeh term of Equation|3.18| leads to 


rO n 

/ - 7 dr = -Blnr|»,. 


r-0 


Be- 


dr = BEAR'^ 




Jax r 

rO B^ In 
4, r 


dr = -B'AnAxlnrA^ + 0.55^ lnV|^ , and (3.21) 


dr = B‘^\nAxEi(T) 


B^ 


In re 


(3.17) 

(3.18) 


(3.19) 

(3.20) 


dr, (3.22) 


where Ei{x) is the exponential integral. The integral remaining in Equation |3.22| is equal 

to [r 3^3(1,1,1; 2, 2, 2; -r) - 0.5 ln(r) (ln(r) + 2r(0, r) + 2^em)] 

where pEq{ai, ...,ap; 6i,..., bq] z) is the generalized hypergeometrie funetion and is the 
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Euler-Mascheroni constant ( |Euler[ 17401. This new insight for the transmission spectrum 
equations (Equations |3.18p?22| ) allows further developments of Equation |3.17| using the 
following series expansions: 


• Ei{x)\^^^ = \nx + 'jEM - x + 0{x^). 

• [t 3 ^ 3 ( 1 ,1,1; 2, 2, 2; -x) - 0.5 ln(a;) (ln(a;) + 2r(0, x) + 2'yEM)] = ^+0(x). 

• [r 3 ^ 3 ( 1 ,1,1; 2, 2, 2; -x) - 0.5 ln(x) (ln(x) + 2r(0, x) + 2^em)] = ^21^ + 

0(x-5). 

The use of this last series expansion is appropriate if x > 5, i.e., in optically-aetive spectral 
bands where Aa > 5 beeause absorbers/diffusers prevent the light transmission at Rp^. 

In aetive speetral bands, we ean rewrite Equation |3.18| as 

^ = —B InO + i? In Aa + i?(ln0 + Jem) — BEi{Ax) — In Aa InO + B^ In^ 74 a + 0.5i?^ In^ 0 

-0.5^2 In^ Aa + 52 In 74A(ln 0 + jem) - In AxEi{Ax) - 0.5B^ In^ 0 + ^ 

= {jem + In^A - E,{Ax)) {B + B^ IiiAa) + B^ ( 3 . 23 ) 

Because Ei{Ax) 1 and jem 3> B(B 1 because of assumption al), we ean 
finally rewrite Equation |3.23| as 

c = (51n74A)^ + 2(1 + i?7^M).Bbi74A + (2S7EM)- (3.24) 

The last step is to use Bjem <C 1 to write 

1 + c ~ {Bjem + B\nAx + l)"^. (3.25) 

Therefore, we obtain the following solution to Equation |3.16| 

VeffW - S( 7 £;m TIhAa) and (3.26) 

riveff) ~ ^ 0.5615 (3.27) 
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—note that a first order approximation on c leads to |/e//(A) = B In Ax and T^y^ff) = 1. 

Equations |3.26| and 3.27 summarize the way in whieh a planet’s atmospherie proper¬ 
ties (ui, T, and p) are embedded in its transmission speetrum (Equation |3.3[ ). Coneeptually, 
Equation |3 .261 tells us the height (expressed in planetary radius) where the atmosphere ean 
be eonsidered transparent. As an example, if Ax is 10^ then the atmosphere ean be eonsid- 
ered transparent at 9 (~ In(10^)) seale heights above the referenee radius. Equation |3.27| 
shows that the slant-path optieal depth at the effeetive height is the eonstant (Eig- 

3-2[ panel C)—this extends previous numerieal observations that 


ure 


' eq ~ 


0.56 in some 


eases ( [Leeavelier Des Etangs et al.H2008| ). 

The appropriateness of Equations |3.26| and |3.27| is emphasized in Eigure [3^ that shows 
the relative deviation on the effeetive height between numerieal integration of Equation |3.16| 
and our analytieal solution (Equation |3.26| ). Eor Earth {B Ri 0.1%), the relative errors in the 
aetive speetral bands will be below 0.1% whieh eorresponds to an uneertainty on h^ffiX) 
and Rp{X) below 10 meters—whieh is well below the observational preeision we ever hope 
to aehieve. 


3.2.2 Dependency for Collision-Induced Absorption 


We extend here the results of Seetion [3.2.l| to extinetion proeesses with a oc like CIA. 
Eor sueh proeesses, Equation [3.15| ean be rewritten as 


T 


{y, A) ~ Axe where < 


y = 
Ax = 
B = 


^ / ^p,o 

and (3.28) 


/+! 


H 




is the temperature- and species-dependent absorption coefficient. For example, for 
CIA I = 1 and the ai^i^T = Ki{T)nl where Ki(T) depends solely on the temperature (Bo 


rysow[ 2002). Now that we obtain the same form for r as in Seetion [3 .2.1 1 (Equation |3 .1 5| ) 
we ean apply the same derivation leading to Equations |3.26| and |3.27[ Note that B and Ax 
have different formulations from the ones derived in Seetion [3.2.1] —although the general 
formulation of Equation|3.28|for proeesses with a oc eneompasses the ease of Rayleigh 
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logio(B) 

Figure 3-3: Relative errors on the effective atmosphere height (heff{X)) using the analytic 
solution to Equation |3.3| (or Equation |3.16[ ) for extinction processes with a cross-section 
independent of T and p (such as Rayleigh scattering). Relative errors for the zeroth-order 
derivation heff{X) = {z : t{z,X) = 1} correspond to the dotted lines, and the ones for our 
derivation heff{X) = {z : t{z, X) = are the thick lines. 
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scattering (/ = 0 ). 


3.2.3 Dependency for Molecular Absorption 


We extend here the results of Seetion 3.2.1 to moleeular absorption. We apply the same 
strategy as for CIA by showing that the slant-path optieal depth (Equation |3.4[ ) ean be for¬ 
mulated as r(|/, A) = —around heff{\) and for most A—and then relate to the 


derivation in Seetion 3.2.1 For moleeular absorption, the eross-seetion ean be expressed as 


cr,(A,T,p) = ^i,,(T)/„(A - A,,„T,p), 


(3.29) 


where Si^j and fij are the intensity and the line profile of the line of the atmospherie 
speeies. Eaeh quantity ean be approximated by 


nj' 


^ 5,,, (T,e/) and 


m=0 


Xi,j,T,p) ^ Aij{X,T) 


p + aij{X,T) 

p^ + h,{X,Ty 


(3.30) 

(3.31) 


where = 3 is suffieient to interpolate the line intensity dependeney to T (Gamaehe 


et al. 1990). Aij, aij, and bij are the new parameters we introduee to model the variation 
of the line with p at fixed {A,T} (Figure [?4 , panel B)—as an example, at low pressure 
fi,j{X - Xij,T,p) ^ Aij{X,T)aij{X,T)/hij{X,T) where Aij{X,T)aij{X,T)/bij{X,T) 
is the amplitude of the Doppler profile of at {A, T}. The seeond term of Equation 3.31| is a 
dimensionless rational funetion with a zero, —aij, and a pair of eomplex eonjugate poles, 
—details on rational funetions and their properties are provided in Appendix]^ The 
positions of the zero and the poles in C induee four regimes of speeifie dependeney of fij 
on T and p (Figure [3^ C). 

1. Doppler regime: while p < aij, fij is independent of p (i.e., fij ^ Aijaij/bij) 
beeause neither the zero nor the poles are aetivated. In terms of distanee to the 
line eenter (uj), the Doppler regime dominates when {u — uj) < 7 ^, where 7 ^ = 
{u : dllnfvli, = O} (/y and 7 y are the Voigt profile and its FWMH, respeetively). 
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2. Voigt-to-Doppler transition regime: while p'^ < bij, fij oc because only the 
zero is activated (i.e., p > aij). In particular, fij = Aij{p + aij)/hij. In terms of 
distance to the line center, this regime dominates when {u — Uj) < 'yv 

3. Voigt regime: while ~ bij and p > atj, fij behaves as the rational fraction 
introduced in Equation |3 .311 because the zero and the poles are activated. In terms of 
distance to the line center, this regime dominates when (z/ — i/j) ~ yy. 

4. Lorentzian regime: while and p 3> aij, fij oc p~^ because one zero and 

two poles are activated. In terms of distance to the line center, this regime dominates 
when {u — uj) > yy. 


While ttij and bij govern the regime of the line profile, Aij models its variation with 
temperature; Aij oc T"', where w is the broadening exponent. In the Lorentzian regime, w 
ranges from 0.4 to 0.75, while it is -0.5 in the Doppler regime ( |Seager[ [20I0| ). 

Using Equations |3.30| and |3.3I| and the ideal gas law, the absorption coefficient can be 
formulated as 


a{X,T,p) 


'^ni{T,p)ai{X,T,p) 

i 

VP -y rp- 

i m=0 


(332) 


y2s,,(T„f)Aij{\,T) 


p + at i(X, T)_ _ 

p2 + bij{X,T) 


Eor most A, the extinction process is dominated by one line. Therefore, the overall 
dependency of the absorption coefficient can be formulated as 


«(A, T, p) = 

where k = {A, T, and A^, a^, and 6^ are model parameters introduced to fit the absorp¬ 
tion coefficient variation in the T — p — X space with A being fixed. All of these parameters 
are known a priori from quantum physics and/or lab measurements ( [Rothman et al.[[2009| ). 
In particular, for most wavelengths, a^, and 6^ are XjS'jj(A)(T)Ajj(A)(A), aij(A)(A), 
and j(a)(A), respectively, where j(A) refers to the line that dominates at A. We show 
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,200K^^^ i,j,v-v ,200K^^^ T [Kl i.j.v-v ,1200K 

’j’ i/„x’ L''J ’j’ i/„\’ 



Figure 3-4: A line profile (fy) depends on the pressure (p) as a rational function at fixed 
temperature (T) and frequency (z/). (A) Dependency of fi, on T and p, at fixed z/, that 
shows the four domains of different dependency regimes of f,y on p whose boundaries are 
T-dependent (gray dot-dash lines). The black dot-dash lines represents the position of the 
slices in the {T — p — f} space used to highlight that f^, behaves as a rational function 
of p, at {T, z/} fixed. Planels B and C present the slices at T = 1200K and T = 200K, 
respectively. These slices show that f^, behaves as a rational function of p with a zero and 
a pair of conjugated zeros (Figure [A^ . In particular, the absolute value of the zero is less 
than the poles’, as underscored by the sequential transition from the following dependency 
regimes oc cx p^, oc and oc p~^ with increasing p —the dotted and the dashed lines 
represent a slope of 1 and 2, respectively. (D) Dependency of the absorption coefficient ( 0 ,^) 
on p, dXT = 200K. The exponent of the ai, dependency on p increases by one compared 
to dependency on p. The exponent increase by one because of aj^’s additional zero at 
p = 0 that originates from the number density. 
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in Figure[3^D that the absorption coefficient behaves like fij (Equation |3.31[ ) but with an 
additional zero at p = 0, which originates from the number density. 

Using the assumptions al, a2, and a3, Equation [34| becomes 


. OA rP^+^-PA 

t{z, A) = 2A«; / ^ dx, where 

Jo P + 


p = poexp{—z'/H) 


z' = z\z,x) 


(3.35) 


‘iRpfi 


+ z 


Transmission spectroscopy probes a limited range of atmospheric layers at each wave¬ 


length (see, e.g., Eigure 9-12(b) from |Elachi and van Zyl[|2006] ). Therefore, only a limited 
part of the dependency on T — p of the dominant line at A is recorded. We use this property 
to extend our demonstration assuming that each wavelength records only one regime of 
dependency (Eigure [3^ . By doing so, we show that the slant-path optical depth (Equa¬ 
tion 3.4) can be approached by r(p, A) = A\e~yl^ in the probed atmospheric layers (i.e.. 


around heff{X ))—and for most A. This approach extends the use of Equations |3.26| and 
3.27 to molecular absorption based on the derivation performed in Section[3.2.1[ 


Doppler regime 

When the dominant line at A behaves as a Doppler profile (p a^) in the atmospheric 
layers probed around /ie//(A), the absorption coefficient depends on the pressure only 
through the number density (i.e., alike the Rayleigh-scattering case, Section |T2.1| ). There¬ 
fore, Equation |3.35| can be rewritten as 

t{z, A) = (3.36) 


Voigt-to-Doppler transition regime 

In this regime, p^ <b^, therefore, Equation|3.35|can be rewritten as 


t{z,X) = ^y2'KRpfiHA^-Apfye 

Ou~ 



Po c-RH 


(3.37) 


This regime encompasses the three following subregimes: p a^, p ^ a,^, and p S> a^. 
The first subregime corresponds to the Doppler regime (Equation |3.36[ ). Equation [3.35| is 
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rewritten for the seeond and third subregimes, respeetively, as 


t(z,X) = {pqc and 

V 


(3.38) 

(3.39) 


Voigt regime 


This regime refers to the general formulation of the problem, i.e., when the Doppler and the 
Lorentzian behaviours affeet the line profile with comparable magnitudes. This formulation 
does not simplify the integral in Equation 3.35 Therefore, we rewrite Equation|3.35|for the 


transition case—i.e.,p^ ~ 6^ 3> a 


t{z,\) = ^^2^lRpflHA^ 


VK 


Poe 


-z/H 


(3.40) 


Lorentzian regime 

In this regime, > —bij and p S> —aiy; therefore, the absorption coefficient is mostly 
pressure-independent (Eigure [3^ panel D). As a result, such a regime is not expected to 
be recorded, under the assumptions a2 and a3. 


3.2.4 Summary and Discussion 

Now that we have gone through all the different extinction processes, we can find out 
what formulation of the effective height is generally true and identify the key parameters 
behind a planet’s transmission spectrum. We demonstrate that the slant-path optical depth 
(Equation|3.4[) is of the form 


A) = Axe 


(3.41) 
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for most A under the assumptions al, a2, and a3. As a result, the effective atmospheric 
height can be expressed as 


^e// — RpflBipfEM + In Aa). (3.42) 

We show the appropriateness of our demonstration in Figure [3^ using the numerical sim¬ 
ulation of the transmission spectrurrj^ of an Earth-sized planet with a isothermal and iso- 
compositional atmosphere—same abundances as at Earth’s surface. It confirms the pivotal 
role of 7 £;m in transmission spectroscopy as 


T P 

' eg ~ ^ 


(3.43) 


for a significant fraction of the active spectral bins (~99%). In addition, we note that a 
large fraction of the active bins (~70%) recorded a {oc P^}-dependency of r. This means 
that transmission spectroscopy preferably records transitions from the Voigt regime 
to the Doppler regime; because 7 ^ is small, the spectral bins are more likely to be on 
a line wing, rather than close to the line center. The fact that transmission spectroscopy 
records preferably the Voigt-to-Doppler transition regime is important because this regime 
embeds independent information about the pressure (Equations |3.39| and |3.39| ). 

We show in Eigure [T^ that for 50% of the optically active bins, our formulation for heff 
is still adequate for Earth, i.e., for a planet with atmospheric temperature and mixing ratios 
strongly dependent on the altitude^ We emphasize in Eigure [3P that water is the main 
cause for the deviation of Teg’s distribution from ~ 0.56 shown in Eigure [3^ The removal 
of water from Earth’s atmosphere leads to Teq ~ 0.56 for 80% of the active bins. The effect 
of water on Teg’s distribution is due to the strong variation of the water mixing ratio with the 
altitude close to Earth’s surface. In particular, the water scale height is smaller than the local 


atmospheric scale height, which invalidates the demonstration in Section 3.2.3 Therefore, 
heff deviates from the formulation RpfiB{'yEM + IhAa). As discussed above, it is not 
because hgff deviates from RpflB{'yEM + In Aa) that the specificity of the dependencies of 


^For details on our transmission spectrum model please refer to Section 


4.3.2 


'^Steps towards a generalization of the derivations performed in Section 3.2 are introduced in Appendix[B] 
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Figure 3-5: Two-dimensional histogram of Teg (A) as a function of log]^o(^A) for an 
isothermal-isocomposition Earth. The color scale is linearly proportional to the count in 
each bin. ~99% of the active spectral bins have a r^g ~ 


102 









a transmission spectrum to its parameters is lost—we are just currently unable to provide 
the general derivation in such a general atmosphere. Note that the significant effect of 
the scale height of a specie—such as water—on a transmission spectrum is favorable for 
habitability assessment (see the paragraph at the section end). 


Insights into the Dependencies Ax and B record the dependency of a planet’s transmis¬ 
sion spectrum on atmospheric properties in ways that vary based on the extinction process 
and the regime recorded at A (Table|3]T]). While B is solely affected by H, Ax is affected 
in independent ways by the following parameters: 

• The scale height, H, affects Ax via for all extinction processes, which is a 
geometry factor from the light path. 

• The number densities of the main atmospheric absorbents, Ui, affects Ax propor¬ 
tionally for molecular absorption—recall A^po = ni,o5'jj(A)(T)y4jj(A)(A)—and for 
Rayleigh scattering, while their squares affect Ax for CIA. These different depen¬ 
dencies on Tii result in stronger constraints on n* when assessed primarily from CIA 
signal. 

• The reference pressure, po, affects Ax for molecular absorption in the Voigt-to- 
Doppler regime (Aa cx: Pq^ if p ~ and Aa cc pj if p S> a^). Fortunately, 
this regime is recorded by most spectral bins and, hence, allows constraining the 
atmospheric pressure. 


The temperature, T, affects A a in two ways: (i) through A^ due to the line-intensity 
dependence on T (via a third-order polynomial) and (ii) through the line-profile de¬ 
pendence on T (T’"-dependency of the line broadening). In addition, processes like 
CIA are also dependent on T, in ways that are known a priori from quantum me¬ 
chanics (e.g.,|Borysow[|200^. 


Although the temperature and the number densities are both embedded in A^, fortunately 
they affect A^ in different ways. A^ oc n* independently of A, while it depends on T in a 
A-dependent way through the line intensity and profile. Note again that those dependency 
are known—or can be known —a priori from quantum mechanics or lab measurements. 
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Figure 3-6: Two-dimensional histogram of Teg (A) as a function of logio(^A) for Earth— 
temperature-pressure-mixing ratio profiles from |Cox and Pilachowski| ( |2000 ). The color 
scale is linearly proportional to the count in each bin. ~50% of the active spectral bins 


have a r^q ~ 


exp 


-lEM 
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Figure 3-7: Two-dimensional histogram of Teg (A) as a function of logio(^A) for Earth 
with no water in its atmosphere—temperature-pressure-mixing ratio profiles from|Cox and 


Pilachowski (2000). The color scale is linearly proportional to the count in each bin. ~80% 


of the active spectral bins have a Teg ~ exp 


—lEM 
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In summary, the key parameters of a transmission speetrum—the atmospherie seale 
height, the number densities of the main absorbents, the temperature, and the pressure— 
affeet it in unique ways. Therefore, it is theoretieally possible to retrieve eaeh of these 
parameters independently. 


Table 3.1: Dependeney of Ax and B to the regimes reeorded 


Extinction process —Recorded regime 

Aa 

B 

Rayleigh scattering 

x/2'kRp^qH 

H 

^p,0 

Collision Induced Absorption 

PJRPI Y,iK.(T)nl 

H 

2i?p,o 

Molecular absorption 

Doppler regime (p < a A 

x/27rRpflH A^f^Po 

H 

Rp,o 

Voigt-to-Doppler transition regime {p^ < bR) 
if p ~ Ok 

if p > 

x/ TiRpfiH Ak^Po 

H 

1.5-Rp,o 

H 

2r?p,o 

Voigt regime 

if p'^ b^:$> al 

\/2T^RpflH A.^^j=pq 

H 

Rp,o 


Discussion The analytieal derivations performed in this seetion required using the three 
following assumptions: (al) the extent of the optieally aetive atmosphere is small eom- 
pared to the planetary radius (z <C Rp), (a2) the atmosphere ean be assumed isothermal 
[dzT{Rp{X)) ~ 0], and (a3) the atmosphere ean be assumed isoeompositional [(i 2 Xj(i?p(A)) ~ 
0]. We diseuss below how Seetion [T^ s eonelusions are mostly unaffeeted by the relaxation 
of the assumptions a2 and a3—al being justified. In partieular, we explain eoneeptually 
that the effeet of eaeh key parameter remains unique while relaxing these assumptions, 
meaning that MassSpec ean be applied to any exoplanet atmosphere, theoretieally. In prae- 
tiee, MassSpeRs applieation requires the transmission speetra to be of high-SNR, suffieient 
speetral resolution, and with an extended speetral eoverage. 

The analytieal derivations beeome more eomplex if assumptions a2 and a3 are relaxed. 
We provide in Appendix]^ the first steps towards a generalization of the previous 
derivations. Below, we deseribe eoneeptually how these relaxations affeet the previous 
results. 
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(i) A non-isothermal atmosphere translates primarily into an altitude-dependent scale 
height {H = H{z)] —therefore, it is required to model the planet’s atmosphere with a 2 ;- 
dependent scale height in atmospheric retrieval methods (see, e.g., Section [4.3.3| ). However, 
the temperature can still be self-consistently retrieved from a planet’s transmission spec¬ 
trum because it affects in its own way the slant-path optical depth profile, r(r, A), through 
the extinction cross section profile, cXi (T(r'),p(r'), A). 

For molecular absorption, T affects Uj (T(r'),p(r'), A) through the line intensities and 
profiles in ways that are known a priori from quantum mechanics and/or lab measurements. 
T affects mainly the line intensities through the total internal partition sum (TIPS) and the 
Boltzmann populations for molecular absorption. The TIPS describes the overall popula¬ 
tion of the molecule’s quantum states, and is solely dependent on the molecular structure 
(i.e., the species) and the local temperature. The effect of temperature on the TIPS can 


be appropriately approximated by a third order polynomial ( [Gamache et al.[ |1990[ ). On 
the other hand, the temperature effect on the population of a molecule’s individual state 
depends solely on the temperature and the energy of the state. Therefore, the extinction 
cross-section depends on the temperature in ways that vary with the wavelength, in oppo¬ 
sition to the TIPS. Finally, the temperature affects the line broadening as oc T"', where w 
is the broadening exponent and ranges from 0.4 to 0.75 in the Lorentzian regime (classical 
value: 0.5), while it is -0.5 in the Doppler regime ( |Seager[ [2010| ). 

The overall effect of a change in T cannot be compensated/mimicked by other atmo¬ 
spheric parameters because it is specific—and a priori known. For example, although a 
local decrease in temperature could be compensated at the zeroth order by an increase of 
the local number densities: 


• The increase in number densities required to compensate the change in the line inten¬ 
sities will be inconsistent with the increase required to mitigate the change in local 
pressure. While a molecule’s local number density will have to scale as a third order 
polynomial to compensate line intensity changes, it should scale as 1/T to compen¬ 
sate the change in pressure. 

• The compensation enabled by an increase in number densities is wavelength-independent 
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while the effect of T on the absorption coefficient is strongly wavelength-dependent. 


(ii) A non-isocompositional atmosphere translates into different number-density scale 
height for each component. Alike in the case of temperature changes, changes in compo¬ 
sition affect the slant-path in specific ways, preserving MassSpec’s capability to constrain 


independently the parameters of the mass equation (Equation 3.2). A change in number 
density is species specific and, hence, cannot be compensated by global atmospheric quan¬ 
tities such as the temperature, scale height, or reference pressure. 

In some case, molecules such as water may require the use of an individual scale height 
for proper retrieval. The use of species-specific scale heights will likely he of primary 
importance to assess the hahitahility of planets, as a significantly smaller scale height 
for the water number density could indicate the presence of water surface reservoir. For 
such applications, comparisons with the scale height of other molecules will be required. 
Molecules that are expected to present constant mixing ratio throughout the atmosphere 
could then be considered as independent markers to extract the pressure scale height. The 
best marker candidate is carbon dioxide because (i) it is a chemically-stable molecule that 
is usually well-mixed in a planet’s atmosphere and (ii) it presents numerous strong ab¬ 
sorption bands that enable CO 2 detection at low abundance (down to ~ O.lppm). In the 
extreme cases where water and carbon dioxide are depleted from a planet’s atmosphere 
due to condensation, nitrogen and/or hydrogen—which are known to be chemically sta¬ 
ble at temperatures below which CO 2 condensates—would be the dominant species and, 
hence, the primary marker for the pressure scale height. 


3.3 Proof of Concept with HD 189733b 

We demonstrate in Section [T^ that, in theory, the transmission spectrum of an exoplanet 
constrains independently its atmospheric properties (temperature, pressure, and composi¬ 
tion) and, hence, its mass. In practice, the determination of an exoplanet’s mass and at¬ 
mospheric properties using solely transmission spectroscopy requires a high SNR over an 
adequately-sampled extended spectral coverage (e.g., from 0.2 to 5 pm with A/AA A 300), 
which will be available with future facilities (see Section[4.3.4[). This means that with cur- 
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rent facilities, MassSpec cannot yield independent constraints on the temperature, compo¬ 
sition, and pressure of a planet’s atmosphere solely from its transmission spectrum. There¬ 
fore, complementary information is required for immediate applications—see Equation|3.2[ 


We show here that MassSpec can be applied immediately to some hot Jupiters. Their 
temperature can be determined by emission spectroscopy. Their mean molecular mass is 
known a priori (H/He-dominated atmosphere: p ~ 2.3). And, their scale height can be 
directly derived from transmission spectroscopy via their “Rayleigh-scattering slope”. In 
spectral bands where Rayleigh scattering dominates, the effective planetary radius (Equa- 
tion |3.12[ ) can be rewritten using Equations |3.9| and [T9| as 


aH = 


di?p(A) 

dlnA ’ 


(3.44) 


because the cross section is of the form crsc(A) = cro(A/Ao)“ with a = -4 (Eecavelier 


Des Etangs et al. 2008). In other words, the slope of the effective planetary radius in the 
In A-space relates to the local scale heighj^ 


Therefore, using Equations|3.2|and[3.44|the planet mass can be derived from 


Mp = 


AkT[Rp{\)f 

dlnA 


(3.45) 


MassSpec's estimate of HD 189733b’s mass is 1.15 Mjup, based on estimates of T (i^ 
1300 K), di?p(~ 0.8/im)/dlnA -920 km), and i?p(~ 0.8/im) (f» 1.21i?j„p) derived 
from its emission and transmission spectra (Madhusudhan and Seager[ 2009 [ Pont et al. 


2008 2013). MassSpec's estimate of HD 189733b’s mass is in excellent agreement with 


the mass derived from RV measurements (1.14±0.056Mj„p Wright et al. 2011) for this 
extensively observed Jovian exoplanet. 


^We refer to Figures lT^ and |B-2| for examples of Rayleigh-scattering slopes for isothermal and non- 
isothermal atmospheres, respectively. 
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3.4 Preliminary Discussions 


3.4.1 The Importance of Accurate Extinction Cross Section Databases 

Atmospheric retrieval methods based on spectroscopy—like MassSpec —require accurate 
extinction cross sections (Equation |3.4| ). Conceptually, atmospheric retrieval methods solve 
the inverse problem of determining the conditions (i.e., temperature, pressure, composi¬ 
tion) of the medium probed from its transmission spectrum knowing how the transmission 
of light is affected by a medium’s conditions. Hence, our capability for characterizing 
exoplanetary atmospheres dependents on the accuracy of our knowledge concerning the 
dependency of a medium’s optical properties on its conditions. That is why we advocate 
for devoting significant efforts to generate accurate extinction cross section databases that 
cover various atmospheric conditions (i.e., temperature, pressure, composition). 


3.4.2 Complementarity of MassSpec and RV 

Here we introduce MassSpec's sensitivity to the planetary system parameters and show 
MassSpec's complementarity with the RV method to yield planetary mass measurements. 
For that purpose, we derive scaling laws to identify the main dependency of each method’s 
signal on the properties of the observed planetary system, such as the planet’s semi-major 
axis. 


Intensity of the RV and Transmission Signals 

The signal targeted by the RV method, the RV shift (A'*), can be expressed as 




G 


sm? 


{Mp -\- M^)a 


(3.46) 


where M*, Mp, a, f, and e are the host star’s mass and the planet’s mass, orbital semi¬ 
major axis, inclination, and eccentricity ( Murray and Correlal 2010|. (We do not discuss 
here the effects of i and e.) On the other hand, the signal in transmission depends solely 
on the area of the opaque atmospheric annulus, 2nRp^oheff{^), and the host star spectral 


no 









radiance Bx(T^) (where Bx is the Planek funetion and T* the star’s effective temperature). 
Using Equation |3.7[ we ean write that heff{X) oc H —for aetive moleeular bands Ax > 10^ 
therefore heff{X) = nH with n > 6. H = kT/jig where T ean be approaehed by the 
planet’s equilibrium temperature at first order, 

T,, = T,{Rjaf-^ (/'(I - , (3.47) 


where i?* is the star’s radius and a, /', and are respeetively the semi-major axis, a pa¬ 
rameter for the heat redistribution in the planet’s atmosphere]^ and the Bond albedo (Seager 


20101. By rewriting the planet’s surfaee gravity as g = ATiGppRp/S {pp is the planetary 


density), we ean summarize how the system parameters affeet the transmission spectrum 
and the RV signals (Table [X2] ). 


Table 3.2: Dependeney of a signal intensity in transmission spectroseopy and RV to the 
system parameters. 


Signal 

Planetary parameters 

Stellar parameters 

Transmission speetrum 


Bx{T,)%Rki-^ 

RV shift 

Mpa-^-^ 

M-0.5 


Table |3.2| highlights that signals in transmission are more intense for low-density plan¬ 
ets and atmospheres and bright or large stars, while the RV shift is amplified by massive 
planets and low-mass stars. In partieular, bright and large stars inerease the atmospherie 
temperature, henee the atmospherie seale height, for a fixed planetary density, atmospherie 
eomposition, and semi-major axis. Note that MassSpec will also be relevant for gas giants, 
espeeially for those whose star’s aetivity prevents a mass measurement with RV (e.g., the 


hottest known planet, WASP-33b Collier Cameron et al. 2010). 


Sensitivity of a Transmission Spectrum SNR 

In order to derive the aetual sealing law of transmission speetrum SNR—and, henee, 

MassSpeRs sensitivity—one has to aeeount for the observational parameters. The over- 

= 1/4 if the heat deposited by the stellar radiation is uniformly distributed. /' = 2/3 for a tidally- 
locked planet without atmospheric advection. 


Ill 

















all significance of a signal in transmission scales as 


SNRst oc 




(3.48) 


where the numerator and the denominator relate, respectively, to the number of photons 
bloeked by the planet atmosphere and the uneertainty on the baseline. The latter eor- 
responds to the square root of the total number of photon emitted by the host star and 
eolleeted out-of-transit in the same spectral band (Poisson process), respeetively. t is the 
observation time, A and r] are the teleseope’s colleeting area and total optieal throughput, 
and d is the host star’s distance to Earth. Equation|3.48|ean be rewritten as 


SNRst \ (3.49) 


Eor a given planet (Rp^, Mp, p, a, and T fixed), Equation[3.49|can be rewritten as 


SNRst oc 


x/tApBxjT^) 

R^,d 


(3.50) 


Equation |3.50| highlights the interest of M dwarfs for obtaining high-SNR transmission 
speetrum (or transit light-curve) for fixed planet properties (e.g., T^q ~ 300 K when seareh- 
ing for habitable planets). Equation |3.50| shows that for a given planet, the signal signif- 
icanee scales as x/Bx(T^)/{R-kd). Eigure 3^ shows the ratio x/Bx(T^)/R-k —normalized 
for a Sun-like star—as a funetion of the stellar effeetive temperature (T* — i?* relation based 
Reid and Hawley], 2005; Torres et al. 2010). Eor stars with earlier speetral types than 


on 


M2V, the signifieanee is independent of the host-star type. However, the signifieanee in- 
ereases substantially for M dwarfs. We take advantage of this favorable properties of late M 
dwarfs to show in Chapter]^ that MassSpec will enable the eharaeterization of Earth-sized 
planets within the next decade, possibly leading to the identification of the first habitable 
planets. 


112 



























Figure 3-8: The significance of in-transit signals increases for planets transiting M dwarfs. 
The figure shows the ratio Bx{T^) /i?*—normalized for a Sun-like star—as a function of 
the stellar e ffective temperature (T* — i?* relation based on Reid and Hawley[[2005 ; Torres 


et al. 


2010). x/B\{T^)/Ri, scales as the overall significance of an in-transit signal such as a 


transmission spectrum, for fixed planetary properties (Equation |3.50[ ). For stars with earlier 
spectral types than M2V, the significance is independent of the host-star type. However, 
the significance increase significantly towards late M dwarfs. 


3.4.3 Possible Insights Into Planetary Interiors 


Mass and radius are not always sufficient to obtain insights into a planet’s interior. MassSpeRs 
simultaneous constraints on a planet’s atmospheric properties and average density may help 
to break this degeneracy, in some cases. A precision on a planet mass of 3 to 15%, com¬ 
bined with the planetary radius can yield the planetary average density and hence bulk 
composition. Even with a relatively low precision of 10 to 15%, it is possible to infer 


whether or not a planet is predominantly rocky or predominantly composed of H/He (Sea- 


ger et al.[[2007[[Fortney et am2007[ ). With a higher planet mass precision, large ranges of 
planetary compositions can be ruled out for high- and low-mass planets, possibly revealing 
classes of planets with densities intermediate between terrestrial and giant planets with no 
Solar System counterpart ( [Rogers and Seagerj [201 0[ [Kipping et al.[ [2013[ ). Typically the 
bulk density alone cannot break the planet interior composition degeneracy, especially for 
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planets of intermediate density. However, measurement of atmospheric species may add 
enough information to reduce some of the planet interior composition degeneracies—e.g., 
the rejection of H/He as the dominant atmospheric species constrains significantly the bulk 
composition, independently of the mass uncertainty. 

3.4.4 Enhancing Capabilities 

We show that transmission spectroscopy can constrain the mass and the atmospheric prop¬ 
erties of an exoplanet. Therefore, any dataset that can independently enhance the con¬ 
straints on the mass and/or the atmospheric properties of an exoplanet can enhance MassSpec's 
capabilities via a joint analysis. In particular, RV measurements and emission spectra 
would be great assets to provide complementary constraints on a target’s mass and its atmo¬ 
spheric properties, respectively. High-SNR emission spectra will be acquired in alternation 
with transmission spectra and are particularly sensitive to temperature. 
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Chapter 4 


Future Prospects 




In this Chapter, we assess the potential of the methods introduced in Chapters|^and|^when 
applied to observations from space-based facilitie^ We introduce short-term goals such 
as mapping for the first time an exoplanet’s atmosphere in 3D and weighting the hottest 


exoplanet known with the Spitzer Space Telescope {Spitzer, [Werner et al.j |2004| ) and the 
Hubble Space Telescope (HST), respectively. We discuss long-term goals such as obtain¬ 
ing time-dependent 3D maps of exoplanetary atmospheres and assessing the habitability 
of Earth-sized planets with the James Webb Space Telescope (JWST; launch date 2018, 
Clampin 2010) and EChO-class mission^ i.e. within the next decade. 


4.1 Telescopes Noise Models 


We introduce here the optical performance models used to assess the potential of our meth¬ 
ods with current and future telescopes. Conceptually, such models are mainly effected by 
(1) the collecting area, (2) the total optical throughpuj^ and (3) the instrument’s duty cycle 
for a given target. 


‘Work published in part in Science, Volume 342, Issue 6165, pp. 1473-1477 (2013), see 
Seager (2013|l and in Experimental Astronomy, see 


de Wit and 


Parmentier et al. (2014|l. Work performed in part to 


support observation proposals: two for the Spitzer Space Telescope —PI: de Wit, rejected, and PI: Lewis, 
accepted, see Lewis et al. (2013ai—and one for the Hubble Space Telescope —PI: Crossfield, under review. 

'Both methods can be applied to ground-based observations. 

^Exoplanet Characterisation Observatory (EChO) was a M3 mission candidate of the European Space 
Agency ( Tinetti et al.| 2012| ). While EChO was not funded, such a mission could still revolutionize our 
understanding of exoplanet atmospheres. Hence, we point out here the potential of our methods when applied 
to the observations of an EC/iO-class mission— ‘"EChO”, for brevity. 

^The total optical throughput of a telescope is the ratio of photons collected by its primary mirror to 
electrons read over a spectral resolution element of its detector. 
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HSTrWFC3 


We model the optical efficiency of HST and the Wide Field Camera 3 (WFC3) as follow. 
We use a telescope effective area of 4.5 square meter, a global optical efficiency of 0.14 
0.25 to 0.95 fim and of 0.20 from 1.00 to 1.75 /im—including optical throughput and 
detector quantum efficiency. We derive the duty cycle on a target-by-target basis to account 
adequately for the detector saturation—e.g., ~ 40% for the brightest target of our sample, 
WASP-33b. We validate our model with real data. In particular, our model predicts a RMS 
of 110 ppm for 103-sec bin of WASP-43, which is in agreement with the RMS of actual 
observation is ~ 100 ppm (data from GO 13467, PI Bean). 

We recall here that HST is in a low Earth orbit (LEO), hence it does not allow continu¬ 
ous monitoring for most targets, but rather ~ 48-min viewing windows. 


SpitzerflRAC 

To date, no model of Spitzer/lRAC's optical efficiency has been developed in the context 
of this thesis. The mapping capabilities introduced in Section [432] are based on scaling 
the RMS on previous observations with the square root of number of observation to be 
performed. 


JWST/NIRSpec 


NIRSpec (Rauscher et ah, 2007aI is the Near-Infrared Spectrograph for the James Webb 


Space Telescope (JWST, see |Clampml|20I0| ). The purpose of NIRSpec is to provide low (R 
= 100), medium (R = 1000), and high-resolution (R = 2700) spectroscopic observations 
from 0.6 to 5 p,m. We focus here on the medium resolution mode because MassSpec 
requires a sufficient spectral resolution. Eurthermore, a larger spectral dispersion enables 
the observation of brighter stars. However, the medium and high resolution modes require 
the necessity to use individually three grisms to obtain a full spectrum, hence, decreasing 
the effective observation time by three. 
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Figure 4-1: JWSTINIRSpec optical performance for the medium spectral resolution mode 
(R = 1000). The red line shows the total optical throughput for each spectral resolution 
element. The green line shows the flux fraction going to the brightest pixel of each spectral 
resolution element. 


JWST/NIRSpec Total Optical Throughput We determine the total optical throughput 


of JWST/NIRSpec according to Boker and Tumlinson (2010). We model the detector pixel 
efficiency as a plateau at the 75% level with a drop of 3% at the pixel edges (G. Cataldo, 
private communication). We integrate the spectrograph A-dependent PSF over the pixel 
grid to estimate the resolution element sensitivity. By doing so, we account for distortions 
for wavelengths under 1 pm using the 1 pm PSF-size. We present the throughput budget 
summary for the NIRSpec medium-resolution mode in Figure [47T| Figure [47T] also shows 
our estimate of the flux fraction going to the brightest pixel of each resolution element 
which we will use further to derive the saturation time of the detector. 

JWST/NIRSpec Noise Budget We derive the budget noise from the instrumental perfor¬ 
mance shown in Figure [47T| First, we estimate the electron flux produced on the brightest 
pixel of the detector for a given star. Then, we derive the saturation time from the electron 
flux and a pixel well capacity of 60,000 e“ to ensure linear response. We derive the vari¬ 
ance of a frame, (Tii7’(A), including a readout noise of 6 e“ px“^ rms and a dark current of 


0.03 e (px s) ^ and assuming the read mode to be MULTIACCUM-2xl (Boker and Tum- 


linson[ 2010). We derive the duty cycle based on a readtime of ~0.53 sec, which, together 


with crii7’(A), leads to the signal-to-noise ratio (SNR) for a given observation time, t. 
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JWST/NIRSpec Noise Model Comparison We find JWST/NIRSpec observations will 
be photo-noise dominated, unlike [Deming et al.| ( |2009[ ) that suggested that the intrapixel- 
sensitivity indueed flux variations (see Seetion [2.3.3[ ) would be the main souree of noise. 
Reeent estimates of NIRSpec’s deteetor performanee show that its pixel-phase deviates 
mainly from a plateau at the pixel edge with a relative drop below 2% due mainly by eross- 
talk (G. Cataldo, private eommunieation)—i.e., the information is mainly transferred, not 


lost. These variations are more than an order of magnitude less than presented in Deming 


et al. (2009). Furthermore, even for eurrent faeilities signifieantly affeeted by pixel-phase 


like the Spitzer Space Telescope, the observations are within 10 to 20% of the photon noise 
after systematie eorreetion. In addition, we take into aeeount the neeessity to ehange the 


grisms to obtain a full speetrum (Boker and Tumlinson 2010). For praetieal purposes 
ineluding stability and baseline follow-up, we eonsider that one transit is observed in a 
unique grism. For that reason we further sealed the SNR by 1 /a/S to inelude the inherent 
sharing of the integration time between the three ehannels of NIRSpec's medium-resolution 
module. 


EChO-Class Mission 


While EChO was not funded, sueh a mission eould still revolutionize our understanding of 
exoplanet atmospheres. For that reason, we ehose to point out the potential of our methods 
when applied to the observations of an EChO-cXd&s mission. In order to model a future 


EChO-clz-SS mission, we use EChO's noise model introdueed in Barstow et al. (2013al. 


In partieular, we use a teleseope effeetive area of 1.13 square meter, a deteetor quantum 
effleieney of 0.7, a duty eyele of 0.8, and an optieal throughput of 0.191 from 0.4 to 0.8 
pm, 0.284 from 0.8 to 1.5 pm, 0.278 from 1.5 to 2.5 pm, 0.378 from 2.5 to 5 pm, 0.418 
from 5 to 8.5 pm, 0.418 from 8.5 to 11 pm, 0.326 from 11 to 16 pm. 


Future-generation telescope 

We use a sealed-up version of EChO to model the 20-meter spaee-based teleseope. In 
partieular, we use a speetral resolution of 1000, a deteetor quantum effleieney of 0.7, a duty 
eyele of 0.8, and an optieal throughput of 0.4 from 0.4 to 16 pm. 
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4.2 Mapping Potential 


4.2.1 On the Benefits of Multi-Wavelength Observations 


The spectral range used for light-curve observations determines the layers probed in the tar¬ 
get’s atmosphere, as its optical depth is A-dependent. Hence, if the SNR is high enough it is 
possible to map exoplanets’ atmospheres in 3D with multi-wavelength observations. Such 
observations also mitigate the correlation between the system parameters and a planet’s 
brightness distribution. The system parameters are not wavelength-dependent, hence their 
effects on the occultation shape can better be disentangled from the effect of the wavelength- 
dependent brightness distribution. 

Applications of our method to observations in the visible can also provide us with 
insights into a planet’s cloud properties (e.g., Section |23] ). Recently, |Madhusudhan and 
Burrows| ( 2012[ ) provide a theoretical framework for interpreting geometrical albedos from 
phase curves, which is is indicative of the scattering and absorptive properties of the atmo¬ 
sphere. Our method could, hence, complement such frameworks with 2D constraints on 
the geometrical albedo. 


4.2.2 Towards the First 3D Map of an Exoplanet’s Atmosphere 


As a contribution to the accepted Spitzer proposal ID #10103 ( [Lewis et~aL| |2013a| ) we 
demonstrate that complementary observations of the hot Jupiters HD 189733b and HD 209458b 
can yield to the first 3D maps of exoplanetary atmospheres. In particular, we estimate that 
a total of nine and seven eclipses respectively in Spitzer/IRAC 3.6 and 4.5-pm channels 
would enable to map HD 189733b’s atmospheric layers probed at those wavelengths and 
twenty-one and sixteen eclipses of HD 209458b respectively at 3.6 and 4.5 pm would yield 
the same results. 

The end of the complementary observations is scheduled before 2015, hence we expect 
the first 3D maps of exoplanet atmospheres within a year. Note that their resolution will 
be limited both in longitude/latitude and in altitude because of a limited SNR and spectral 
resolution, respectively. 
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We demonstrate the feasibility of mapping in 3D HD189733b’s and HD209458b’s 


atmospheres as follow: 

1. We determine the eurrent RMS achievable on archived occultation observations of 
our targets in the channels of interest—200 ppm and 175 ppm per 1-min bin in 
IRAC’s 3.6 and 4.5-/im channels respectively for HD 189733b and 300 ppm and 250 
ppm for HD 209458b. These RMS estimates account also for the correlated noise. 


2. We estimate the amplitude of the deviation in occultation ingress/egress due to the 
targets’ brightness patterns. We expect deviation with an amplitude of ~ 180 ppm 
(see Figure|47^and compare with FigurejlHO]). 


3. We estimate the number of eclipses required to secure a lOa detection of these de¬ 
viations in ingress/egress. We choose to aim for a high significance in order to 
mitigate the correlation between the brightness distribution and the system param¬ 
eters (see Section [T2| ). We estimate the significance of the deviation as follow. The 
ingress/egress duration of our targets is 25 minutes, hence a total of Nb ~50 bins cap¬ 
ture an average deviation (Da) of half the deviation amplitude (Da ~90 ppm). We 
thus need to achieve a RMS of ~70 ppm per 1-min bin—i.e., 70 Da/a*y/Nb with 
cr = 10 the significance level chosen. To do so, we need a total of nine, ~ (200/70)^, 
and seven eclipses, ~ (175/70)^, at 3.6 and 4.5 /im for HD 189733b’s and twenty- 
one and sixteen—~ (300/70)^ and ~ (250/70)^, respectively—at 3.6 and 4.5 fim 


for HD 209458b. 

4. Finally, we demonstrate that the RMS requested above would be sufficient to retrieve 
meaningful constraints on the brightness distributions of our targets. For that pur¬ 
pose, we simulate synthetic light curves and show that we could retrieve the bright¬ 
ness maps used to produce the synthetic light curve. As an example, Figure|4^ 
shows the significance of deviations from a uniformly-bright disk that we expect to 
detect for both targets at 3.6 and 4.5 /rm. Figures |4-3a| and [4^3b show the simulated 
and retrieved brightness maps, respectively. Figure [4^ shows that the numbers of 
eclipses requested yield a SNR sufficient to retrieve a brightness map consistent with 
the synthetic map used to simulate the synthetic data. This demonstrates the feasibil- 
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Figure 4-2: Synthetic observations of HD 189733b in Spitzerl\K\C 4.5-/im channel for the 
brightness map shown in Figure [47^ The red line represents the best fit using a conven¬ 
tional model (i.e., assuming HD 189733b to be a uniformly-bright disk) and the green line 
represents the best fit using a non-uniform dayside brightness model. The deviation of the 
occultation shape from the occultation of a uniformly-bright disk is shown in the bottom 
pannel. 

ity of our proposal. We also investigate what constraints on the location and extent 
of bright “hot-spots” would be achievable, as those constraints provide insights into 
the underlying physical processes shaping these features. We show that the comple¬ 
mentary observations will enable: 

(a) Localization of the peak in the dayside brightness distribution with an uncer¬ 
tainty ~2-3° in longitude and ~6-10° in latitude. If the hottest regions of the 
planet were substantially shifted from the equatorial region (e.g. 45° in latitude) 
the precision with which we could localize the peak in the dayside brightness 
would be decreased to ~3.5-6.0° in longitude and ~ 12-20° in latitude. 

(b) A precision better than 17° on the latitudinal extent of the hot-spots and better 
than 8° on the longitudinal one, i.e., a precision relevant for comparisons with 
three-dimensional hot-Jupiter GCMs. 
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Figure 4-3: Feasibility of mapping HD 189733b and HD 209458b using SpitzerllKKC 
3.6 and 4.5/im channels, (a) shows one of the synthetic map used to simulate synthetic 
Sp/zter/IRAC’s datasets of HD 189733b (see Figure|4^ and to estimate the mapping ca¬ 
pabilities with seven eclispes. [(b)] retrieved brightness map that is consistent with [(a)] and 
hence demonstrate the feasibility of our proposal—i.e., of mapping in 3D HD 189733b’s 
and HD 209458b’s atmospheres. 
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4.2.3 3D Map of the Hottest Jupiters in the Visible 


As a contribution to Dr. 1. Crossfield’s HST proposal to map the hottest planetary atmo¬ 
spheres, we applied a methodology analogous to the one deseribed at the previous Section. 
The main differences with the methodology adopted for the Spitzer proposal are that we 
investigate the entire list of planets known to determine which would be the best targets to 
perform mapping and we use the HSTfWFC3 model introduced in Section [4T| Note that 
HSHs orbital eonfiguration implies that the observation of an oeeultation ingress or egress 
requires at least two HST orbits, one to perform the pointing and a seeond to perform the 
seienee—in some eases, it might be neeessary to perform an observation before and after 
the seienee in order to get the light-eurve baseline. 


Finding the Most Suitable Planets to Map in the Visible 


Finding the best targets to map with a given teleseope requires determining the number 
of occultations neeessary to reaeh a sufficient SNR on the oeeultation ingress/egress. We 
eonsidered the exoplanets listed in Wright et al. (20111 and Schneider et al. ( 2011| ) and 
proeeed for eaeh of them as follow: 


1. Determine the RMS per l-min bin over the HSTfWFC3 speetral range based on their 
host-star properties. 

2. Based on the target’s brightness temperature, we derive its eelipse depth. If bright¬ 
ness temperature is not yet eonstrained, we derive its equilibrium temperature (Equa- 
tion [3471 ). 

3. We derived the number of eelipses required to yield a ratio eelipse depth to RMS per 
l-min bin of 30, whieh ensure a sufficient preeision of the ingress/egress. We note 
that we use an extra faetor to aeeount for the dependenee on the planet’s brightness 
temperature of the brightness contrast ratio between a hot spot and the background 
temperature. Figure [4^ shows that the signal in oeeultation ingress/egress of colder 
planets is enhanced by an increased eontrast ratio, and henee partially eompensate 
for their lower eelipse depth in HSTfWFC3 speetral band. 
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TefftX] 

Figure 4-4: Dependence on the brightness temperature of the brightness contrast ratio be¬ 
tween a planetary hot spot and the background atmosphere in HSTfW¥C3 spectral band. 
The brightness contrast is shown for a modeled 10% temperature difference between the 
hot spot and the background. 
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Table 4.1: Ten most favorable exoplanets to map with HSTfW¥C3 


Target 

Number of eclipses 

WASP-33b 

8 

WASP-18b 

18 

WASP-76b 

18 

WASP-19b 

36 

WASP-12b 

40 

WASP-103b 

58 

WASP-77 Ab 

74 

WASP-82b 

93 

WASP-3b 

102 

WASP-14b 

120 


We present in Table |4.l| the ten best targets to perform exoplanet mapping with HSTfW¥C3. 
We sort these targets based on the number of eclipses required to secure a sufficient SNR 
to perform eclipse mapping with HSTfW¥C3. In Dr. I. Crossfield’s proposal, we pro¬ 
posed to target WASP-18b and WASP-19b for scientific considerations that are not part of 
the present work but also because: (1) WASP-33b orbits a Delta Scuti variable, hence its 
pulsations would mitigate the mapping performance—particularly in the context of HSTs 
non-continuous monitoring—and (2) we have currently very view observational constraints 
concerning WASP-76b. For future references, we also note that mapping WASP-12b would 
require a particular attention as it is distorted by its host (Li et al., 2010)—see Section[2^ 


3D Maps of WASP-18b and WASP-19b with HSTfWFCS 


As for [Lewis et al. (2013a), we demonstrate that the number of eclipses mentioned in Ta- 
ble |4.1| for WASP-18b and WASP- 19b would be sufficient to retrieve meaningful constraints 
on the brightness distributions of our targets. For that purpose, we simulate synthetic light 
curve—including partial phase curve—and show that we could retrieve the brightness map 
used to produce the synthetic light curve. We also show that these observations would 
enable: 


1. The localization of the peak in the dayside brightness distribution with an uncertainty 
~10° in longitude and latitude. 
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2. Constraints on the extent and the intensity distribution of the hot-spots with a preci¬ 
sion better than 20° on the extent in longitude and in latitude. 

The HST proposal partially described above is currently under review for HST cycle 

22 . 


4.2.4 Time Variability of Exoplanets’ Amospheres 3D Structures 


As a contribution to Parmentier et al. (2014), we investigate the mapping capability of an 
EChO-c\d&s mission based on the EChO's model introduced in Section lTTl We also discuss 
here JWST’s mapping capabilities. 

We show in Figure [43] the mapping potential of an EChO-class mission. In particu¬ 
lar, Figure |4^ shows the map retrieval of a synthetic brightness map of the hot Jupiter 
HD 189733b harboring a hot spot with a temperature contrast of 30% and located on the 
northern hemisphere. The retrieval assumes observations in a spectral bin of resolution 
~ 20—here, in the 4.3/im band of carbon dioxide. Such a hot spot could be formed by 
the presence of patchy clouds or chemical differences between the poles and the equator. 
With one secondary eclipse, EChO would detect the presence of latitudinal asymmetry in 
the planet’s brightness distribution. With ~ 10 and (resp. ~ 100) secondary eclipses, the 
temperature contrast will be measured with a precision of 300 K (resp. 100 K) and the 
latitudinal location of the hot-spot will be determined with a precision of 10° (resp. 3.5°). 

An EChO-c\ass mission would thus be able to yield 3D maps of exoplanets as such 
2D maps would be available at different optical depths, for the most favorable targets. In 
addition, we find that for EChO\ 

• Observation programs dedicated to 3D mapping in the infrared will require a mini¬ 
mum of 150 hours of observation. 


Observation programs aiming to constrain cloud properties will only provide longi¬ 
tudinal constraints (similarly to Demory et al. 2013| ) and will require a minimum of 
250 hours of phase curve observation. It will not be possible to constrain latitudinally 
the cloud coverage because the contrast of feature in reflection (clouds) is ~10 times 
less than in emission (hot or dark spots), in terms of relative brightness. 
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• EChO would be able to capture temporal variabilities in an planetary atmosphere 
if these variations have either large amplitudes (e.g., hot spot oscillation from one 
tropic to another) or large time scales (i.e., > 10 orbital periods to enable capturing 
with a sufficient SNR the extreme brightness distribution of the variations). 

JWST’s large collecting area yields to a significant reduction of eclipses required to 
perform 3D mapping. JWST would yield a 3D map with a single occultation for most 
targets for which an EChO-cl^&s mission could yield 3D mapping. Therefore, the obser¬ 
vations of consecutive occultations will provide unprecedented insights into the physics of 
those planet’s atmosphere in the form of constraints on the time variability of their 3D at¬ 
mospheric structures. However, JWST will not yield cloud maps because its spectral band 
do not cover the visible. Still, JWST could provide strong inference on cloud coverage via 
the detections of dark spots in the infrared—if those spots can be associated to large albedo 
using complementary facilities. 


4.2.5 From Maps to Insights into Atmospheric Physics 


The analysis method introduced in Chapter]^ constrains the brightness distribution of a tar¬ 
get, not its underlying physics. A way around that would be to use brightness models (see 
Section [Z3.2| ) with parameters carrying a physical meaning. As an example, |Heng and 


Workman (2014) derived key analytical formulations for the global structure of planetary 


atmospheres—under the shallow-water assumption—which are governed by up to three 
dimensionless parameters related to the atmospheric conditions. Therefore, the implemen¬ 
tation of such formulations in our mapping method could be used to map exoplanets while 
gaining additional insights into the physicals processes at play in the target’s atmosphere 
via the determination of meaningful brightness-model parameters. 


'^The discussion about such an implementation was initiated by Dr. K. Heng at the Exoclimes III confer¬ 
ence. 
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Figure 4-5: Mapping potential of an £'C/?0-class mission. Synthetic distribution of bright¬ 
ness temperature for the hot Jupiter HD 189733b (“Simulated map”) and the retrieved map 
as a function of the number of occulation (1, 10, and 100) observed in a spectral bin of 
resolution 20—here, in the 4.3/im band of carbon dioxide. 


128 































4.3 Mass Spec’s Potential 


MassSpec’s applications require a high SNR over an extended adequately-sampled spee- 
tral eoverage (e.g., from 0.2 to 5 pm with A/AA A 300) to eonstrain independently an 
exoplanet’s mass and atmospherie properties. We assess in this Seetion the full potential 
of MassSpec using future spaee-based faeilities. But first, we show below that although 
MassSpec’s applieations are eurrently limited (see Seetion[3^ it eould already yield sig- 
nifieant seienee using HSTfW¥C3, sueh as the mass of the hottest planet known. 


4.3.1 Short-Term: Weighing the Hottest Planet Known 


MassSpec’s applieation is eurrently limited to giant planets (Seetion [33] ). Although most 
giant planets’s masses ean be determined from RV measurements, for some their host-stars’ 
aetivities prevent a mass measurement (e.g., the hottest planet known, WASP-33b, |Collier 


Cameron et al. 20101. We show here how a reasonable amount of HST’s time ean yield a 
mass measurement for WASP-33b. For that purpose, we proeeed as follow: 

1. We estimate the RMS on the transit depth in 0.1 /xm-bins in the HST AVFC3 UV/Visible 


ehannel using the noise model introdueed in Seetion 4.1 We assume a eontribution 
of eorrelated noise post-eorreetion for WASP-33’s variability of 259^ and estimate 
the preeision on the transit depth to ~45 ppm per bin. 

2. We simulate WASP-33b’s synthetie transmission speetrum assuming a pure seatter- 
ing slope (Figure [4^ beeause no moleeular absorption is expeeted at the target’s 
extreme temperature (T = 3200 K and R ~ 1.48 Rjup, Smith et aL| 20111. We 
model the speetrum based on Equation|3.45|for a synthetie WASP-33b of 3 Jupiter’s 


mass . 


3. We perturb the theoretieal transmission speetrum using our estimate of the uneer- 
tainty on the transit depth in eaeh speetral bin. 


^We refer to “MassSpec's full potential” for its capability to retrieve the mass and the atmospheric prop¬ 
erties of a planet solely from transmission spectroscopy. 

®We consider that a simultaneous and continuous monitoring of WASP-33 using multi-channel ground- 
based instruments enable to correct for WASP-33 activity over a transit (see, e.g., von Essen et al. 2014[l. 
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4. We fit the synthetie signal to retrieve the seattering slope and derive the relative un- 
eertainty on the WASP-33b’s mass that would originate from the uneertainty on the 
seattering slope. 


.495 
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Figure 4-6: Weighing the hottest planet known with HST. Synthetie transmission speetrum 
of WASP-33b (blaek dots) and their uneertainties (grey bars) in the HSTfWVCi UV/Visible 
band together with the best fit (green). 

We find that the uneertainty on WASP-33b’s mass due to the seattering slope uneertainty 
eould be as low as 6% using five transits observed with HSTTNVd. This preeision is in 
strong eontrast with the eurrent eonstraints on WASP-33b’s mass from RV measurements. 
However, we note that the uneertainty on WASP-33b’s mass would then be primarily af- 
feeted by our uneertainty on its temperature—i.e., on its temperature profile as a vast range 
of pressure levels are probed between 0.2 and 1 /ini[^ In partieular, when deriving the mass 
of a gas giant solely from its seattering slope Equation |3.45[ the relative uneertainty on its 
mass, ean be approaehed at first order as follow 



(4.1) 


^In practice, the Rayleight-scattering “slope” would also provide the target’s atmospheric temperature 
profile (Appendix|B.2|and Figure[B^. 
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where ax is the uneertainty on the parameter X. We note that for WASP-33b the relative 
uneertainty on Rp and /r are negligible eompare to the uneertainty on the temperature profile 
^ < 0.15. Henee, we are eonfident in MassSpec’s eapaeity to yield the mass of the hottest 
planet known with a relative uneertainty of < 20% together with refined eonstraints on its 
temperature profile with less than 20 HST orbits. 


4.3.2 Transmission Spectrum Model 

We model high resolutiouj^ (i? = A/AA = 10®) transmission speetra following Equa- 
tion |3.3[ (i) we ealeulate the slant-path optieal depth for different altitudes and (ii) inte¬ 
grated the eontribution of eaeh projeeted atmospherie annulus to the overall flux drop— 
27 rrdr(l — exp[— r(r, A)]) (Figure [3^ . The extinetion eross seetion aeeounts for moleeu- 
lar absorption, eollision-indueed absorption (CIA), and Rayleigh seattering. We eompute 
the monoehromatie moleeular absorption eross seetions from the HITRAN 2008 database 
( [Rothman et al.[ [2009] ) and model the line profile aeeording to [Liu et ^ ( |2001[ ) to inerease 


the eomputational speed. We use opaeity tables from Borysow (2002) for H 2 —H 2 CIA. We 


determine the Rayleigh-seattering eross seetion, aRp in egs units, from 




ni\^ \Ni{\f + 2 


(4.2) 


where A is the wavelength and n*, Aj(A), and Fi{X) are the number density, the refraetive 
index, and the King eorreetion faetor for the depolarization of the i*^ atmospherie speeies. 


In partieular, we use the refraetive indiees of N2, CO2, CO, CH4, N2O from Sneep and 


Ubaehs(2005 


is necessary to model the radiative transfer process at high resolution to approach adequately the 
absorption lines and their effects on the stellar light. Once the simulated transmission spectrum is modeled it 

can be binned down to the facility’s spectral resolution for comparison. _ 

“^We observe discrepancies between the measured data and the functional forms proposed in Sneep and 


Ubaehs (2005|l for the refractive indices of CO 2 and CO—their equations (13) and (17). Therefore, we use 


the following corrected forms: 


nco ^ - 1 5799.25 , 120.05 

1.1427 X 103 “ (128908.9)2 - A-2 (89223.8)2 - A-2 

5.3334 4.3244 ^ 0.1218145 

(75037.5)2 _ + (67837.7)2 _ ^-2 + (2418.136)2 - A-2 ’ 
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A transmission spectrum simulation requires the computation of millions of absorption 
lines for each atmospheric species and numerous T — p conditions. In addition, we develop 
our retrieval method in a Bayesian framework which requires a large number of transmis¬ 
sion model runs to converge. Therefore, we compute the extinction cross section for each 
component of HITRAN as a function of A, T and p to interpolate later for the required 
conditions. In particular, we generate the extinction cross section 4-D array 

• at a spectral resolution of 10® for 0.4/im to 250/im, 

• for 17 pressure values spread in the logioP[Pa] space from 7 to -3 with a higher 
density around 4—because most of the information is recorded around 1 mbar. 


for 12 temperature values homogeneously spread from 150 K to 700 K. 


We validate our extinction cross section model with Benneke and Seag^ (2012). In 
addition, we validate our transmission spectrum model comparing a synthetic Earth trans¬ 
mission spectrum with |Kaltenegger and Traub (2009) and Hu et al.| (2013). Except for a 
3.2/im-water signaturep^ our simulation is consistent with Kaltenegger and Traub (2009) 
and Hu et al.| ( |2013| ); for that reason, we consider our transmission model appropriate for 
the present study. 


4.3.3 Atmospheric Retrieval Method 


We use an adaptive Markov Chain Monte Carlo (MCMC, see Section 2.3.1 for an intro¬ 
duction) algorithm to retrieve the properties of an exoplanet’s atmosphere embedded in its 
transmission spectrum (Equation |3.3[ ). The key improvement compared to previous stud¬ 
ies (e.g., [Madhusudhan and Seager[ |2009[ [Benneke and Seag^ |2012[ [Eee et al.j |2012| ) is 
the use of a self-consistent set of parameters derived from first principle of transmission 

0.456 X 10^^ 


nco - 1 
1 X 10-8 


= 22851 -r 


(71427)2 - A 


-2 ■ 


(4.4) 


observe an additional water band at 3.2/rm in Kaltenegger and Traub (2009 1 ; it could be due to a 
complementary list used in this study. The HITRAN water absorption lines are uniformly spaced in terms of 
wavenumber, as expected from quantum mechanics. Furthermore, HITRAN’s aim is to provide the spectro¬ 
scopic parameters required to simulate the transmission and emission light for Earth-like conditions. There¬ 
fore, it is unlikely that such a significant water band would not be included in HITRAN. 
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spectroscopy (Section[3^ that uniquely constrains the planetary mass. In other words, we 
use as jump parameters: the temperature, the pressure scale height and the species num¬ 
ber density at a given planetary radius—whose choice does not affect the retrieval method, 
similarly to the “reference radius”. We assume a uniform prior distribution for all these 
jump parameters and draw at each step a random stellar radius based on a Gaussian prior 
assuming here a 1% relative uncertainty on the stellar radius. 


Synthetic exoplanet scenarios 


We assess MassSpec's capabilities for different super-Earths and Earth-sized planets. We 
use GJ 1214b’s ( Charbonneau et ahj 20091 and Earth’s sizes, respectively. We use the “hot 
mini-Neptune”, “Hot Halley world”, and “nitrogen-rich world” atmospheric scenarios in¬ 
troduced in |Benneke and Seager| ( |2012] ) to provide a representative overview of MassSpec’s 
potential. We assume a well-mixed atmospheres and use temperature profiles similar to 


Miller-Ricci et al. (2009)—in particular, T{p >1 mbar) ^ 300 K. 


We assume a total in-transit observation time of 200 hrs (Doming et al. 2009) and a 
MlVp^ host star located at 15 pc, except for the Earth-like planet observed with JWST 
for 200 hrs around a M7V at 15 pc. We note that for Earth-sized planets, we investigate 
in greater details the nitrogen-world scenario (“Earth-like” planet) because (1) the first 
exoplanets to be confirmed habitable are likely to be hot desert world ( Zsom et ahj 2013| ), 
(2) in most cases hydrogen should have escaped from an Earth-sized planet atmosphere, 
and (3) nitrogen-dominated atmospheres are less favorable for transmission spectroscopy 
due to their larger mean molecular mass (Equation |3.49| ) hence MassSpec's capabilities 
derived from those are more conservative. 


'a MIV star is a representative host in terms of transmission-spectrum significance (Section 


3.4.21. 
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4.3.4 Long-Term: Characterizing Planets as Small as Earth 

We find that future space based facilities designed for exoplanet atmosphere characteriza¬ 
tion will also be capable of mass measurements for super-Earths and Earth-sized planets 
with a relative uncertainty as low as ~ 2%—a precision that has not yet been reached using 
RV measurements, even for the most favorable cases of hot Jupiters. 

The pool of planets accessible to MassSpec will extend down to Earth-sized planets 
with the future observatories. We estimate that with data from JWST, MassSpec could 
yield the masses and atmospheric properties of mini-Neptunes, super-Earths, and Earth¬ 
sized planets up to distances of 500 pc, 100 pc, and 50 pc, respectively, for M9V stars and 
200 pc, 40 pc and 20 pc for Ml V stars or stars with earlier spectral types (Eigure [4^ . Eor 
an EChO-class mission, the numbers would be 250 pc, 50 pc and 13 pc, and 100 pc, 20 pc 
and 6 pc, respectively. 



Eigure 4-7: Boundaries of MassSpec'?, application domain for 200 hours of in-transit ob¬ 
servations. Using JWST, MassSpec could yield the mass and atmospheric properties of 
super-Earth and Earth-sized planets up to the distances shown by the black dashed, and 
dotted lines, respectively. Similarly, the maximum distances to a mini-Neptune, a super- 
Earth, and an Earth-sized planets using EChO's observations are represented by the blue 
solid, dashed, and dotted lines, respectively. The green dotted line refers to the case of an 
Earth-sized planet observed with a 20-meter space telescope. The grey area shows the stars 
too bright for JWSTINIRSpec in the R=1000 mode (J-band magnitude < 7). 
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Earth-sized Planets 


We find that MassSpec could determine the masses and atmospheric properties of Earth¬ 
sized planets using JWST or an EChO-class mission (Figure |4^ . For example, even for 
case of Earth-like planets, MassSpec could constrain their masses with a relative uncer¬ 
tainty of < 8% using JWST (Figure |4^ —if these are transiting late M dwarfs within 15 pc 
of Earth. In the future era of 20-meter space telescopes, sufficiently high quality transmis¬ 
sion spectra of Earth-sized planets will be available (Ehrenreich et al.[ 2006 [ Kaltenegger 


and Traub 20091. By using MassSpec, such facilities could yield the mass of Earth-sized 


planets transiting a MIV star (or stars with earlier spectral types) at 15 pc with a relative 
uncertainty of ~ 5% (Figure 4-9| ). For M9V stars, such an observatory enable the character¬ 
ization of Earth-sized planets up to 200 pc and for Ml V stars or stars with earlier spectral 
types, up to 80 pc (Figure[4^. 


Super-Earths 

We find that MassSpec would yield the mass measurements of super-Earths transiting an 
MIV star at 15 pc with a relative uncertainty of ~ 2%(~ 10%)[~ 15%] for 
hydrogen(water)[nitrogen]-dominated atmospheres (Figures [4-10[ |4-11| and |4-12| ), if ap¬ 
plied to 200 hrs of such a target. The larger significance of the mass measurements obtained 
for hydrogen-dominated super-Earths results from higher SNR of their transmission spec¬ 
tra (Equation |3.49[ ) —the small mean molecular masses of their atmospheres reduce their 
extent. For the same super-Earths with hydrogen(water)-dominated atmospheres, EChO's 
data would yield mass measurements with a relative uncertainty of ~ 3%(~ 25%) (Fig- 
;4-13 and 4-14|), respectively—for a nitrogen world in the same configuration it will 


ures 


not be possible to constrain its mass (Figure [4^T5] ). 


Giant Planets 

JWST's and EChO ’s observations of planets larger than super-Earths could provide high- 
significance mass and atmospheric constraints with a limited number of transits. Using the 
scaling law (Equation|3.49|), we estimate that for giant planets within 100 pc half a dozen 
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Figure 4-8: MassSpec's application to the synthetic transmission spectrum of an Earth-like 
exoplanet transiting a M7V star at 15 pc as observed with JWST for a total of 200 hrs 
in-transit. (A) Synthetic data and the best fit together with the individual contributions of 
the atmospheric species. (B) Normalized posterior probability distribution (PPD) of the 
atmospheric species number densities at the reference radius. (C) Normalized PPD for the 
scale height. (D) Normalized PPD for the pressure at deepest atmospheric level probed by 
transmission spectroscopy. (E) Normalized PPD for the temperature. (F) Normalized PPD 
for the exoplanet mass. The diamonds indicate the values of atmospheric parameters used 
to simulate the input spectrum and the asterisks in the panel A legend indicate molecules 
that are not used to simulate the input spectrum. The atmospheric properties (number 
densities, scale height, and temperature) are retrieved with significance yielding to a mass 
measurement with a relative uncertainty of ~ 8%. 
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Figure 4-9: MassSpec's application to the synthetic transmission spectrum of an Earth¬ 
sized planet with a nitrogen-dominated atmosphere transiting a MIV star at 15 pc as ob¬ 
served with a future-generation 20-meter space telescope for a total of 200 hrs in-transit. 
The panels show the same quantities as on Figure [4^ The atmospheric properties are re¬ 
trieved with high significance yielding to a mass measurement with a relative uncertainty 
of ~ 5%. Note that the significant observation of the Rayleigh-scattering slope combined 
with the lack of H 2 — H 2 CIA feature at 3 microns (Figure [4-10[ ) yields to the retrieval of 
nitrogen as the dominant atmospheric species. 
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Figure 4-10: MassSpec's application to the synthetic transmission spectrum of a hydrogen- 
dominated super-Earth transiting a MIV star at 15 pc as observed with JWST for a total of 
200 hrs in-transit. The atmospheric properties (number densities, scale height and temper¬ 
ature) are retrieved with high significance yielding to a mass measurement with a relative 
uncertainty of ~ 2%. Note the significant difference between the number density PPDs of 
hydrogen, water, carbon dioxide, and methane and those of ozone and nitrogen (B). The 
latter two gases were not part of the synthetic atmosphere. Ozone and nitrogen are not 
detected, because no constraints on their mixing ratios can be made. Hydrogen-dominated 
planets are targets that are particularly favorable for MassSpec because their extended at¬ 
mosphere leads to high-SNR transmission spectra. 


138 



























Normalized PPD Relative transit depth 



Figure 4-11: MassSpec’s application to the synthetic transmission spectrum of a water- 
dominated super-Earth transiting a MIV star at 15 pc as observed with JWST for a total of 
200 hrs in-transit. The panels show the same quantities as on Figure [4^ The atmospheric 
properties (number densities, scale height, and temperature) are retrieved with significance 
yielding to a mass measurement with a relative uncertainty of ~ 10%. 
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Figure 4-12: MassSpec’s application to the synthetic transmission spectrum of a super- 
Earth with a nitrogen-dominated atmosphere transiting a MIV star at 15 pc as observed 
with JWST for a total of 200 hrs in-transit. The panels show the same quantities as on 
Figure [4^ The atmospheric properties are retrieved with high significance yielding to a 
mass measurement with a relative uncertainty of ~ 15%. 
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Figure 4-13: MassSpec's application to the synthetic transmission spectrum of a super- 
Earth with a hydrogen-dominated atmosphere transiting a MIV star at 15 pc as observed 
with EChO for a total of 200 hrs in-transit. The panels show the same quantities as on 
Figure [4^ The atmospheric properties are retrieved with high significance yielding to a 
mass measurement with a relative uncertainty of ~ 3%. Note that EChO's capabilities for 
mass measurements are enhanced by its large spectral coverage that yields the Rayleigh- 
scattering slope, which is particularly valuable for mass and atmospheric retrieval. 
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Figure 4-14: MassSpec's, application to the synthetic transmission spectrum of a super- 
Earth with a water-dominated atmosphere transiting a MIV star at 15 pc as observed with 
EChO for a total of 200 hrs in-transit. The panels show the same quantities as on Fig¬ 
ure 4-8 The atmospheric properties are retrieved with high significance yielding to a mass 
measurement with a relative uncertainty of ~ 25%. 
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Figure 4-15: MassSpec’s application to the synthetic transmission spectrum of a super- 
Earth with a nitrogen-dominated atmosphere transiting a MIV star at 15 pc as observed 
with EChO for a total of 200 hrs in-transit. The panels show the same quantities as on 
Figure [4^ The planet mass is not retrieved because the data quality does not yield the 
atmosphere composition and temperature—although the scale height and the signatures of 
water, methane, and carbon dioxide are retrieved. 
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of transits would secure mass measurements with a relative uncertainty < 5%. An EChO- 
class mission would be particularly useful to characterize giant planets as its wide spec¬ 
tral coverage would allow to measure their Rayleigh-scattering slope at short wavelengths. 
MassSpeEs applications to giant planets will be particularly important for giant planets 
whose star’s activity prevents a mass measurement with RV (e.g., the hottest known planet, 


WASP-33b |Collier Cameron et al.||2010| ). In addition, independent high-significance mass 
measurements of targets accessible by RV will also be beneficial to compare the capabilities 
of both techniques and to assess their future complementarity. 


4.3.5 Discussion 


How Many (Habitable) Planets Are Within Reach? 


We present in Table |4.2| estimates of the number of exoplanets that would be characteri- 
zable with MassSpec, or any other atmospheric retrieval methods based on transmission 
spectroscopy. These estimates assume a stellar density of 0.14 stars per pc^ and consider 
that 80% of those are M dwarfs. We derive the planet’s transit probability (~ R^/a) as¬ 
suming them to be at the center of their host’s habitable zone (a = (i?*/i? 0 )^(T*/T 0 )'^). 
We use a planet occurrence rate of 0.15 for Earth-sized planets in their host’s habitable 
zone and 1 for other planet types ( Dressing and Charbonneau] 20131. The number within 
brackets refer to the estimate assuming the most optimistic occurrence rate of HZ planets 
from Zsom, in prep. 

We note that about 75% of the planets to be accessible for in-depth characterization orbit 
a host with a spectral type later than M5V. We develop this point in Figure [4%^ and show 
the estimated fraction of characterizable planets as a function of the host star’s type. As it 


was pointed out in Section 3.4.2 late M dwarfs are favorable for transmission spectroscopy 
because of their large ratio of radiance over projected area (Figure |3^ . For this reason, 
atmospheric retrieval methods are applicable to larger distances for late M dwarfs than 
for stars with earlier spectral types (Figure |4^ . Hence, the number of favorable hosts 
is considerably larger as it scales with the cube of the distance. In addition, the transit 
probability is larger for stars with later spectral types. 
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Table 4.2: Estimated number of exoplanets within reaeh for atmosphere eharaeterization. 
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Figure 4-16: Estimated fraction of exoplanets characterizable via transmission spec¬ 
troscopy as a function of the host star’s type. 


This new insight for the future exoplanet characterization highlight that: 

(i) If Earth-sized planets are found around the habitable zone of M dwarfs before 2018, 
one could constrain their atmospheric properties and mass with future facilities such as 
JWST. In other words, the first habitable distant worlds could be identified within the 
next decade. 


(ii) If found by 2018, there will be too many favorable planets to be characterized with 
the available facilities. As an example, the sole characterization of each of the ~ 200 
possibly-habitable Earth-sized planets would require a continuous use of JWST over its 
lifetime (5 years). Not only will JWST have other scientific goals for the field of exo¬ 
planetology, but it is also not dedicated to planetary science. Hence, it will be key for 



M9V M8V M7V M6V M5V M4V M3V M2V M1V MOV 

Host star’s type 
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the community to pre-select targets of interest—which requires to discover these well in 
advance of JWST's launch. 


(i) Detecting the Earth-Sized Planets to he Identify Hahitahle in the Next Decade 

JWST’s optical performance could yield the characterization of about 200 Earth-sized plan¬ 
ets in the habitable zone of their host stars, primarily late M dwarfs. All of them are yet to 
be detected. Their early detection (prior to 2018) is key to ensure that JWST would moni¬ 
tor enough of their transits to secure a sufficient SNR. In practice, half a dozen Earth-sized 
planets in the habitable zone of their host stars could be properly characterized over JWSTs 
lifetime. To date, one mission is dedicated to this aim: SPECULOOS {Search for Habit¬ 
able Planets Eclipsing Ultra-cool Stars). SPECULOOS is a European Research Council 
mission that will begin observing the coolest M dwarfs in 2016. The mass of the com¬ 
panions found will not be constrained by RV because of the faintness of their host stars. 
However, MassSpec's application to their JWST spectra would yield both their masses and 
atmospheric properties (Eigure [4^ . Hence, MassSpec is likely to play a pivotal role in the 
assessment of their potential habitability. 


(ii) Optimizing the Scientific Return of JWST 

Time prioritization will be key for the future of exoplanetology. The number of targets 
favorable for characterization and the time it would require will be too large to allow for a 
systematic characterization due to our limited capabilities. That is why we suggested to aim 
for community efforts (1) to discover targets of interest well in advance of JWSTs launch 
and (2) pre-select the targets for which JWST's time will be dedicated over its life-time. To 
help mitigate the drawback of dedicating JWSTs time to a limited amount of targets, we 
suggest to consider more modest missions such as EChO. The smaller aperture of EChO 
would enable it to observe brighter stars (i.e., early-type and close-by stars), hence EChO'?, 
and JWST'? time could be respectively prioritized on super-Earths and Earth-sized planets 
for M9V stars closer than 25 pc and for MIV stars (or stars with earlier spectral type) 
closer than 10 pc (Eigure [4^ ). Similarly, EChO'? and JWST'? time could be respectively 
prioritized on giant planets and super-Earths for M9V stars closer than 125 pc and for Ml V 


146 








stars (or stars with earlier speetral type) eloser than 50 pe. EChO would be partieularly 
useful to determine the mass—and atmospherie properties—of giant planets beeause its 
wide speetral eoverage would allow to measure their Rayleigh-seattering slope at short 
wavelengths. 

The same argument eoneerning the synergy between JWST and an EChO-class mission 


in terms of mapping eapabilities ean be found in Seetion 4.2.4 


Clouds Will Not Overshadow MassSpec 


Clouds are known to be present in exoplanet atmospheres (Demory et al. 20131 and to af- 


feet transmission speetra beeause atmospherie layers below the eloud deek are not probed 
( [Benneke and Seag^ |2013| [Barstow et al.| , |2013b| ). However, MassSpec is not rendered 
ineffeetual by elouds beeause elouds are deeper than the lowest pressure probed by trans¬ 
mission speetroseopy. In other words, there will always be atmospherie information avail¬ 
able from transmission speetroseopy, in theory. In praetiee, speetral features would be 
deteetable in the presenee of elouds but with a redueed signifieanee (eompare panels A of 
Figures [4-17| to [4H9| ): the higher the eloud deek, the larger is the uneertainty on the mass 
estimate due to the redueed amount of atmospherie information available. We show the 
effeet of eloud deeks at 100, 10, and 1 mbar in Figures [4-17[|4-18l and 4-19[ respeetively. 
Note that 100 mbar eorresponds to Earth’s and Venus’ elouds’s pressure level and 1 mbar 


is the lowest pressure where thiek elouds are expeeted (Howe and Burrows 2012). Cloud 
deeks at 100 and 10 mbar affeet marginally MassSpec'?, eapabilities beeause a limited frae- 
tion of the speetral bins probe as deep as than 10 mbar in transmission. However a 1-mbar 
eloud deek affeets signifieantly atmospherie retrieval results, inereasing the uneertainty on 
the mass estimate by a faetor of ~ 2 over the one derived in the eloud free seenario. The 
inereased uneertainty in the mass estimate results mainly from an inereased uneertainty on 
the seale height and the mean moleeular mass estimates. The reason is that (1) the preei- 
sion on seale-height estimates are redueed beeause the signatures of atmospherie speeies 
are truneated and (2) additional atmospherie seenarios are now possible beeause numerous 
eombinations of moleeular signatures eould be masked by a high-eloud deek. Disentan¬ 
gling between sueh seenarios requires a suffieient speetral resolution. 
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Figure 4-17: MassSpec’s application to the synthetic transmission spectrum of a super- 
Earth with a water-dominated atmosphere presenting a cloud deck at 100 mbar (such as 
Venus) transiting a Ml V star at 15 pc as observed with JWST for a total of 200 hrs in-transit. 
The panels show the same quantities as on Figure [4^ The atmospheric properties are 
retrieved with high significance yielding to a mass measurement with a relative uncertainty 
of ~ 10%. MassSpec's capabilities are not affected by the 100-mbar cloud deck because 
transmission spectroscopy does not probe deeper than 10 mbar (Figure[4^ 11 [ panel D). 
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Figure 4-18: MassSpec’s application to the synthetic transmission spectrum of a super- 
Earth with a water-dominated atmosphere presenting a cloud deck at 10 mbar transiting 
a MIV star at 15 pc as observed with JWST for a total of 200 hrs in-transit. The panels 
show the same quantities as on Figure [4^ The atmospheric properties are retrieved with 
high significance yielding to a mass measurement with a relative uncertainty of ~ 10%. 
MassSpec's capabilities are marginally affected by a 10-mbar cloud deck because a limited 
fraction of the spectral bins probe deeper than 10 mbar. 
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Figure 4-19: MassSpec’s application to the synthetic transmission spectrum of a super- 
Earth with a water-dominated atmosphere presenting a cloud deck at 1 mbar (lowest pres¬ 
sure level where clouds are expected Howe and Burrows 20121 transiting a MIV star at 
15 pc as observed with JWST for a total of 200 hrs in-transit. The panels show the same 
quantities as on Figure [4^ Despite the presence of high and thick clouds, the atmospheric 
properties are retrieved with sufficient significance to yield a mass measurement with a 
relative uncertainty of ~ 20% (i.e., twice larger than in the cloud-free scenario). 
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Number Densities vs Mixing Ratios in Atmospheric Retrieval Method 

Species’ mixing ratios are not adequate parameters for atmospheric retrieval method. The 
reason is that the uncertainty on mixing ratios encompasses the uncertainty on the num¬ 
ber densities and the pressure. We demonstrate that pressure and number densities are the 
atmospheric parameters embedded in a planet’s transmission spectrum. This implies that 
pressure and mixing ratios are correlated in the context of transmission spectroscopy (Fig- 
ure |4-20[ A). Hence, the posterior probability distributions of mixing ratios are the results of 
the convolution of the number densities’ and pressure’s PPDs (Figure [4-20[ B). Therefore, 
the significance of a molecular detection based on the mixing ratio will be less significant 
than if based on the number densities, in particular because number densities are more 
constrained than the pressure. 




Figure 4-20: The use of species’ mixing ratios in atmospheric retrieval method. (A) 
Marginal posterior probability distribution (69% and 95% confidence intervals, respectively 
in red and blue) of maximum pressure probed and H 2 O mixing ratio shows the correlation 
between both parameters. The correlation shown between the pressure and a species mix¬ 
ing ratio translates the fact that the key parameter of a planet’s transmission spectrum are 
the number densities, not the mixing ratios. Therefore, the uncertainty on mixing ratios 
combines the uncertainty on the pressure and the number densities. (B) Water’s mixing 
ratio PPD (dot-dash line) and water’s number density PPD (solid line), shifted along the 
x-axis for comparison. The larger extent of the water mixing ratio PPD results from the 
combination of the uncertainty on the number density and pressure. 
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No False Molecular Detections 


MassSpec’s applications do not lead to false moleeular deteetions—assuming that all lines 
and bands of a moleeule are known. As emphasized previously, the applieation of MassSpec 
requires a planet’s transmission speetrum of suffleient quality (SNR, speetral resolution, 
and eoverage) to identify the major moleeular speeies. Onee these eriteria are met, the 
speeifieity of the speeies’ signature in extinetion (shown in panels A of Figures |4-8| to [AT9 1 
prevents from the false positive deteetion of atmospherie speeies. As an example, we try to 
retrieve ozone that was not part of the synthetie atmospheres and it is deteeted in none of 
our retrieval. 
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Chapter 5 


Summary and Conclusions 


The detailed characterization of extra-solar planets is the next milestone for the field of 
exoplanetology. The quality of the data available to characterize exoplanets is constantly 
improving. In that context, we approached the question “What new insights into the at¬ 
mospheric and interior properties of exoplanets will future high-quality data bring within 
reach?”. For that purpose, this thesis focused on transiting exoplanets, which are of special 
interest due to the wealth and diversity of data made available by their favorable orbital 
configuration. We introduced two new characterization methods. The first method aims to 
map transiting exoplanets’s atmospheres based on their light curves, with a particular focus 
on their secondary eclipses. The second method aims to determine transiting planet masses 
and atmospheric properties solely from transmission spectra, i.e. the starlight filtered by a 
planet’s atmosphere during transits. 

For both methods, we introduced the underlying concepts and theory, emphasized the 
possible caveats, and presented the Bayesian frameworks developed to analyze available 
datasets and assess their potential with future observatories. Below we summarize our main 
results and conclusions for each method. In particular, we provide a takeaway message for 
each method, a reminder of what is possible at the time this thesis was written, and what the 
future has to offer. We conclude with general recommendations for the field of exoplanet 
characterization. 
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5.1 Mapping Exoplanet Atmospheres 


Takeaway message: We found that a planet’s shape, brightness distribution, and system 
parameters are underdetermined by its light curve (i.e., correlated). This means that as¬ 
sumptions on any of these properties affect the others and yield biased estimates and also 
underestimated error bars. For example, we show that assuming HD 189733b to be uni¬ 
formly bright in the Spitzerf\KhC 8 /rm channel leads to an overestimation of the stellar 
density (and hence, the planetary density) by 5%. For that reason, we developed a method 
to map exoplanets based on global analyses using datasets complimentary to transits, such 
as radial velocity measurements, to relax assumptions. 

Present: We mapped the hot Jupiter HD 189733b’s emission in the infrared (8 /im) in 
2D and constrained the characteristic of its hot spot. We mapped Kepler-7b in the visible 
and contributed to the first map of cloud coverage on a planet outside of the Solar System. 

Future: Exoplanet mapping will provide unprecedented insights into the physics of 
exoplanetary atmospheres via their temperature and composition patterns. By 2015, we 
will obtain the first 3D maps of exoplanet atmospheres. By 2020, we should be able to 
assess the time variability of 3D structures in exoplanet atmospheres. 


Summary of key points: 

• Mapping exoplanets is an underdetermined problem because the target’s brightness 
distribution is underdetermined by its light curve, hence multiple solution can fit 
a given dataset. Therefore, mapping methods have to ensure that all the possible 
solutions are revealed. 

• Mapping exoplanets also involves degenerate solutions because four of the systems 
parameters (namely, e, uj, b, and p*) and the target’s projected shape at conjunc¬ 
tions can also affect the occultation shape, hence those are correlated with its bright¬ 
ness distribution. To mitigate the effects of this degeneracy, it is crucial to avoid 
unassessed assumptions on one of the above properties as it could bias the estimates 
of the others. For example, we show that assuming HD 189733b to be uniformly 
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bright in 5piYzer/IRAC-8/im channel leads to an overestimation of the stellar density 
(and hence, the planetary density) by 5%—which is at 6 a of the actual estimate. We 
advocate for relaxing assumptions in global analyses of occultation, transits, phase 
curves, and RV measurements to consistently probe the correlated parameter space of 
the light-curve fitting problem. We note that the effects of assumptions—such as “the 
planet can be modeled as a uniformly-bright disk”—are adequate to fit light curves 
with low SNR, typically when the ratio eclipse of depth to photometric precision on 
a ~l-min bin is 10. 

• We applied our mapping method to the hot Jupiter HD 189733b. We detected the 
deviation of its occultation shape from the one of a uniformly bright disk at the 6 a 
level in Sp/tzer/IRAC’s 8/rm channel. We demonstrated that this deviation emerges 
mainly from a large-scale brightness structure in HD 18733b’s atmosphere—we re¬ 
jected HD 189733b’s shape as a possible contributing factor based on its transit resid¬ 
uals. HD 18733b’s brightness distribution indicated a hot spot shifted east of the 
substellar point within the atmospheric layers probed at 8 /rm. In addition, we in¬ 
vestigated the practical effects of underlying model assumptions on our inferences 
due to the correlation between the system parameters and the brightness distribu¬ 
tion. Notably, we found that the more complex HD 189733b’s brightness model, the 
larger the eccentricity, the lower the densities, the larger the impact parameter and 
the more localized and latitudinally-shifted the hot spot estimated. We outlined the 
influence of HD 189733’s RV measurements on our inferences obtained solely from 
SpzYzer/IRAC 8-/rm photometry. In particular, we observed redistributions of the 
probability density for the posterior probability distributions (PPDs) obtained with 
complex brightness models, which result from the rejection of solutions involving 
^/esinoj > 0.15 by HD 189733’s RV data—favored by the photometry for complex 
brightness models. This reemphasized the necessity of global analyses to consis¬ 
tently probe the correlated parameter space of exoplanets’ light curves and results in 
a new upper limit of HD 189733b’s orbital eccentricity, e < 0.011 (95% confidence). 

• We applied our mapping method to the hot-Jupiter Kepler-7b as a contribution to 
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Demory et al. (20131. We found that the large geometrie albedo of Kepler-7b in 


Kepler's bandpass originates mainly from the western seetor of its dayside. In partie- 
ular, we found Kepler-7b has a bright seetor extending from the western terminator 
of its dayside to ~ 30° east of the substellar point with a mean relative brightness of 
78 ± 4 X 10^ ppm. 


• We highlighted that multi-wavelength observations enable 3D mapping of exoplanet 
atmospheres and mitigate the degeneraey between the brightness distributions and 
the system parameters. Multi-wavelength observations probe different optieal depths 
under the same orbital eonfiguration, henee yielding 3D mapping while reeording the 
same information about the system parameters. 


• We showed how applying our method to eomplementary observations of HD 189733b 
(and HD 209458b) eould yield to the first 3D map of exoplanet atmospheres within 
a year (as a eontribution to the Spitzer proposal #10103, PI Dr. Nikole Lewis, see 
Lewis et aH|2013a| ). 

• We determined the best known exoplanets to perform eelipse mapping in the visible 
using HST/WFC3 are WASP-18b, WASP-19b, and WASP-76b (TablejAT]). 

• We showed that future spaeed-based observatories like JWST and EChO-class mis¬ 
sions eould yield time-dependent 3D maps of distant worlds. 


• We pointed out the synergy of JWST and EChO-class missions to eonstrain the eloud 
eoverage of an exoplanet. 


Finally, we diseussed how the implementation of atmospherie eireulation models 
(e.g., Heng and Workmanl 2014[ ) in our mapping method eould yield direet eon- 
straints on the atmosphere’s eireulation parameters while mapping it. 
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5.2 Insights into Exoplanets from Transmission Spectroscopy 


Takeaway message: We demonstrated from first prineiples that a planet’s mass is a key 
parameter behind its transmission speetrum, akin to the planet’s atmospherie properties 
(temperature, pressure, and eomposition). Henee, for eonsisteney, atmospherie retrieval 
methods have to aeeount for the planetary mass. On this basis, we developed a new retrieval 
method ealled MassSpec. 

Present: We applied a simplified version of MassSpec to the hot-Jupiter HD 189733b 
and retrieved a mass in exeellent agreement with the mass estimated from RV measure¬ 
ments. We showed how the same proeedure eould determine the mass of the hottest planet 
known, WASP-33b, with less than 20 HST orbits. 

Future: We demonstrated MassSpec'?, eapability to determine the masses and atmo¬ 
spherie properties of planets as small as Earth with future observatories. We estimated that 
JWST eould eharaeterize the properties of ~ 200 Earth-sized planets in their hosts’ habit¬ 
able zones, if sueh planets are diseovered. The vast majority of these hosts will be late M 
dwarfs, whieh may prevent a mass determination via RV measurements. Eor that reason, 
it is likely that MassSpec will play a pivotal role in the identifieation of the first habitable 
exoplanet, whieh eould oeeur within the next deeade. 


Summary of key points: 

• To date, mass eonstraints for exoplanets are predominantly based on radial veloeity 
measurements, whieh are not suited for planets with low masses, large semi-major 
axes, or those orbiting faint or aetive stars. We presented an alternative method to the 
RV method, MassSpec, to determine the mass of an exoplanet solely from transmis¬ 
sion speetroseopy. MassSpec also provided a solution to the problem of redueed at¬ 
mospherie retrieval quality based on a planet’s transmission speetrum when its mass 
is inadequately eonstrained (see Barstow et al.[ 2013a| ) by eonsistently determining 
(i.e., simultaneously) a planet’s mass and its atmospherie properties. 


• We prove MassSpec'? feasibility analytieally by identifying the key parameters of 
and their effeets on an exoplanet’s transmission speetrum. In partieular, we demon- 
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strated that a transmission spectrum depends in unique ways on the atmospheric 
scale height, reference pressure, temperature and number densities of the main at¬ 
mospheric absorbents—which constrain the mean molecular mass. The uniqueness 
of these dependencies enables the independent retrieval of each of these key param¬ 
eters that in turns leads to the planet’s mass. 


We found that the Euler-Mascheroni constant ('fEM, see |Euler[|1740[ ) plays a pivotal 
role in the equations of transmission spectroscopy. In particular, 7EM is central in 
the identification of the key parameters of transmission spectrum. Eurthermore, we 
showed that it is also key to improving the computational efficiency for modeling 
transmission spectra via solving hf,ff{X) = [z : t{z, A) = 


• We tested MassSpec on real data using HD 189733b’s transmission spectrum and 
showed an excellent agreement between the masses retrieved by MassSpec and from 
RV mass measurements. 


• We discussed the complementary nature of RV measurements and transmission spec¬ 
tra to constrain planetary masses. To do so, we expressed the sensitivity of each 
signal’s significance to the system parameters. 

• We reemphasized the benefit of late M dwarfs for transit observations. We showed 

that the significance of in-transit signals scales as and is thus increased 

up to a factor of three when the host-star of a given planet is a late M dwarf, instead 
of, e.g., a Sun-like star. In that context, we also discussed the characterization of 
habitable earth-sized planets around late M dwarfs in the next decade. 

• We demonstrated MassSpec'?, capacity to yield the mass of the hottest planet known, 
WASP-33b, with a relative uncertainty of < 20% with less than 20 HST orbits. 

• We demonstrated numerically that future space-based facilities designed for exo¬ 
planet atmosphere characterization will also be capable of mass measurements for 
super-Earths and Earth-sized planets with a relative uncertainty as low as ~ 2%—a 
precision that has not yet been reached using RV measurements, even for the most 
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favorable cases of hot Jupiters. In particular, we showed that JWSTs observations of 
unfavorable cases, such as Earth-sized planets with nitrogen-dominated atmosphere 
(“Earth-like planets”), could yield mass measurements with a relative uncertainty of 
~ 8%—if the targets are transiting late M dwarfs within 15 pc of Earth. We esti¬ 
mated that with data from JWST, MassSpec could yield the mass of mini-Neptunes, 
super-Earths, and Earth-sized planets up to distances of 500 pc, 100 pc, and 50 pc, 
respectively, for M9V stars and 200 pc, 40 pc and 20 pc for MIV stars or stars with 
earlier spectral types. Eor EChO-cXd&s missions, the numbers would be 250 pc, 50 
pc and 13 pc for M9V stars and 100 pc, 20 pc and 6 pc for MIV stars or stars with 
earlier spectral types. 

• We estimated JWST could determine the mass and the atmospheric properties of 
~2000 Earth-sized planets, 10% of which would be in their host’s habitable zones. 
Eor an EChO-class mission, this numbers would be 40. This highlights how close we 
are to identifying habitable exoplanets. 

• We discussed MassSpeEs input to gain insight into planetary interiors. Although 
MassSpec provides constraints on the mass and radius of a target, those are not al¬ 
ways sufficient to obtain insights into a planet’s interior. MassSpec's simultaneous 
constraints on a planet’s atmosphere and bulk density may help to break this degen¬ 
eracy, in some cases. 

• We emphasized that joint analyses of transmission spectra with any dataset that can 
independently constrain a planet’s mass and/or atmospheric properties can enhance 
MassSpec’s capabilities. In particular, RV measurements and emission spectra would 
be great assets to provide complementary constraints on the mass and the atmo¬ 
spheric properties, respectively. 

• We showed that, in principle, MassSpec is not rendered ineffectual by clouds because 
clouds are deeper than the lowest pressure probed by transmission spectroscopy. In 
other words, there will always be atmospheric information available from transmis¬ 
sion spectroscopy, although high-altitude clouds may increase the uncertainty on the 
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planetary mass estimates by a faetor of ~ 2. 


5.3 General Recommendations for Future Progress 

5.3.1 On the Need for Accurate Extinction Cross Section Databases 

Like most methods aimed at determining the properties of distant objeet, the methods in- 
trodueed in this thesis require aeeurate extinetion eross seetions (Equation |3.4[ ). On the 
one hand, 3D mapping of exoplanetary atmospheres neeessitates the determination the at- 
mospherie layer from whieh the deteeted flux originates. On the other hand, atmospherie 
retrieval methods—like MassSpec —solve the inverse problem of determining the eondi- 
tions (i.e., temperature, pressure, eomposition) of the medium probed from its transmis¬ 
sion speetrum. In both oases, the key to a data interpretation is to know a priori how 
the transmission of light is affeoted by a medium’s oonditions (i.e., temperature, pressure, 
eomposition). That is why, we advoeate for devoting signifioant efforts to generate aoeu- 
rate extinetion eross seotion databases that oover various atmospherie eonditions in order 
to prevent them from limiting our ability to eharaeterize exoplanets. 

5.3.2 How to Optimize the Scientific Return of JWST 

The number of targets favorable for eharaeterization and the time it would require will 
exeeed the time allotted to exoplanets with JWST. This is why we suggested to aim for 
eommunity efforts (1) to diseover targets of interest well in advanee of JWSTs launeh with 
missions like SPECULOOS and (2) pre-seleet the targets for whieh JWST's time will be 
dedieated over its lifetime. To help mitigate the drawbaeks of dedieating JWST's time to 
a limited number of targets, we suggest to eonsider more modest missions sueh as EChO. 
We pointed out the synergy between JWST and an EChO-class mission to map exoplanets’ 
atmospheres, eonstrain their eloud eoverage, and determine their atmospherie properties 
and masses. 
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Appendix A 


Introduction to Rational Functions 


A rational function, /(.), is a function that can be expressed as the ratio of polynomials, 
P{.)/Q{.). The zeros and poles of /(.) are the zeros of P(.) and the zeros of (5(.)> respec¬ 
tively. Rational functions are used in various fields of science and engineering (signal pro¬ 
cessing, acoustics, aerodynamics, structural dynamics, electronic circuitry, control theory, 
etc.) as mathematical representation of the relation of the inputs (i(.)) and outputs (o(.)) 
of a system, called the transfer function of the system. The transfer function of a system 
{H (s)) is the linear mapping of the Laplace/Fourier transform of its inputs {I{s) = C{i{.))) 
to the Laplace/Fourier transform of its outputs (0(s) = £(o(.))), i.e., 0(s) = H{s)I{s) 
where s is a spatial or temporal frequency if i and o are functions of space or time, respec¬ 
tively. As an example, the transfer function of imaging devices is the Fourier transform 
of the point spread function (PSF)—because the theoretical input leading to the PSF is a 
spatial impulse (i.e., a Dirac function) in the field of view and the Fourier transform of an 
impulse is one, therefore C{PSF) = H{s). 

The transfer function of a system relates directly to the differential equations used to 
represent mathematically the system. Therefore, the poles and zeros of a transfer function 
characterize the behavior of a system. For example, the differential equation of a second 
order system is 

X + 2(UnX + u'^x = 0, (A.l) 

where ( and Un are the damping ratio and the natural frequency of the system, respectively— 
e.g., for a mass-spring-dashpot system C, = cj(l\/km) and Un = \/k/m where m, c, and 
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Figure A-1: Typical behavior of rational functions. Left: Transfer function of a mass- 
spring-dashpot system with a natural frequency of 10 Hz and a damping ratio of 1.5. The 
transfer function is independent of the frequency at frequency lower than the system’s nat¬ 
ural frequency. At frequency larger than the natural frequency, the two poles of the transfer 
function are activated and the dependency of the transfer function on the frequency changes 
from oc to oc a; ^ (i.e., the exponent decreases by two). Right: Transfer function with 
a zero at 0.1 Hz and a pair of conjugated poles at ±10 Hz. The effect of the zero is to 
increase the exponent of the transfer function dependency on the frequency by one, while 
the effect of the pair of poles is to decrease it by two. 


k are the mass, the spring constant, and the damping coefficient, respectively. The transfer 
function of a second order system is 


H{s) 


K 






(A.2) 


where K is the system gain (K = 1 for the mass-spring-dashpot system) and s = juj where 
j is the imaginary unit (see Figure [AH] for more details). Therefore, at frequencies low 
compared to a;„, H{ju) ^ K, meaning that the response of a mass-spring-dashpot system 
to a low-frequency input, i{t) = A exp(ju;t), is the input, i.e. o(t) = i(t). At high frequen¬ 
cies, H{juj) K, K/{juj)'^ = —K/u'^ meaning that the response is in opposition of phase 
with the input and its amplitude is proportional to 

The position of a transfer function poles and zeros in the complex field (C) define the 
domains of different dependency regimes for the transfer function (Figure[AH]). In the 
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neighborhood of a zero, the exponent of the transfer function dependency to its variable 
increases by 1, while it decreases by 1 in the neighborhood of a pole. Hence, the overall 
shape of transfer function relates directly to its formulation. Therefore, from simulations 
or measurements, it is possible to derive the appropriate formulation to represent mathe¬ 
matically a system using a rational function. 
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Appendix B 


Towards a Generalization of 


Transmission Spectrnm Eqnations 


In this Appendix, we aim to extend the analytieal derivations introdueed in Seetion[3^ 
In partieular, we address here the ease of a non-isothermal atmosphere. The ease of a 
non-isoeompositional atmosphere is derived analogously. Thus, we do not provide the 
derivation here. 

B.l General Formulations for Atmospheric Quantities 

B.l.l Temperature 

The temperature profile of a planet’s atmosphere ean be formulated as 


n 



(B.l) 


where To is the temperature at and the (3h are the eoeffieient of the Maelaurin series 
expansion of the temperature profile. Note that (3^ = 1. For the rest of the demonstration, 
it is suitable to write T{z) as follow 


T{z) = ToMz), 


(B.l) 


165 




because it allows to rewrite T{z) as a function of the zeros of friz), oj : 

T(z) = To n (B.3) 

{Vuj:fT{(zj)=0} 

Note that the smallest positive zero of friz) define the domain in which the series ex¬ 
pansion is physically adequate—i.e., from 0 to ojq, ({ccq : uq < u,\/u > 0}). In addition, 
friz) cannot physically have zeros with a multiplicity greater than 1. 


B.1.2 Pressure 


The pressure profile of a planetary atmosphere can be derived from the equation of hydro¬ 
static equilibrium, 


dln(p) 


~dz, 

1 


Hofriz 


-dz, 


(BA) 

(B.5) 


where Hq is the atmospheric scale height at Rp^. Because friz) has no zero with a mul¬ 
tiplicity greater than 1, f^^iz) can be expressed as a sum of irreducible rational functions 
(see Appendix[A|): 




{Va;:/T(t<j)=0} 


R. 


{z-u)' 


(B.6) 


where R(^ are the residues of frplO in the context of rational function expansion in a 
irreducible rational function. In particular. 


R 


UJ 


1 

/r(^) ’ 


(B.7) 


where /^(.) is the derivative of /^ !(.), because /r(-) has no zero with a multiplicity greater 
than 1. 
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The second term of Equation|B.4|can then be integrated by parts: 


1 




{'iu]:fT{ui)=0} 




dz = 




{'iuj:fT{uj)=0} 


- log(z - Uj) 


+ Co, (B.8) 


where Cq is the integration constant. As a result, one can express the pressure profile as 


p{z) = Cl {z — u) (B.9) 

{Vui:fTiu]}=0} 

= C 2 JJ 

{Vu]:fT{uj)=0} 

where Ci and C 2 are integration constants. The integration constants are determined using 
the identity p{z = 0) = Pq, where Pq is the pressure at Rp^. For example, C 2 = Pq. 
Therefore, the generalized formulation of the pressure profile is 

= po n • (B.ii) 

{Va;:/T((^)=0} 


B.1.3 Number Densities 


The number-density profiles can be derived from the ideal gas law 

pV = nRT, (B.12) 

where V is the volume of gas, n is the gas number density, and R is the ideal gas constant 
(~ 8.314J.K“^mol“^). In particular, the generalized number-density profile can be derived 
from: 


n{z) = no — — , 
T Po 


(B.13) 
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where no is the atmospheric number density at Rp^. Therefore, using Equations |B.3| and [B.ll 
one obtains 


{'iuj:fT{oj)=0} 




(B.14) 


B.1.4 Slant-Optical Depth 

The slant-optical depth profile (Equation|3.4|) can now be written as 


t{z,X) = 2noao{X)T^pQ 


n 


X 




Fiuj) 


dx, (B.15) 


where Xoo = is the boundary of the series-expansion domain in the x- 

space and 


assuming the cross section to be of the forn^a{X) a Eor example, for Rayleigh 

scattering and CIA, {<5 = 0, e = 0} and {5 ~ —l,e = 1}, respectively. 

We find no solution to the integral in Equation |B .IS] using standard mathematical func¬ 
tions. However, we find analytical formulations for the three following simplified cases: 

1. the case where friz) has a single zero, uq {uq E M’*'), 

2. the case where friz) has two zeros, wq, (coq E M’*' and ui E M“), 


3. the case where friz) has three zeros, uq, uji,U 2 (cuq E M"*", uji,U 2 G M ). 


'For molecular absorption, the absorption coefficient (Equation 3.331 takes the form 


a{z,X) = ar^i^r. 


771=0 


n 

{V(.j:/t(i.j)=0} 


O-S) 


771—1 — 


Pjz) + I 


(B.17) 


p'^{z) + bi^' 

The form of a(A, T,p) provides an analytical form for t{z, A) for each of the 4 regimes introduced in Sec¬ 


tion 


3.2.3 via an integration by parts analogous to the derivation in Section B.1.4 However, the transmission 


spectrum equation (Equation |3.5| ) cannot be expressed in terms of standard mathematical functions because 
these analytical forms for t{z, A) have multiple zeros (see SectionB.2 1 . Therefore, we do not address here 
the case of molecular absorption. 
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The general formulation for the slant-optieal depth profile is 
r(z,A) = r(0,A) J] ~ 

{Va;:/T(a;)=0} 

|B.1.4[ 1 T{z) is a first-order polynomial 

In this ease, the temperature deereases linearly with the altitude, i.e.. 


Mz) = 1--. 

Uo 


(B.19) 


Henee, 


friz) = - 


1 

Wo 


(B.20) 


and. 


t{z,X) = 2n{z)aoiX)TQPQ 


It ean be shown that 


1 - 


X 


2(-Rp,o + z)uq 




dx. (B.21) 


-'I 


dx = 




(B.22) 


and. 


2^1 (a, 6 ; c; 1 ) 
2 Fi(a, 6 ; c; 0 ) 


r(c)r(c-a-6) 

r(c-a)r(c-b) 


(B.23) 


where 2-^1 (a, c; z) is the Gaussian hypergeometrie funetion and r( 2 ;) is the Euler integral 
of the seeond kind: the Gamma funetion. Therefore, Equation|B.21|ean be written as 


t{z,X) = + z)oJoao{X)T^pln{z) ° 

= ^27r(i?p,o + ^)woao(A)To^p^no-^^ + 1 _ _ 1 )) 


r (5 + i-^(e-i))’ 


(B.24) 


1 - - ) (B.25) 

^ 0 , 
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|B.1.4[2 T{z) is a second-order polynomial 


In this case, the temperature profile is a seeond-order polynomial with two zeros, one pos¬ 
itive and the other negative—to be physieally plausible. 


Henee, 


and. 





UJl-UQ 

UQLJl 

^ 0-^1 

UJQLJl 


(B.26) 


(B.27) 


t{z,X) = 2n{z)ao{X)T^pQ 





_ gX _ 

2{Rpfi+z)uji 


Fioji) 

dx, (B.28) 


where 


n{z) 


no 




(B.29) 


It ean be shown that 





and. 


Fi{a;bi,b 2 ;c] zi,0) = 2 i^i(a, ^ 2 ; c; l) 2 Fi(a, 6 i; c - 62 ; 2 ^ 1 ) 

Fi(a; 6 i, 62 ;c; 0 , 0 ) = 1 


(B.32) 


where Fi(a, b; c; z) is the Appell hypergeometrie funetion. Therefore, EquationfB .281 ean 
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be written as 


t{z,X) = 


^y27^{Rpfl + z)uoao{X)T^pln{z) 


X 2F1 


{I 1 _ A _L 
^2,-L ^ ^ 


6 + 



t^n^i(£-i) \ 



i^Ot^l(e-l) \ 

-ffo(“o-“i) y 


1 a;ot<Ji(e—1) . oJq 

2 IIo(uJo—uJi) ’ 


|B.1.4[ 3 T(z) is a third-order polynomial 


(B.33) 


In this case, the temperature profile is a third-order polynomial with three zeros, one posi¬ 
tive and two negative—to be physically plausible. In addition, one of the zeros is particu¬ 
larly small (|a; 2 | z). 


Hence, 


and. 




fri^o) = 

— UJq +U1QU11+U10U12 —UJi_U12 

U)QUJ\U)2 

f^iuji) = 

— OJ^+UllOJo+CJlOJ2—‘^0‘^2 

LJ0UJliJ2 

f'rp{U}2) = 

— OJ^+OJ20Jl+OJ20JO—‘^l‘^0 

LJ0UJiiJ2 


(B.34) 


(B.35) 


t(z,X) 


2n(z)ao(X)T^Po 




2(^Rp^Q “ 1 “ z^tUi 


FiiOi) 

dx, 


(B.36) 


where 


n{z) 



(B.37) 


It can be shown that 

xF,(F + ii-A-B;F + |;§i.3i). 


We can use Equation|B.38|for the present simplified case as we assumed that \u: 2 \ 2 ;. 
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Therefore, Equation|B.36|ean be written as 


_ 2n{z)ao{X)T^pl 2 F{u}2)+1 / 1 

2F(i^2)+1 °° ^ ‘^2^ r(F{oJo)+F{ui2)+^) 

X2F1 (F(a; 2 ) + i -FM F{uo) + F^u^) + §; |) • 


B.2 Transmission Spectrum Equations 


The transmission speetrum (Equation |3.5| ) ean now be rewritten as 

(1 + y) (^1 - dy, (B.40) 

where y = z/Rp^ and ce' = Ui/Rpfl. 

We find no solution to the general integral introduced in Equation |B.40| using standard 
mathematical functions. However, we find that the above integral can be expressed using 
standard mathematical functions if /r(-) has less than two zeros. If /t(-) has no zero, the at¬ 
mosphere is isothermal. The isothermal case is discussed in detail in Section [3.2.1[ If /r(-) 
has a zero, the temperature is a linear function of the altitude. We introduce the derivation 
for such a case below. Note that the case of a linear temperature profile encompasses the 
case of a dry-adiabatic atmosphere in hydrostatic equilibrium addressed in Section [B3| 

If /r(.) has a zero, the integral introduced in Equation|B.40| takes the form 



c 

2 


(1 + 1/) 1 - e 




dy, where 


F 


r(0,A) 
5-1 + 


(B.41) 

£—1 

Hq f!p (ujq) 


The integral introduced in Equation|B.41|can be rewritten using the following identity: 


/(l + 2/)(l-e D)®)d|/= |(|/ + 2) 

(B.42) 

+ f (|,C'(l - + constant, 

because C, D, E, and y G in the context of transmission spectroscopy. Therefore 
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Equation|B .41 1 can be rewritten as 


c = oJq + 2uJq 

-2‘^Af (r(i)-r(i,A,)) 
+ 2 fAf (r (A) - r (A A,)). 


We recall that 


(B.43) 


r(a;) — r(a:, a) = * dt 

= 7(x,a), 

the lower incomplete gamma function. In addition, 

a""^r(a;,a) = dt 

El — x(^0j)y 


(BAA) 


(B.45) 


the generalized exponential integral. For the present derivation, we introduce ei_a;(a), the 
“lower” generalized exponential integral, which we define as 

rl ^-at 

ei_^(a) = J dt. (B.46) 

Using the lower generalized exponential integral, we can rewrite Equation|B.43|as 


c = UJ. 


a 


+ 2UJn 


2^ + CJn 




/2 


F 






(B.47) 

(B.48) 


because Wq 1 and ^4^ > 5—i.e., optically-active spectral bands, see the definition of 


Rpfi. Therefore, we obtain from Equation 3.16 


2/e//(^) 




(B.49) 


because (^F — ei_|^(y4A)j 1. 
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Figure B-1: Convergence of Ti^jeff) to e for non-isothermal atmosphere. 


As a result, we provide an analytical form for the effective planetary radius (Equa- 
tion |3.3[ ): 


R,(\) ~ Bp,o+y (f-e,_.(A)), 




r : r(r, A) = Aa 





(B.50) 

(B.51) 


We note that as expected 


lim Aa 

_F—>00 



F 




(B.52) 


Indeed, as F —cx) the temperature profile tends to an isothermal. Hence, one retrieves 
the same formulation as the one found in Section [3.2.1| for Teq (see Equation 3.21). We 
emphasize this convergence in Figure[BA] —note that for Aa > 5, T(yeff) is independent 
of Aa. 
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B.3 The Case of Dry Adiabatic Lapse Rate 


We focus in this Section on the formulation for atmospheric layers where the tempera¬ 
ture profile is solely affected by an adiabatic lapse rate. An adiabatic lapse rate refers 
to the negative/positive change in temperature experienced by a parcel of air as it moves 
upwards/downwards by convection and expands/shrinks without exchanging heat with its 
surroundings. The balance between potential energy and kinetic energy can be written as 


1 


Cp dT 


dp, 


(B.53) 


P 


where Cp is the specific heat at constant pressure and p is the gas density. In the context of 
hydrostatic equilibrium (dp = —pg dz), Equation|B.53|can be rewritten as 



(B.54) 


where =Tlji is the (dry) adiabatic lapse rate. 

Cp 

Equations |B.3[|B.lTj and |B.14| can thus be rewritten for the case of an adiabatic atmo¬ 
sphere in hydrostatic equilibrium. Eirst, Equation|B.3|is now 



(B.55) 


Therefore, 



(B.56) 


(B.57) 


Hence, we obtain for the pressure and number-density profiles 



(B.58) 


(B.59) 
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where Equation |B.58| is the barometrie formula for a non-zero lapse rate. 

Finally, the slant-path optieal depth and the transmission speetrum equation ean be 
written as 


t(z,A) = r(0,A)(l 

Rp{A) ~ Rpfi + 


r LR 
-1-0 
ToHo 


LR. 


^O^LR 


Tn - H.T 


QJ- LR 


Tn 


, and 
Ho^lr 


H.T 


To 


or LR 


To-HqTpjI 


(r(0,A)) 


(B.60) 

(B.61) 


where 


I T r( 

r(0,A) ~ J2TT{Rpfi +z)^aoiX)no:^ 

V i- LR i ( o 


To 


HoLlr' 


r(i + 


To 


HoTlr 


y 


(B.62) 


assuming the extinetion proeess to be Rayleigh seattering, i.e., = 0, e = 0} (Equations |B.15| 

and B.41). Using the sealing law funetion for the Rayleigh-seattering eross seetion, cTsc(A) = 
<to(A/Ao)“, we show in Figure B-2 the effeet of a and on Rp{X). Figure; 


B-2 


shows that 


Rp{X) does not present a Rayleigh-seattering linear slope—see Seetion [33| and Figure[4^ 
[6] —but rather a slope that gets steeper towards larger wavelength as these relate to larger 
atmospherie seale height—i.e., deeper and hotter atmospherie layers. 



Figure B-2: Rayleigh “slope” for non-isothermal atmospheres. 
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Appendix C 


Simple Validation Tests for Transmission 
Spectrum Models 


We introduced here how the analytical derivations of Section [J!2| and Appendix]^ can be 
used to validate transmission spectrum models. In particular, the transmission spectrum of 
a synthetic atmosphere in hydrostatic equilibrium with a dominant scatterer has a simple 
analytical formulation. Hence, it is suitable case for transmission model validation. 

(i) If for such an atmosphere T{z) = Tq, the slant-path optical depth (Equation |3.4| ) can 
be formulated as 


t{z,\) ~ (T(A)noe ^27r{Rpfl + z)H, (C.l) 

see Equation |3.14[ and the effective planetary radius as 

Rp{\) ~ i?p,o + H (^-fEM + In (noa{X)y/2nRpfiH'^'^ , (C.2) 

see Equation |3.12| if noao{X) ^y27^RpfiH 3> 1—for validation tests, set nocr(A) > 

10, i.e., ensure the atmospheric layer at Rp Q to be optically thick. In other words, by setting 
up the profile of the atmospheric absorption coefficient to 

a{r) = (C.3) 
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where ao,r, and H are respeetively the absorption eoeffieient at the referenee radius (Rp^, 
e.g., the synthetie planet’s surfaee), the radial distanee, and the atmospherie seale height, 
one should be modeling a transmission speetrum of the form 


AF 

~Y~ 


Rpfi + H (jem + In {ao^y27iRpfiH)) 


R. 


(C.4) 


(ii) If for sueh an atmosphere T(2;) oc z—i.e., T{z) = To(l — ^), the slant-path optieal 
depth (Equation |3.4[ ) ean be formulated as 

~ J27r(i?p,o + ^)a;oao(A)no ^ il- — ] \ (C.5) 

U2 + V ^0/ 

see Equation[B.25[ and the effeetive planetary radius as 
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{r{0,X)) 

Pi' ^0 

^ Ho' 


rfi + 

2 ^ Ho' 


(C.6) 

(C.7) 


see Equation [B.50[ if r(0, A) S> 1—for validation tests, set r(0, A) > 10, i.e., ensure the 
atmospherie layer at Rp^ to be optieally thiek. In other words, by setting up the profile of 
the atmospherie absorption eoeffieient to 


a{z) 


tto 



Uo 


- 1 + 


^0 


(C.8) 


where ao, 2 ;, and Hq are respeetively the absorption eoeffieient at the referenee radius (Rp^, 
e.g., the synthetie planet’s surfaee), the altitude, and the atmospherie seale height at Rp^, 
one should be modeling a transmission speetrum of the form 


AF 

~Y~ 


Rp,o + 


i^oHo 

‘^o — Ho 
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— e_Eo_ 
‘^Q-Ho 
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R. 






(C.9) 


178 

























Appendix D 


Modeling Transmission Spectra by 
Solving Rp{X) = {r : T(r, A) = 


In Section [3^ and Appendix]^ we show the pivotal role of the Euler-Maseheroni eon- 
stant (7EM, |EuIer||1740[ ) in transmission speetroseopy—see Equations |3.43[|B.5T] and 
Here, we introduee a computationally efficient way to model transmission spectra. Solving 


B.52 


Rp{X) = {r:r(r, A) = e '’'®"},or (D.l) 

heffiX) = {z : t{z, A) = , (D.2) 


is computationally more efficient than a direct numerical integration of Equation |3.3| as 
it solely requires solving Equation |D.l| or |D .21 knowing the dependency of each species’ 
extinction cross-section on the pressure and the temperature parameters—e.g., A^, a^, and 
(Equation|3.34[). These parameters are described by quantum physics or they can be 


measured in the lab (e.g., Rothman et al. 20091. This is why we advocate in this thesis that 
accurate extinction cross section databases covering various atmospheric conditions (i.e., 
temperature, pressure, composition) will be essential in the future to interpret consistently 
exoplanet transmission (and emission) spectra (Section[3.4.1|). 
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Error on [H] 

Figure D-1: Distribution of the error on heff{X) (expressed in scale heights) when esti¬ 
mated from heff{X) = {z : t{z, X) = for an Earth-sized planet with a isothermal 

and isocompositional atmosphere with the same abundances as at Earth’s surface. The er¬ 
ror is below 3% for 99.7% of the lines. An error below 3% corresponds to an error on the 
simulated effective height below 250 meters. 


Validation with Isothermal and Isocompositional Atmospheres 

We show here that it is adequate to model a planet’s transmission spectrum via Rp{X) = 
{r : r(r, A) = with an isothermal and isotermal atmosphere. In particular, we 

model here an Earth-sized planet with a isothermal and isocompositional atmosphere with 
the same abundances as at Earth’s surface. Eigure [D^ shows the resulting error on h^ffiX) 
is below 3% of the scale height for 99.7% of the active spectral bins. An error on heff{X) 
below 3% corresponds to an error of the simulated effective height below 250 meters— 
which is sufficient to advocate for the adequacy of modeling transmission spectra by solv¬ 
ing heff{X) = {z : t{z, X) = 


Validation with Earth’s Atmosphere 


We show in Eigure [D^ the error on heff{X) when modeling Earth’s transmission spectra 
using heff{X) = {z : t{z, X) = The error on heff{X) is below 18% for 99.7% of 
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Error on [H] 


Figure D-2: Distribution of the error on /ie//(A) (expressed in scale heights) when es¬ 
timated from heff{X) = {z : t{z,X) = for Earth—temperature-pressure-mixing 

ratio profiles from Cox and PilachowsEl(2000|. The error is below 18% for 99.7% of the 


lines. An error below 18% corresponds to an error on the simulated effective height is below 
1500 meters. The error is larger than for an isothermal-isocomposition Earth (Eigure [D^ 
mainly because water’s mixing ratio drops significantly in the troposphere (as highlight by 
the comparison of Eigs[3^and 3-7). 


the active lines. An error on heff{X) below 18% corresponds to an error of the simulated ef¬ 
fective height below 1500 meters—which is sufficient to advocate that modeling transmis¬ 
sion spectra based on solving heff{X) = {z : t{z, X) = is adequate. We note that 

while Earth’s case deviates significantly from the assumptions underlying our derivations 
(atmospheric temperature and mixing ratios strongly dependent on the altitude), modeling 
its spectrum using Tgq ~ remains adequate. The main origin for the error on h^ffiX) 

comes from the strong water maxing ratio—as shown in Eigures [3^ and |3-7[ 

In summary, modeling transmission spectra via Rp{X) = {r : r(r, A) = is ad¬ 

equate if the vertical scale for temperature/composition changes are larger than the local 
scale height, as emphasized by Equation B.52 in Appendix]^ 
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